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^ Abstract  

This  note  develops  a new  method  for  the  solution  of  EMP  in- 
er act ion  problems.  Basically  it  involves  expanding  the  solu- 
ion  in  terms  of  its  singularities  in  the  Laplace  transform  or 
omplex  frequency  (or  s)  plane.  In  the  time  domain  each  term 
omes  from  an  inverse  transform  of  the  corresponding  term  in  the 
singularity  expansion.  Finite  size  objects  with  well  behaved 
media  have  only  poles  in  the  finite  s plane  for  their  delta 
function  response.  These  factor  into  terms  involving  the  class- 
ical natural  frequencies  and  modes  but  in  addition  bring  out 
factors  which  we  call  coupling  coefficients  as  well  as  the  pos- 
sibility of  higher  order  poles  besides  simple  poles,  but  still 
of  finite  order  in  the  finite  s plane.  If  the  incident  waveform 
ir  has  singularities  in  the  finite  s plane  the  response  can  be  gen- 
x erally  split  into  an  object  part  (containing  object  po  ** 5 and  a 
waveform  part  containing  the  waveform  singularities.  The  object 
j poles  directly  give  amplitudes,  frequencies,  damping  constants, 

J and  phases  for  the  damped  sinusoidal  waveforms  seen  so  commonly 
Z in  EMP  tests  using  pulsed  waveforms.  There  is  some  latitude  in 
the  calculation  of  coupling  coefficients  and  some  difficulties 
^ are  discussed.  ___  j- — j 
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Foreword 


Approximately  this  last  summer  I started  looking  at  this 
singularity  expansion  method , Based  on  physical  arguments  the 
general  form  of  the  extension  and  how  the  terms  factored  for 
finite  objects  to  exhibit  dependence  on  separate  variables  of 
the  BMP  interaction  problem  soon  became  apparent  to  me.  In  par- 
ticular some  physical  observations  generalized  from  many  EMP 
tests  in  time  domain  led  to  the  pole  expansion  concept  in  fre- 
quency domain.  This  forms  the  starting  point  for  what  follows 
in  this  note. 

Next  chronologically  I started  some  discussion  going  on 
this  subject,  in  particular  with  some  of  my  colleagues  at  North- 
rop Corporate  Labs  in  Pasadena.  I certainly  wish  to  thank  them, 
in  particular  Dr,  L.  Marin,  Dr.  K-  S.  H.  Lee,  Dr.  R.  W.  Latham, 
and  Dr.  F.  Tesch  (who  is  now  with  Dikewood)  for  many  very  stim- 
ulating conversations  about  this  technique.  They  certainly 
helped  me  test  and  refine  some  of  the  concepts.  They  are  also 
presently  working  on  reports  to  further  refine  the  method  and 
calculate  some  example  problems.  In  particular  Drs.  Marin  and 
Latham  have  rather  far  advanced  some  analytic  solutions  of  the 
magnetic-field  integral  equation  for  finite  size  perfectly  con- 
ducting objects  in  terms  of  the  singularity  expansion. 

In  September  there  was  a meeting  in  Pasadena  with  some  sig- 
nificant attention  given  to  this  subject.  I would  like  to  thank 
everyone  who  came  to  that  meeting  for  the  stimulating  discussion 
on  this  subject.  On  various  occasions  both  at  this  meeting  and 
on  other  occasions  I have  had  occasion  to  discuss  this  matter 
with  various  people.  In  particular  I would  like  to  thank  Prof. 
R.  J.  Garbacz  of  Ohio  State  U.,  Prof.  C.  T.  Tai  of  U.  of  Michi- 
gan, Prof.  S.  W.  Lee  of  U.  of  Illinois,  Prof.  C.  Taylor  of  U.  of 
Miss.,  and  Dr.  A.  Poggio  of  Cornell  Aeronautical  Lab.  Some  of 
these  people  are  already  beginning  studies  on  various  aspects 
and  specific  problems  concerned  with  the  singularity  expansion 
method . 


...when  suddenly  a white  rabbit  with  pink  eyes  ran  close  by 

her. 

There  was  nothing  so  remarkable  in  that;  nor  did  Alice 
think  so  very  much  out  of  the  way  to  hear  the  Rabbit  say  to  it- 
self, "Oh,  deer!  Oh,  dear!  I shall  be  too  late!"  (when  she 
thought  it  over  afterward  it  occurred  to  her  that  she  ought  to 
have  wondered  at  this,  but  at  the  time  it  all  seemed  quite  nat- 
ural) ; but  when  the  Rabbit  actually  took  a watch  out  of  its 
waistcoat  pocket,  and  locked  at  it,  and  then  hurried  on,  Alice 
started  to  her  feet,  for  it  flashed  aoroes  her  mind  that  she  had 
never  before  seen  a rabbit  with  either  a waistcoat  pocket  or  a 
watch  to  take  out  of  it,  and,  burning  With  curiosity,  she  ran 


across  the  field  after  it,  and  was  just  in  time  to  see  it  pop 
down  a large  rabbit  hole  under  the  hedge. 

In  another  moment  down  went  Alice  after  it,  never  once  con- 
sidering how  in  the  world  she  was  to  get  out  again. 

(Lewis  Carroll,  Alice  in  Wonderland) 
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I.  Introduction 


(•  * 

This  note  takes  off  on  a new  course  in  the  treatment  of  the 
interaction  of  electromagnetic  fields  with  bodies  located  in 
free  space  or  in  other  simple  media,  including  the  effects  of 
the  proximity  of  one  body  with  respect  to  mother  (such  as  a 
body  in  an  EMP  simulator) . 

By  way  of  introduction  some  physical  observations  are  in 
order.  Suppose  one  excites  an  object  such  as  a missile,  an  air- 
craft, a building  with  various  conductor  geometries,  etc.  with  a 
fast  electromagnetic  pulse.  What  are  the  general  characteris- 
tics of  the  resulting  waveforms  for  fhe  various  electromagnetic 
quantities  (such  as  current,  charge,  etc.)  associated  with  this 
object?  Someone  with  much  experience  in  EMP  testing  of  such  ob- 
jects in  EMP  simulators  could  observe  that  an  extremely  common 
characteristic  of  such  waveforms  is  the  presence  of  one  or  more 
exponentially  damped  sinusoidal  oscillations.  This  is  the  case 
not  only  in  the  excitation  of  the  internal  circuitry  but  also 
for  the  body  geometry  as  well.  This  electromagnetic  resonance 
phenomenon  is  particularly  pronounced  for  long  and  comparatively 
slender  conductors  such  as  the  main  body  resonance  of  a missile 
or  the  body  and  wing  resonances  of  various  aircraft.  These 
damped  sinusoids  in  the  response  are  observed  not  only  in  EMP 
interaction  studies  (both  experimental  and  theoretical)  but  also 
in  fast  pulse  time-domain- type  radar  scattering  studies. 

Not  all  aspects  of  the  electromagnetic  response  of  objects  i 

look  like  damped  sinusoids.  Parts  of  the  time  domain  response 
may  look  like  the  incident  waveform,  or  perhaps  its  time  deriva- 
tive or  time  integral.  It  would  also  appear  that  in  some  cases 
even  more  complex  types  of  responses  occur. 

Looking  at  the  forms  of  the  observed  responses  one  might 
ask  if  there  is  some  way  that  these  observable  features  of  the 
waveform  can  be  found  directly  when  one  calculates  the  object 
electromagnetic  response.  Can  the  amplitude,  frequency,  damp- 
ing constant,  phase,  etc.  of  each  damped  sinusoid  be  directly 
calculated?  How  do  they  depend  on  the  incident  wave?  Can  other 
kinds  of  response  which  give  waveforms  related  to  the  incident 
waveform  be  directly  calculated?  ( 

The  purpose  of  this  note  is  to  begin  to  characterize  elec- 
tromagnetic interaction  with  objects  in  terms  of  quantities  di- 
rectly identifiable  with  various  characteristics  of  resulting 
interaction  waveforms.  Some  characteristics  are  associated  with 
the  object  characteristics  including  the  presence  of  neighboring 
objects.  Other  characteristics  are  associated  with  the  waveform 
of  the  incident  field.  Yet  others  are  associated  with  the  spa- 
tial distribution  of  the  incident  fields,  such  as  specified  by 
direction  of  incidence  and  polarization.  What  is  in  effect  ac- 
complished here  is  a decomposition  of  the  interaction  problem 

/’ 
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into  various  quantities  which  depend  on  different  variables  of 
the  problem.  The  dependence  of  the  interaction  on  different 
variables  can  then  be  separately  considered  resulting  in  a con- 
siderable simplification  in  understanding  how  the  resulting 
electromagnetic  interaction  can  vary  over  all  possible  varia- 
tions of  the  parameters  of  a particular  problem  being  considered. 
This  effectively  extends  the  complexity  of  the  object  geometries 
one  may  be  willing  to  consider  for  detailed  calculations. 

Having  identified  what  appear  to  be  exponentially  decaying 
sinusoids  in  typical  interaction  experimental  data  one  might  use 
thiB  as  a clue  toward  finding  a mathematical  representation  of 
the  electromagnetic  interaction  which  has  these  terms  as  part  of 
the  decomposition.  Consider  the  Laplace  transform  of  the  var- 
ious waveforms.  (The  two-sided  Laplace  transform  is  used 
throughout  the  note.)  The  Laplace  transform  of  an  exponentially 
damped  sinusoid  gives  a pair  of  complex  conjugate  poles  in  the 
complex  s plane  of  the  form  l/(s  - s^)  and  l/(s  - Sajj)  where  s 
is  the  Laplace  transform  variable,  where  sai  ° Sct2  with  the  bar 
above  a quantity  indicates  complex  conjugate,  and  where  and 
02  are  sets  of  indices  to  label  the  pedes  being  considered.  If 
one  could  find  these  poles  with  their  coefficients  in  the  com- 
plex s plane  from  an  integral  equation  or  other  form  of  the  so- 
lution (such  as  from  an  eigenfunction  expansion  of  the  solution) 
then  not  only  would  he  have  a representation  of  part  of  the  fre- 
quency or  Laplace  domain  solution,  but  also  of  a part  of  the 
time  domain  waveform  (damped  sinusoids)  as  well.  Such  poles  in 
the  complex  s plane  are  termed  natural  frequencies  of  the  object 
since  they  are  frequencies  for  which  the  object  can  have  a re- 
sponse (in  Laplace  or  complex  frequency  domain)  with  no  excita- 
tion in  the  form  of  an  incident  wave.  If  the  body  is  excited  at 
a natural  frequency  then  its  response  is  infinite  at  that  com- 
plex frequency. 

Suppose  one  were  to  take  a solution  for  some  interaction 
problem  expressed  in  the  Laplace  domain  either  explicitly  or  im- 
plicitly (such  as  in  the  form  of  an  integral  equation) . Fur^ 
thermore  suppose  one  wishes  to  convert  this  into  a time  domain 
solution.  This  can  be  done  using  the  inverse  Laplace  transform 
integral,  a contour  integral  in  the  complex  plane.  This  contour 
can  be  deformed  in  the  complex  s plane,  passing  it  through  re- 
gions where  the  response  is  an  analytic  function  of  s.  On 
reaching  singularities  such  as  poles  and  branch  points  the  con- 
tour can  be  deformed  around  the  poles  and  branch  cuts  to  obtain 
terms  associated  with  each  separate  pole  and  branch  cut  and  per- 
haps a contribution  from  a portion  of  the  contour  for  |s|  *►  «. 

Our  basic  approach  then  in  expressing  the  solution  to  electro- 
magnetic interaction  problems  is  to  express  it  as  a sum  of  such 
terms  in  both  Laplace  (frequency)  domain  and  time  doatain.  Dif- 
ferent types  of  terms  will  have  different  properties  and  we  wish 
to  understand  these  properties  in  detail  so  as  to  take  advantage 
of  the  decomposition  of  the  interaction  problem  into  its  various 
parts. 
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Not  only  are  the  natural  frequencies  of  the  object  of  con- 
cern because  there  can  be  other  singularities  in  the  response. 
Associated  with  certain  kinds  of  object  geometries  one  may  also 
need  to  consider  branch  cuts  in  the  complex  s plane.  An  example 
of  such  a body  is  an  infinite  length  perfectly  conducting  cir- 
cular cylinder  for  which  the  branch  cuts  can  be  associated  with 
cylindrical  Hankel  functions.  This  note  is  mostly  concerned 
with  cases  that  the  object  response  does  not  have  branch  cut 
contributions  associated  with  the  object  characteristics.  How- 
ever the  same  approach  as  for  the  case  of  poles  may  be  used  by 
including  terms  for  the  branch  cut  integrals  in  the  general  ex- 
pansion for  the  inverse  Laplace  transform  integral.  Back  in 
Laplace  domain  the  individual  branch  cut  terms  can  be  found  by 
first  subtracting  all  the  pole  contributions  and  then  treating 
what  is  left.  Since  the  only  important  contribution  at  the 
branch  cut  is  the  change  in  the  function  across  it,  then  in  cal- 
culating the  branch  cut  integral  one  could  use  this  change  to 
define  an  appropriate  tern  only  associated  with  this  change 
(along  the  entire  branch  cut) . Such  problems  are  not  considered 
in  this  note  but  it  is  determined  that  for  a class  of  objects  of 
interest  there  are  no  branch  cut  contributions  associated  with 
the  object  geometry. 

The  incident  waveform  will  also  typically  have  singulari- 
ties in  its  Laplace  transform  and  there  are  terms  in  the  object 
response  which  correspond  to  these.  The  object  response  can 
then  be  split  to  some  extent  into  terms  associated  with  the 
geometry  and  other  electromagnetic  characteristics  of  the  object 
on  the  one  hand,  and  the  incident  waveform  characteristics  (in- 
cluding how  they  couple  to  the  object)  on  the  other  hand.  In 
many  cases  the  waveform  contributions  will  be  through  rather 
simple  singularities  such  as  simple  poles  associated  with  expo- 
nential waveforms . 

Having  decomposed  the  object  response  into  various  terms  in 
this  manner,  one  can  then  see  how  accurately  the  first  so  many 
terms  describe  the  complete  object  response  in  various  of  its 
characteristics  in  both  frequency  and  time  domains.  It  is  not 
apparent  that  all  interaction  problems  can  be  most  conveniently 
characterized  using  the  natural  frequencies  and  modes  and  other 
singularities.  However,  for  highly  resonant  structures  it  ap- 
pears to  considerably  simplify  the  comprehersion  of  the  impor- 
tant features  of  electromagnetic  interaction  with  the  structure. 
Other  techniques  will  continue  to  be  valuable  and  waveform  and 
frequency-response  detailed  calculations  will  still  be  needed. 

In  sotr\e  cases  the  waveforms,  for  example,  will  be  useful  in  de- 
termining how  many  natural  modes  etc.  are  needed  for  the  object 
and  over  what  range  of  the  object  parameters  the  incomplete  sum 
of  modes  is  adequate. 

This  note  first  considers  the  general  form  of  these  solu- 
tions which  separates  out  various  aspects  of  the  object  response 
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Then  based  on  finite  matrix  formulation  of  integral  equations  in 
a form  such  as  used  for  numerical  solutions  various  properties 
of  the  solution  can  be  found  based  on  the  analyticity  properties 
of  the  matrix  elements  and  vector  components.  For  instance,  the 
conclusion  that  only  poles  of  finite  order  appear  in  the  object 
response  for  finite  size  objects  applies  not  only  to  perfectly 
conducting  objects  but  ones  of  finite  conductivity  as  well. 

Numerous  topics  are  then  briefly  discussed  to  point  out  many 
areas  for  further  investigation.  Finally  an  appendix  discusses 
the  special  "natural  frequency"  s « 0 and  another  appendix  works 
out  the  object  response  for  a perfectly  conducting  sphere  as  an 
illustration  of  natural  frequency,  mode,  and  coupling  coeffici- 
ent calculation  and  indexing. 

In  addition  to  what  is  discussed  in  this  note  many  refine- 
ments of  the  singularity  expansion  method  are  possible  and  var  J 

ious  extensions  of  the  results  vrould  seem  possible.  In  solving  I 

specific  boundary  value  problems  with  this  approach  some  other 
general  results  may  be  suggested  by  the  data,  thereby  focusing 
attention  on  the  proof  or  disproof  of  these  conjectures  and  the 
consideration  of  other  boundary  value  problems  which  betuer  ex- 
hibit the  same  general  results  or  test  their  validity. 


While  this  note  considers  the  solution  of  electromagnetic 
boundary  value  problems  in  terms  of  natural  frequencies  and 
other  singularities  the  technique  can  be  applied  to  experimental 
data  as  well.  From  frequency  or  time  domain  data  one  should  be 
able  to  determine  natural  frequencies  a^d  modes  and  other  sin- 
gularity characteristics  by  extending  the  data  to  the  remainder 
of  the  complex  s plane  using  Laplace  transform  techniques  or 
even  analytic  continuation.  Of  course  there  are  numerical  error 
problems  as  in  other  data-  reduction  processes  and  this  will  also 
require  quantitative  understanding. 

Since  this  is  a new  approach  to  the  solution  of  electromag- 
netic boundary  value  problems  there  is  clearly  much  work  to  be 
done  to  fully  understand  its  many  ramifications. 
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II . The  Form  and  Some  General  Characteristics  of  the  Singular- 
ity Expansion  for  Some  of  tKe  Simpler  Cases 

Let  uo  now  look  at  some  of  the  advantages  nssoci  cod  with 
the  singularity  expansion  method  because  of  the  form  of  the  ex- 
pansion in  simple  cases. 


A 


First  we  write  the  complex  Laplace  transform  variable  in 
terms  of  its  real  and  imaginary  parts  for  notational  purposes  as 


s = ft  + iw  (2.1) 

An  arbitrary  function  of  time  f(t)  which  could  be  in  general  a 
scalar,  vector,  or  tensor  of  arbitrary  rank  has  a Laplace  trans- 
form (bilateral  or  two  sided)  assuming  the  integral  of  |f(t)| 
over  any  finite  interval  exists  defined  byi 


CO 

L[f(t)J  = f(s)  = f(t)e"stdt 

•'—CO 


(2.2) 


s 

where  the  integration  is  taken  on  the  real  t axis.  This  is  the 
two  sided  Laplace  transform  (indicated  by  a tilde  - over  the 
function)  where  f(t)  is  required  to  have  a behavior  such  that 
f(s)  exists  and  is  analytic  in  some  strip  ft_  < Re[s)  < ft+  in  the 
s plane.  The  inverse  transform  (where  f(t)  is  of  bounded  varia- 
tion nearr  t)  is  given  by 

-ft  +i» 

L*-1  If  (s)  ] = f(t)  = / ° f(s)estds  (2.3) 

*wi  'ft  -i» 
o 

where  the  limits  to  +i~  can  be  interpreted  in  a Cauchy  principal 
value  sense  and  where  ft_  < ft0  < ft+  unless  f (t)  is  discontinuous 
at  t in  which  case  the  inversion  gives  [f(t-)  + f(t+))/2.  In 
our  cases  of  interest  f(t)  «=  0 for  t < t0  in  which  case  ft+  *»  *» 
and  the  transform  effectively  reverts  to  a one  sided  Laplace 
transform  with  the  lower  limit  as  to.  Typically  also  ft-  » 0 as 
long  as  f (t)  does  not  grow  as  fast  gs  an  increasing  exponential 
for  t *►  +«.  Thus  we  normally  have  f(s)  an  analytic  function  of 
s for  Re[s]  > 0,  the  right  half  of  the  s plane,  and  the  inver- 
sion integral  is  defined  along  Re[s)  = ftQ  > 0. 

The  essence  of  the  singularity  expansion  method  involves 
evaluating  f(s)  (which  may  be  surface  current  density  or  various 
other  electromagnetic  quantities)  by  evaluating  f(s)  in  terms  of 
the  left  half  plane  singularities  (Re[s3  < 0) . Express  the  time 
domain  form  f(t)  in  terms  of  these  same  sXngularities  as  would 
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be  done  by  deforming  the  contour  for  the  inversion  integral 
(equation  2.3)  into  the  left  half  plane  and  splitting  the  inte- 
gral int^  parts  associated  with  each  singularity.  Note  auto- 
matically that  since  we  are  only  concerned  with  f(t)  real  for 
real  t then  the  Laplace  transformed  f(t)  has  some  symmetry  which 
can  be  found  by  splitting  the  transform  integral  into  real  and 
imaginary  parts.  Denoting  the  complex  conjugate  by  a bar  - over 
the  quantity  we  have 

f(s)  - f~( s)  (2.4) 

Singularities  are  then  automatically  symmetrically  placed  with 
respect  to  the  R axis  except  that  branch  cuts  can  be  moved 
around  as  long  as  the  branch  points  stay  symmetrically  placed 
with  respect  to  the  R axis.  For  convenience  we  constrain  the 
branch  cuts  to  also  be  symmetrically  placed  with  respect  to  the 
R axis  so  that  equation  2.4  always  holds  except  of  course  right 
at  the  singularities.  Having  found  the  term  associated  with  one 
singularity  we  then  automatically  have  the  result  for  the  conju- 
gate singularity.  Of  course  this  does  not  help  us  for  those 
singularities  on  the  R axis.  When  we  index  the  terms  associated 
with  each  singularity  with  a set  of  integers  we  can  also  adopt 
the  convention  of  a sign  reversal  on  one  of  the  integer  indices 
corresponding  to  conjugate  positioned  singularities.  For  such 
an  index  positive  integers  can  be  associated  by  convention  with 
o>  > 0 (upper  half  plane)  a zero  integer  for  w = 0 (the  R axis) 
and  a negative  integer  for  w < 0 (the  lower  half  plane) . Alter- 
nately another  symbol  can  be  introduced  to  indicate  which  of  a 
conjugate  pair  is  meant. 

While  there  are  various  forms  of  incident  waves  that  one 
might  use  we  choose  the  commonly  used  plane  wave  for  our  exam- 
ples. If  one  wishes,  more  complex  field  distribution  can  be 
found  by  superposition  of  many  plane  waves.2  For  our  present 
purposes  wo  consider  an  incident  plane  wave  propagating  in  the 
ex  direction  (independent  o£  s)  wi£h  electric  field  polarization 
in  some  combination  of  the  62  and  03  directions.  Here  e is  used 
for  a unit  vector.  The  three  unit  vectors  for  our  plane  wave 
are  all  mutually  orthogonal  and  form  a right  handed  system  of 
unit  vectors  as 

®1  X ®2  * ®3  (2.5) 

As  shown  in  figure  2.1  this  plane  wave  ^s  incident  on  some  ob- 
ject of  finite  linear  dimensions.  Let  r denote  the  observer  po- 
sition and  r*  coordinates  on  the  object.  Then  by  ap  object  of 
finite  dimensions  we  can  require  jr'|  < rQ  for  all  r'  where  rp 
is  some  finite  radius.  (Note  that  all  dimensions  are  rational- 
ised MKSA  throughout  the  note.)  Typically  the  coordinate  origin 
(r  * 5)  would  be  chosen  near  or  even  inside  the  object  If  we 
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Tt%  AND  5^  ARE  ALL 

MUTUALLY  ORTHOGONAL  UNIT  VECTORS 

WITH  FIXED  DIRECTIONS. 

E?  , H^nc,  AND  ^ ARE  ALL 
MUTUALLY  ORTHOGONAL  AT 
EACH  TIME  AND  POSITION  IN  SPACE. 


REFERENCE  AXIS 
(PERPENDICULAR  TO  ip 
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FIGURE  2.1  PLANE  WAVE  INCIDENT  ON  OBJECT 
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Refine  some  reference  axis  (as  shown  in  figure  2.1)  then  given 
ei  we+can  choose  e2  parallel  to  £he  plane  which  is  parallel  to 
both  ex  and  the  reference  axis;  e3  will  be  perpendicular  to  this 
plane.  If  the  object  has  an  axis  of  symmetry  this  would  nor- 
mally be  chosen  as  the  reference  axis.  In  spherical  coordinates 
as  in  one  the  appendices  the  z axis  (or  9 » 0,  ir)  would  be 
chosen  as  the  reference  axis  for  convenience. 

Our  incident  plane  wave  is  assumed  to  propagate  in  free 
space  with  a propagation  constant 


Y = ik  - | 


(2.6) 


where  the  speed  of  light  in  vacuum  is 


(2.7) 


The  wave  impedance  of  free  space  is 


(2.8) 


The  permeability  of  free  space  is  y0  and  the  permittivity  of 
free  space  is  eQ.  It  is  not.  strictly  necessary  for  the  medium 
to  be  free  space;  and  eQ  can  *>e  regarded  as  parameters  of  the 
large  volume  (ideally  infinite  in  size)  of  the  medium  in  which 
the  object  of  interest  is  placed.  However  for  some  results  it 
may  be  necessary  that  this  medium  have  zero  conductivity  (be 
lossless)  so  zero  conductivity  is  specified  for  the  infinite 
medium  for  considerations  in  this  note. 

The  general  form  of  our  incident  plane  wave  can  now  be 
written  as3#4,5 

e3 

(2.9) 

e2 


where  E0  is  a scale  factor  with  dimensions  volts/m.  The  two  in 
cident  fields  are  related  by 


* [j('-  - ~5~} 

• M,(.  - VK  • -.(•  • W 
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(2.10) 


Sine'*'0  - ^ h x 


iinc<?'t)  " ‘ZoSi  x 

In  Laplace  form  the  incident  plane  wave  is 


Einc^'s)  " Eo**2*s*e2  + ^3  (®)  ®2^ e 


-Yel*r 


5inc(*'s)  " Z2^2(s)®3  " *Vs)®2Ie  Y 1 
o 


(2.11) 


Note  that  for  each  of  the  two  independent  polarizations  we  can 
have  separate  waveform  functions  f2(t*)  and  f3(t*)  where  the  re- 
tarded time  is  defined  by 


t*  = t - 


e.  *r 
__ 


(2.12) 


The  waveform  functions  have  subscripts  which  relate  them  to  the 
polarization  vector  for  the  electric  field.  For  the  magnetic 
field  2 and  3 are  interchanged  with  a sign  reversal  in  one  ca^e. 

No  evanescent  waves  are  allowed  for  our  present  purposes,  so  ex 
(even  when  associated  with  the  Laplace  form)  is  taken  as  a real 
unit  vector  although  for  some  purposes  it  need  not  be.  This 
plane  wave  can  be  expanded  in  terms  of  the  vector  wave  function 
for  various  other  coordinate  systems  such  as  spherical  coordi- 
nates* (which  is  used  in  appendix  B)  and  cylindrical  coordinates.5 

For  convenience  one  may  introduce  a unit  dyadic  plane  wave 
for  Laplace  domain  purposes  as 


(1  0 0\  -ye, *r 
0 1 Ole  A 
0 0 1/ 


(2.13) 


(with  bx,  02  indicating  in  this  case  a pair  of  indices  each 
ranging  from  1 to  3)  which  in  the  time  domain  is  a dyadic  delta 
function  plan^  wave  as 


f 


1 


m 


(* 


(2.14) 
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The  subscript  1 indicates  propagation  in  the  ex  direction.  Mul- 
tiplying (in  a dot  product  sense)  this  unit  dyadic  p?ane  wave  by 
a waveform  and  polarization  gives  our  incident  plane  wave  as 


^Lnc(?'s)  " EoI?2(s>®2  + *3(8)®31  ' *1 
^inc(?'8)  “ F^2(8)®3  " l3<8>®21  * *1 


(2.15) 


In  the  time  domain  convolution  is  also  needed  which  we  might  in- 
dicate by 


m Eolf2<t)S2  + f3<t,S3»* 


Sinc<*.«  * f2<f2<t>*3  - f3<t>®2)* 


Note  that  the^dyadic  plane  wave  "contains”  both  longitudinal 
waves  in  the  ex  part  and  transverse  waves  in  the  e2  and  03  parts. 
Only  the  transverse  waves  satisfy  Maxwell's  equations  in  source 
free  media  and  so  the  2 and  3 components  giye  the  most  general 
uniform  plane  wave  propagating  parallel  to  ex. 

The  first  problem  in  the  solution  of  our  electromagnetic 
interaction  problem  is  then  to  find  the  object  response  to  two 
incident  waves 


where  these  can  be  taken  as  either  electric  or  magnetic  fields 
because  of  the  way  they  interchange  with  one  another  (with  a 
sign  reversal  in  one  case) . Knowing  the  response  to  these  two 
waves  (taken  as  both  electric  or  both  magnetic)  then  £he  re- 
sponse to  a general  plane  wave  incident  in  direction  ex  can  be 
formed  in  Laplace  domain  simply  by  re introduction  of  the  wave- 
form functions  f 2 and  £3  and  other  scale  factors  just  as  they 
appear  in  thg  incident  plane  wave  as  above.  In  time  domain  the 
response  to  U2  and  U3  can  be  separately  convoluted  with  f2  and 
f 3 and  scale  factors  reintroduced  to  obtain  the  general  solution. 

What  then  is  the  object  response  to  3x  and  to  $2  taken  for 
convenience  as  the  electric  field  components  to  correspond  di- 
rectly to  pol'/jrization?  Here  all  the  waveform  characteristics 
are  factored  out  leaving  two  problems,  each  of  which  has  an 


* 1 


• 1 


(2.16) 
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"incident  field  function”  which  has  no  singularities  in  the  en- 
tire s plane.  These  are  termed  entire  functions;  they  are  ana- 
lytic for  all  finite  s although  in  thir  case  they  each  have  an 
essential  singularity  at  <».  These  are  the  simplest  plane  waves 
for  our  purposes. 

There  are  many  kinds  of  objects  which  one  might  considc*. 
They  can  be  of  finite  size  or  infinite  size  in  various  shapes. 
They  can  be  composed  of  various  media  arranged  in  various  dis- 
tributions. In  this  note  we  concentrate  our  attention  on  a cer- 
tain kind  of  such  objects  and  work  out  some  general  results  for 
this  class  of  objects.  Other  kinds  of  objects  can  also  be  con- 
sidered in  this  way  and  some  comments  are  made  regarding  such 
other  classes  of  objects.  Ag  shown  in  figure  2.1  we  consider 
finite  bodies  described  by  |r'|  < rQ  for  all  points  on  the  body. 
Furthermore  for  some  of  the  considerations  and  the  example  in 
appendix  B we  take  this  object  as  perfectly  conducting  which  re- 
duces our  considerations  to  the  body  surface  which  has  surface 
charge  and  current  densities  as  quantities  of  primary  interest 
for  describing  the  electromagnetic  interaction.  This  body  is 
not  necessarily  continuous  but  may  be  composed  of  several  sepa- 
rate parts. 

In  a later  section  of  this  note  it  is  demonstrated  that 
current  density  (either  surface  or  volume)  for  a finite  object 
has  its  delta  function  response  corresponding  to  any  singulari- 
ties in  the  finite  s plane  expressed  in  terms  of  poles  of  finite 
order.  We  identify  these  poles  as  sa  where  a is  some  index  set 
which  indicates  which  pole  is  meant.  This  result  is  of  funda- 
mental importance  for  our  general  representation  of  the  solution 
to  the  interaction  problem.  We  have  a series  with  (s  - so)~na 
for  n<x  *=  1,  2,  3,  ...  as  a factor  in  each  term.  What  we  need 
now  is  the  rest  of  the  expression  at  each  of  these  poles  to  com- 
plete the  expansion  of  the  two  delta  function  plane  wave  re- 
sponses in  Laplace  domain.  If  there  is  more  than  one  order  of 
pole  at  8 » sa  the  a index  set  can  have  a number  to  designate 
each  term  in  the  expansion  corresponding  to  each  order  pole; 
clearly  this  number  could  be  just  na. 

These  poles  sa  are  the  natural  frequencies  of  the  object. 

By  a natural  frequency  is  meant  a value  of  s for  which  the  ob- 
ject has  a response  without  an  incident  field  exciting  it.  The 
natural  frequencies  are  generally  in  the  left  half  of  the  s 
plane  (Re[s)  < 0)  because  of  energy  loss  (to  radiation  in  the 
case  of  a perfectly  conducting  body)  except  that  natural  fre- 
quencies can  be  on  the  iw  axis  (0  ■ 0)  with  first  order  poles 
for  lossless  situations  such  as  for  interior  cavity  nodes.  How- 
ever in  the  time  domain  such  simple  poles  on  the  iw  axis  corres- 
pond to  undamped  sinusoids  which,  if  they  can  be  excited  by  the 
incident  wave,  must  continue  to  radiate  power  to  infinity  indef- 
initely by  reciprocity.  Since  this  would  violate  conservation 
of  energy  the  residues  of  such  poles  must  be  sero  and  we  drop 
them  from  consideration. 
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If  ihe  object  responds  at  s » sa  with  no  incident  f^eld 
£hen  such  response  must  be  independent  of  polarization  (e2  and 
e3>  and  direction  of  incidence  (ei)  since  they  determine  no  in- 
cident fields  for  this  ideal  case.  This  leads  to  the  concept  of 
a natural  mode.  By  a natural  mode  we  mean  a current,  charge, 
field,  etc.  distribution  associated  with  a body  self  oscillation 
at  a complex  natural  frequency  sa.  We  designate  natural  modes 
by  $a(F)  (*')  *or  A vector  quantity  F and  va'p)  (r‘)  for  a scalar 
quantity  F.  A scalar  quantity  F could  be  a surface  charge  den- 
sity Os  and  a vector  quantity  F could  be  a surface  current  den- 
sity 5S  parallel  to  the  object  surface. 

There  is  of  course  the  question  of  the  uniqueness  of  the 
natural  modes.  Clearly  the  natural  modes  can  be  modified  by  a 
scale  factor,  but  this  is  just  a problem  of  an  appropriate  nor- 
malization of  the  natural  modes.  Depending  on  the  problem  of 
concern  there  may  be  various  appropriate  normalizations.  Ap- 
pendix B considers  the  case  of  a perfectly  conducting  sphere  for 
which  all  the  natural  frequencies  correspond  to  first  order 
poles  (na  * 1 c.ily)  . For  the  perfectly  conducting  sphere  the 
natural  modes  are  types  of  spherical  harmonics  and  here  we  use 
definitions  which  fit  naturally  with  common  usage. 

The  perfectly  conducting  sphere  is  an  interesting  example 
in  that  all  the  natural  modes,  coupling  coefficients,  and  natu- 
ral frequencies  can  be  more  readily  calculated.  Any  general  re- 
sults for  perfectly  conducting  finite  sized  objects  must  be  true 
for  the  perfectly  conducting  sphere.  Thus  the  sphere  (and 
other  shapes  such  as  prolate  and  oblate  spheriods  and  disks)  can 
be  used  to  test  general  results.  Moreover  they  can  be  used  to 
form  a basis  for  conjecture  for  new  general  results.  For  ex- 
ample the  perfectly  conducting  sphere  has  only  simple  poles  (n« 

» 1 only)  which  completely  describe  its  response  to  a delta 
function  plane  wave. 

Aside  from  the  scale  factor  there  is  another  problem  in  de- 
fining the  natural  modes.  This  is  the  possible  degeneracy  of 
the  mc^es.  As  can  be  seen  for  the  sphere  problem  for  example 
(appendix  B)  the  modes  can  be  degenerate  in  which  case  associ- 
ated with  each  natural  frequency  there  may  be  several  modes. 

This  problem  can  easily  be  handled  by  the  a index  set  to  desig- 
nate separate  independent  modes.  There  are  various  different 
ways  to  define  the  different  modes  and  what  is  needed  is  a con- 
venient set  of  modes  with  the  minimum  humber  of  modes  necessary 
to  span  the  space  of  all  possible  distributions  of  the  quantity 
of  interest  associated  with  the  particular  natural  frequency. 

Of  course  the  sphere  has  a high  degree  of  symmetry  and  one  would 
expect  degeneracy  of  the  natural  modes  associated  with  symmetry. 
Bodies  with  a symmetry  axis  will  also  have  a degeneracy  of  the 
form  cos(n$)  or  sin(n$)  as  a factor  in  the  natural  modes  (axcept 
for  n » 0)  giving  at  least  two  independent  natural  modes  (for 
n £ 1)  for  each  natural  frequency.  Note  that  it  is  also 
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possible  (ac  least  in  the  case  of  the  sphere)  to  define  the  nat- 
ural inodes  so  that  more  than  one  natural  frequency  can  hiave  the 
same  natural  mode  or  even  set  of  modes.  In  many  cases  this  de- 
generacy problem  may  be  unimportant,  especially  for  irregularly 
shaped  objects.  i 

Further  on  in  this  note  some  calculational  techniques  for 
finding  the  natural  frequencies  and  modes  from  general  integral 
equation  formulations  of  electromagnetic  interaction  problems 
are  discussed.  In  performing  such  calculations  one  can  observe 
for  each  natural  frequency  being  considered  whether  other  than 
first  order  poles  are  present  and  if  more  than  one  natural  mode 
is  needed  for  that  natural  frequency.  Thus  the  form  of  the  sin- 
gularity representation  can  be  checked  in  a problem  being  worked 
ovt  and  the  results  compared  to  more  conventional  numerical  so- 
lutions at  other  frequencies  in  the  s plane  and/or  in  the  time 
domain. 

Having  the  natural  modes  va(r')  which  depend  only  on  the 
object  coordinates  and  the  natural  frequencies  sa  which  are 
fixeg  complex  numbers  we  next  need  the  coefficients  which  multi- 
ply v0(ra)(s  - sarn«  to  give  the  response  to  our  incident  plane 
waves.  For  the  surface  current  density  on  the  body  we  write 

i*<r\s)  * S (r ' , s)  + J (r',s)  (2.18) 

s s2  s3 

where  p **  2,  3 as  a subscript  designates  the  part  associated 
with  each  polarization  of  the  incident  wave.  Each  part  can  now 
be  written  as 


(r’,s) 

*P 


<r\s) 


(2.19) 


where  Up  is  the  response  of  the  surface  current  density  to  the 
Up  plane  wave  (Laplac^  transformed  delta^function  wave)  taken  as 
tne  incident  electric  field.  Mote  that  Up  is  dimensionless. 

For  cases  with  volume  current  densities  other  normalizations 
would  be  appropriate.  Our  surface  current  density  response 
functions  (two  of  them,  one  for  each  incident  glectric  polariza- 
tion) are  functions  of  the  object  coordinates  r*  and  the  complex 
frequency  s (or  of  time  t when  inverse  transformed)  and  of 
course  depend  on  the  direction  of  incidence  ei. 


Now  we  can  write  the  surface  current  density  response  func- 
tions for  finite  sized  perfectly  conducting  bodies  as 
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up  8 (r* ,s) 


(2.20) 


where  % has  dimensions  of  (time)~na  and  where  the  sum  is  taken 
over  all  the  indices  in  the  index  set  a except  for  p which  is 
one  index  for  n.  Were  a written  as  say  two  or  three  indices 
then  a double  or  triple  sum  would  be  used  in  this  expansion. 
Note  that  the  coupling  to  the  incident  plane  wave  depends  on 
frequency  because  there  is  at  least  a tima  delay  or  advance  in 
when  the  mode  is  "turned  on." 


The  coupling  coefficients  na(al»s)  a*e  entire  functions  of 
s (no  singularities  in  the  finite  s plane)  with  values  at  sa 
which  give  the  proper  pole  coefficients.  Note  that  for  na  > 1 
the  derivatives  of  ha  with  respect  to  s at  sa  take  a pole  of 
order  r»a  and  give  coefficients  to  terms  of  order  na  - 1#  na  - 2, 
etc.  until  a first  order  pole  is  reached.  Then  there  is  some 
flexibility  i$  our  definition  of  na  long  as  at  each  natural 
frequency  na(el*sa)  gives  the  proper  coefficient  to  the  highest 
order  pole  there.  The  lower  order  poles  at  sa  then  can  be  par- 
tially (if  not  entirely)  included  in  the  terms  for  t|je  poles  of 
higher  order  depending  on  the  choice  of  the  form  of  Ha  for  such 
higher  order  poles.  This  points  out  what  might  be  termed  the 
non  uniqueness  of  the  form  of  the  singularity  expansion.  Cer- 
tain features  of  the  singularity  expansion  are  fixed,  but  others 
have  some  flexibility.  One  might  then  ask  what  is  the  best  def- 
inition of  na  consistent  with  the  pole  requirements?  This  might 
involve  such  criteria  as  simplicity  of  the  resulting  functions 
in  frequency  and/or  time  domain,  asymptotic  behavior  for  JsJ  -*•  ® 
so  as  to  avoid  poles  at  infinity  as  separate  terms  which  compli- 
cate the  form  of  the  time  domain  expansion,  etc. 

Also  in  equation  2.20  we  include  ftp ' v s ■ (e^,r',s)  as  an  en- 
tire function  of  s containing  none  of  tne  poles  of  the  response 
in  the  finite  s plane.  This  entire  function  is  connected  with 
the  choice  of  the  fia  ana  has  similar  flexibility  in  its  choice. 
Only  the  resulting  sum  need  be  the  unique  solution  for  the  cur- 
rent density  etc.  We  have  some  flexibility  in  how  we  arrange 
the  terms.  In  what  follows  in  this  section  we  consider  na  and  a 
special  form  for  the  coupling  coefficient  as  cae“8t°  in  develop- 
ing some  of  the  consequences  for  the  singularity  expansion.  The 
additional  entire  function  wp  is  only  included  in  some  of  the 
expressions;  it  is  usually  dropped;  it  is  not  needed  for  the 
perfectly  conducting  sphere  discussed  in  appendix  B.  This  func- 
tion is  further  considered  in  section  3.  For  perfectly  conduct- 
ing finite  sized  bodies  the  perfectly  conducting  sphere  results 
suggest  that  a delay  can  be  factored  out  so  that  a coupling  co- 
efficient ca  can  be  written  as 


\ 


17 
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c,(e,)<» 
a l 


-sto(el> 


How  general  this  result  applies  is  not  presently  clear.  My  con- 
jecture is  that  it  applies  to  perfectly  conducting  finite  ob- 
jects if  not  more  general  objects  such  as  lossy  objects , perhaps 
with  some  nuances  based  on  the  order  of  the  pole  being  considered. 
Here  tc  is  the  time  that  a delta  function  plane  wave  first 
touches  the  object  and  is  given  by 


w 


(2.22) 


which  is  illustrated  in  figure  2.2.  For  the  perfectly  conduct- 
ing sphere  this  result  is  immediately  apparent  as  in  appendix  £. 

Let  us  call  na  and  ca  coupling  coefficients  and  t0  the  turn 
on  time.  In  any  event  the  object  response  for  each  and  every 
mode  is  zero  for  t < 0.  Note  that  if  (as  is  often  the  case)  the 
coordinate  origin  is  inside  the  object  of  interest  then  to  is  a 
negative  time  (an  advance)  . The  allowable  forms  for  fja,  the  re- 
sulting individual  terms  in  the  series  (all  forms  giving  the 
same  sum) , is  a very  important  question  in  the  singularity  ex- 
pansion method.  Any  entire  function  of  s times  (s  - sa)na  for 
example  can  be  added  to  na  without  introducing  any  new  poles  in 
the  finite  s plane.  This  question  is  considered  further  in  the 
next  section. 

Now  that  the  surface  current  density  response  functions  for 
finite  size  objects  are  expressed  in  terms  of  natural  modes , 
other  quantities  can  be  similarly  expanded  through  their  rela- 
tionship to  the  surface  current  density.  This  includes  scattered 
fields.  However,  for  electromagnetic  interaction  questions  we 
concentrate  in  this  note  on  the  surface  current  density  and  sur- 
face charge  density.  From  the  continuity  equation 

V*3+|£«°,  V-j+sp-0  (2.23) 


we  can  find  the  charge  density  from  the  current  density.  When 
we  are  dealing  with  surface  current  density  $s  and  surface 
charge  density  p8  the  divergence  has  to  be  interpreted  as  a sur- 
face divergence  with  the  spatial  derivatives  being  taken  with 
respect  to  two  coordinates  required  to  describe  a position  on 
the  surface.  Thus  we  can  still  write  without  ambiguity 
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FIGURE  2.2  TURN-ON  TIME  FOR  THE  NATURAL  MODES 
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+ sps(r’ ,s) 


0 


(2.24) 


where  the  prime  on  the  7 indicates  derivatives  with  respect  to 
the  object  coordinates 


As  before  split  the  surface  charge  density  as 


os(r’,s)  * ps  (r’,s)  + ps  (r  * , s) 
2 3 


(2.25) 


to  correspond  to  the  two  polarizations.  Write  each  part  as 


Pa  (*',■> 
SP 


e E f (s)U 
o op  p 


(r* ,s) 


(2.26) 


where  the  normalization  using  e0E0  is  again  chosen  to  make  Up 
dimensionless.  Prom  the  surface  current  density  natural  modes 
construct  a set  of  surface  charge  density  natural  modes  as 


( P _ ) . . ( J_)  V 

vo  (r'}  " "aaV’  * va  8 (r’)  (2.27) 


where  aa  is  a scale  factor  with  dimension  of  length  which  we  can 
choose  for  convenience,  such  as  to  allow  some  desired  normalize-* 
tion  condition  on  the  va.  If  desired  the  aa  could  bo  all  the 
same  and  perhaps  chosen  as  some  characteristic  dimension  of  the 
object.  Note  that  some  va  for  the  surface  charge  density  might 
be  identically  zero  for  some  index  sets  a.  This  is  possible  be- 
cause current  density  can  be  split  into  two  parts,  one  with  zero 
divergence  but  nonzero  curl  and  one  with  zero  curl  but  nonzero 
divergence  if  the  current  density  is  confined  to  a volume  of 
finite  dimensions. 

All  the  above  definitions  for  surface  current  and  charge 
density  modes  can  be  directly  extended  to  volume  modes  or  com- 
bined modes  for  volume  and  surface  densities. 

Taking  out  common  factors  the  continuity  equation  allows  us 
to  write  the  surface  charge  density  response  functions  (using 
equation  2.20)  as 
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~<PS>  -* 

Up  (t\a\ 


Iv1-  /■*■  (ps>  ,+tv/s*a\  ”Ra)  /-►  +, 

{2j/1a<®l's,va  (r  >1— ) (+Mp  (el'r's) 


(2.28) 


where  the  splitting  of  na  can  be  done  as  before . All  our  pre- 
vious results  for  the  surface  current  density  then  carry  over  to 
the  surface  charge  density.  The  same  set  of  natural  frequencies 
and  coupling  coefficients  apply  to  both  (and  even  to  scattered 
fields  if  you  like).  Note  that  ftp'P8)  comes  from  Wp(Js)  via 
equations  2.23. 


However,  there  is  what  looks  like  a new  pole  at  & * 0 in- 
troduced into  the  surface  charge  density  response  function.  Us- 
ing a relation  for  separating  poles  at  separate  frequencies  we 
have 


s’1(s  - s r1  - 8~x  (s  - s r1  - s'V1 


i -n  , -n  - -n  +1 

8 - «o)  - 8^(8  - «0,  “ - S”2(S  - Ba)  01 


n -1  -n  i -n  , 

+ •••  + (-1)  a sa  c(s  - s0)’1  + (-so)  % 1 


(2.29) 


By  this  expansion  we  can  write  the  surface  charge  density  re- 
sponse function."  as 


n(Ps)/i.  „ V*  ,_^m-l  c _-m-  ,(ps)  ,i.w_  - v -n0+m-l 

up  vr  ** 


a m»i  a 


+ iEK1  Vel's)va  (r,)  +%  (err',s) 

(2.30) 


where  th*s  second  summation  could  also  be  included  in  a. 

Consider  for  a moment  the  static  response  characteristics 
(b '+  0)  of  a body  of  finite  dimensions.  For  a plane  wave  of 
unit  amplitude  as  s + 0 the  object  response  goes  to  the  static 
limit  in  which  both  current  density  and  charae  density  are  pro- 
portional to  the  field  strength.  For  small  |sj  the  response  is 
negligibly  changed  from  the  s » 0 case.  Thus  there  is  no 
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singularity  at  8 * 0 in  the  s plane  for  either  current  density 
or  charge  drnsity.  As  discussed  before  for  the  current  density 
there  are  no  poles  with  nonzero  residue  on  the  iw  axis  for 
| U | >0  because  of  the  restriction  of  conservation  of  energy  to- 
gether with  reciprocity  when  considering  the  time  domain  re- 
sponse. The  static  response  characteristics  rule  out  poles  in 
the  response  at  s - 0 so  that  all  poles  sQ  lie  in  the  left  half 
pla.se  He[sj  < 0.  Thus  without  loss  of  generality  in  our  expan- 
sions we  can  require 

Re[saJ  < 0 (2.31) 

Cavity  modes  or  any  other  modes  which  do  not  couple  to  the  inci- 
dent field  are  excluded  from  our  consideration. 

Note  that  there  are  static  current  density  and/or  charge 
density  distributions  which  can  exist  on  finite  sized  objects  in 
the  absence  of  any  incident  wave.  However,  such  static  distri- 
butions do  not  couple  to  the  incident  wave  and  can  be  included 
as  an  additive  term  in  the  surface  current  density  and  surface 
charge  density  response,  but  with  no  dependence  on  the  incident 
waveform.  As  the  response  to  the  incident  wave  does  not  depend 
on  these  modes  we  do  not  include  them  i-  our  expansion.  We  can 
call  this  case  the  natural  frequency  a«-  s «=  0 and  this  is 
briefly  considered  in  appendix  A.  This  case  is  like  the  inter- 
nal cavity  modes  which  have  natural  frequencies  on  the  axis 
s » iw;  these  also  do  not  couple  to  the  incident  wave  and  have 
no  dependence  on  the  incident  waveform;  they  can  be  added  into 
the  -results  at  the  end  if  desired. 

Referring  back  to  equation  2.30  for  the  surface  charge  den- 
sity response  functions  note  that  a pole  at  b ■*  0 is  not  allowed 
Thus  we  have  letting  s 0 in  the  limit 

E"not  c ^ <PS) 

(sa)  “ %<«i»0>va  (r,)  + Wp  <V*\0)  - 0 (2.32) 

a ci  “ 

For  cases  that  n(  can  be  made  to  factor  (such  as  for  the  sphere) 
as  in  equation  2.21  then  we  can  write  (also  dropping  ftp) 

^ -n  e + (ps)  + 

jr  c«(ei)v0  !r’>  - 0 <2-33) 

a a 

This  constrains  a relation  among  the  natural  modes  and  their 
coupling  coefficients. 


For  the  case  that  net  **  1 for  all  a (which  Is  th«  oasa  for 
the  perfectly  conducting  sphere,  end  perhaps  for  all  finite 
size,  perfectly  conducting  objects)  and  that  the  coupling  fac- 
tors as  in  equation  2.21,  the  response  functions  for  surface 
current  density  and  surface  charge  density  can  be  written  as 
(dropping  the  W functions) . 


T . -St  .1^  .(2-)  . <| 

5p  S (r',s)  = e °Eca(aj)vas  (r')(s  - sa) 

D (f,s)  - o (r-) 

r a a ct 


(2.34) 


This  is  a rather  simple  looking  result  with  each  term  quite  fac- 
tored. This  also  points  out  the  importance  of  understanding 
under  what  circumstances  (as  for  the  sphere)  the  delay  factoring 
of  the  coupling  coefficient  can  be  used. 

Now  the  waveform  functions  ?p(s)  can  be  reintroduced  and 
multiplied  on  both  sides  of  equations  2.34  or  equations  2.20  and 
2.30.  Recombining  the  polarizations  as  in  equations  2.18,  2.19, 
2.25,  and  2.26  gives  complete  representations  of  the  solutions 
for  surface  current  density  and  surface  charge  density.  This 
would  give  individual  terms  with  frequency  dependence  in  the 
form  f!p(s)np0(ei,s)  (s  - s<x)  -na+m-l  which,  for  the  case  as  before 
that  5 factors,  has  the  frequency  dependence  in  the  form 
e“8^°|p(s)  (s  - sa)“na+R'“1  which  for  na  * 1 reduces  to 
e”stOfp(s) (s  - sa)-l. 

However,  why  stop  here?  The  behavior  of  fp(s)  may  allow  us 
to  conveniently  express  i’  in  terms  of  waveform  singularities, 
just  as  we  have  been  considering  the  object  delta  function  re- 
sponse in  terps  of  object  singularities.  The  general  idea  then 
is  to  expand  fp(s)  in  terms  of  its  singularities  and  separate 
the  waveform  and  object  singularities  into  separate  terms.  We 
might  call  the  resulting  separate  terms  as  the  waveform  part  and 
object  part  of  the  response  for  convenience.  Consider  an  ex- 
ample by  letting 


(s) 


8 t ' 

fp(t)  - e w u (t) 


(2.35) 


so  that  the  waveform  is  a decaying  exponential  with  a simple 
pole  at  sw  (a  waveform  pole) . Note  that  the  commonly  used 
double  exponential  waveform  for  EMP  environments  is  nothing  more 
than  the  sum  of  two  terms  such  as  this. 
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Taking  the  surface  current  density  first  we  have  (for  sw  jt 
sa  for  any  o) 

V (r',s)  i ? (s)$  8 (r',s)  - 2na(el's)vo  <r'Xs-sw) 


*(3.)  4.  i(5s)  4. 

= V 8 (r' ,s)  + $ 8 (rSs) 


(2.36) 


w 


where  the  object  part  is 


~($  ) _ (3  ) -n  +m-1 

■ (r',t)  - 2 2(-l)w  1(«^sM)m?ia(e1.s)v  8 (r')(s-sa)  ° 

° a m*l 

(2.37) 

where  the  second  sum  could  also  be  included  in  a and  where  the 
waveform  part  is 


*(Ss)  * 

"p  <*' •»>  - T 


w 


C 3 j 

hrZ<%  - v'X'5!'”5!  8>(?’»  (2-3e> 

” a 


Note  that  here  and  for  most  of  what  follows  the  W functions  are 
not  included  but  can  be  included  in  a convenient  way  in  the  ob- 
ject-waveform split.  For  the  special  case  of  nB  " 1 this  re- 
duces to 


*<3.»  .* 


<i.) 


4.  o .1 4.  4.  'a'  -1 

^ (r’,8)  - 2w<s a’V  \{el's)K  <?,><8-sa> 


(2.39) 


w a 


For  more  general  waveforms  (but  still  na  - 1)  this  suggests 
a definition  of  the  splitting  as 

3pos)(r',s)  - ^fp(sa)na(elf8)^  s>(r')(t-sa)'1 

(2.40) 
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V B lr’.»>  “ 2.  "a<Vs)va  (r')-E-a--»  ■ ■ 


As  long  as  the  waveform  singularities  are  separate  from  all  sa 
then  no  object  poles  appear  in  the  waveform  part  of  the  response. 
If  the  Sa  do  lie  on  waveform  singularities  then  special  treat- 
ment is  needed  but  the  general  idea  of  equations  2.40  would 
still  seem  appropriate.  For  ntt  more  than  just  1 the  object  part 
can  be  defined  through  a Taylor  expansion  of  ?p  around  each  sa 
giving 


*(3_)  - & f^m"1){s  ) _ ^ _(3  ) + 

^ tr-,s)  s E E-^nrr—  Ws)',a  (r,>(9-90> 

ro  a m-i 

*(Ia)  * »<3.)  i(3.)  . 

? 9 <r\s)  5 $ 9 <r’,s)  - ? 9 (r',s) 


-n0+m-l 


(2.41, 


■w 


sw 


As  a special  case  consider  the  unit  step  waveform  by  letting 
0.  Then  the  waveform  part  of  the  response  we  write  as 


${  s\r',s)  = lIX^s)^  8>(r')(  sa)~1 


(2.42) 


With  na  factored  as  before  and  assuming  n0  =>  1 the  results  re- 
duce to 


-st-  ,<S.) 


X (3  ) ~~Q  X > 

V„  8 (r\s)  = - 3 8 (r’,0) 


w 


s p 

-St. 


% (r’'s)-e  Esa  ca(ei)va  (r')(s-stt)  1 


(2.43) 


where  the  static  surface  current  density  response  is 

S’^NtSO)  - E<-90>'lca'9l!,(f,)<?'> 


(2.44) 


The  step  function  respor.se  of  the  object  is  readily  constructed 
from  the  delta  function  response.  A factor  of  (-sa)"*  is 
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included  with  each  object  pole,  and  then  a new  term  is  added 
which  is  nothing  more  than  the  static  response  time  a step  func- 
tion with  turn-on  time  of  tQ.  Thus  beside  tabulating  the  natu- 
ral frequencies,  modes,  and  coupling  coefficients  of  finite  size 
objects  we  can  tabulate  the  static  response  so  as  to  readily 
construct  the  step  response  of  the  object  in  both  Laplace  and 
time  domains. 


Now  the  neat  thing  about  the  static  response  is  that  we 
need  not  consider  it  in  terms  of  direction  of  incidence  and  pol- 
arization. For  the  static  surface  current  density  r^ponse  we 
can  solve  for  the  magnetic  field  using  the  Laplace  equation  with 
a uniform  magnetic  field  incident  on  the  object.  The  exponen- 
tial factor  in  the  incident  plane  wave  (equations  2.11)  becomes 
irrelevant,  going  to  1.  Thus  we  need  only  consider  three  sepa- 
rate cases  of  the  incident  magnetic  field,  corresponding  to 
three  orthogonal  direction  such  as  the  cartesian  axes  (x,  y,  z) . 
£or  a unit  incident  static  H field  we  have  Js  response  functions 
USx  etc.  for  each  of  the  three  axes  giving  a dyadic  surface  cur- 
rent density  response  function  as 

j<3.>  » -<53>  H.  „ _<3S)  » _<JS>  „ , 

UB  8 (r'l  5 Ua  3 (r’)ex  + Us  3 (r')ey  + Us  8 (r')ez  (2.45) 
x y z 


Since  the  a field  is  in  the  direction  ei  x en  for  p = 2,  3 we 
can  write 


~(3  ) -*-(3  ) 

“fig  8 (r')  [ex*  (e^p)  ]+6g  8 (r ' )(ey«  («1xep)  )+fiB  8 (r ' )[e2*  (exxep)  ] 


(2.46) 


showing  a set  of  direction  cosines  for  three  static  solutions. 
The  static  response  dyadic  for  the  surface  current  density  is 
then  a useful  tool  for  directly  extending  the  delta  function  re- 
sponse to  the  step  function  response.  Cne  could  go  to  other 
types  of  "static"  terms  for  s~2  (ramp)  and  higher  order  such 
waveforms.  The  step  function  waveform  is  a useful  tool  because 
of  its  fast  rise  time  to  a finite  amplitude  making  the  static 
response  correspondingly  important. 

The  surface  charge  density  has  similar  properties  when  the 
waveform  is  reintroduced  as 


( 


( 


( 
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„(PB>  - „<pa>  ^ Jpg)  ^ ,(pa)  .4 

V 8 (r‘,s)  =fD(s)U  8 <r',s)  -V  8 (r',s)  +V  8 (r» ,s)  (2.47) 

P P P PQ  Pw 


As  be£ore  the  idea  is  to  associate  the  object  poles  with  the  ob- 
ject part  of  the  response  Vpo  and  the  rest  with  the  waveform 
part  Vpy.  If  na  takes  on  more  values  than  one  the  double  sum  as 
in  equation  2.30  extends  into  a triple  sum  including  coeffici- 
ents £p(s)  and  the  first  na  - 1 derivatives  all  evaluated  at  sa. 

For  simplicity  just  consider  the  case  of  n<*  ■ 1 with  fja 
factorable  as  above  so  that  we  have  the  simple  surface  charge 
density  response  function  as  in  equations  2.34.  The  object  and 
waveform  parts  of  the  response  may  then  be  written  as 
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(r'rs)-c  °XX<8a> 


p a s a„  ala 
a a a 


e 4 (Pg)  ^ 

c»(e,)v„  8 (r'Hs-sl  1 
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(2.48) 


x "or  c .■*>  . %Ks#  p'"'  p'  a 
(r',s)=e  2^iV"  ca(el)va  (r»)  ^ • 


a a a 


The  response  to  an  exponential  waveform  as  in  equations  2.35 
takes  the  forms 


V (r-,s)=e  JlvV  rrWa 


-1  c , ,-l 

cJe,)vM  (r*)(s-sa)x 


-(c8>  - e’St° 

vp 
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r~ r-  c (e.) v (r‘) 
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(2.49) 


For  sw  » 0 we  have  the  step  response  as 


B)-e"toE 
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(2.50) 


where  the  static  surface  charge  density  is  given  by 
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(2.51) 


~(PS>  - 

Op  8 (r ' ,0) 


E (-.f 1 


-1  c 


8 a 
aa 


+ (Ps)  * 

».«*!>*.  <*'> 


The  step  function  surface  charge  density  response  is  then  also 
constructed  from  the  delta  function  response  In  the  same  way  as 
the  surface  charge  density  with  the  same  form  of  the^results,  at 
least  for  na  « 1 and  a common  delay  factorable  from  ria» 

For  the  static  surface  charge  density  response  we  need  con- 
sider only  the  response  to  three  orthogonal  incident  static  elec- 
tric fields,  directly  analogous  to  the  relationship  of  the 
static  incident  magnetic  field  and  surface  current  density.  For 
a unit  incident  static  E field  we  have  ps  response  functions 
(scalars)  as  USv  etc.  for  each  of  the  three  axes  giving  a vector 
surface  charge  density  response  function  as 


jj'rS'  ,1,,  - ..'"8'?  . . „*rS 


<Pj  ■+  <PJ 

(r ' ) s o 
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<PS>  + 

e + U “e  + U e 
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(2.52) 


Since  the  E field  is  in  the  direction  ep  for  p * 2,  3 we  can 
write 


-(pB)  -►  _JPS)  + + 

Up  8 (r ' ,0)  - 0a  8 (r * ) «ep 


(P.) 


<PJ 


(pJ 


= °B  Iex*V+Us  Jey*V+Us  [ez*V  (2*53) 

x y z 


showing  a set  of  direction  cosines  to  weight  the  three  static 
solutions.  Thus  for  step  response  purposes  it  is  useful  to  tab- 
ulate the  vector  surface  charge  density  response  and  the  dyadic 
surface  current  density  response  which  can  be  multiplied  in  a 
dot  or  inner  product  sense  with  an  appropriate  scale  factor 
times  the  static  field  of  interest  to  obtain  the  static  response. 
This  applied  in  both  frequency  domain  for  small  js|  and  in  time 
domain  for  an  important  term  in  the  step  function  response. 

As  discussed  near  the  beginning  of  this  section  the  fact 
that  we  are  dealing  with  real  valued  time  functions  makes  the 
Laplace  transformed  functions  have  certain  syimnetry  properties 
with  respect  to  the  ft  axis  as  expressed  by  equation  2.4.  Bas- 
ically ?(Q  + iu>)  has  its  real  part  symmetric  with  respect  to  a 
sign  reversal  of  u while  its  imaginary  part  is  antisymmetric 
with  respect  to  a sign  reversal  of  w.  Also  as  discussed  before 
all  object  poles  with  nonzero  coupling  to  the  incident  wave  lie 
in  the  left  half  plane  ft  < 0.  These  two  results  tell  something 


28 


about  the  object  pole  pattern  in  the  complex  s plane.  For  con- 
venience then  let  us  split  the  a index  set  into  three  parts. 

For  w > 0 (the  upper  half  plane)  use  ct+,  for  w ■ 0 (the  ft  axis) 
use  a0,  and  for  to  < 0 (the  lower  half  plane)  use  a-.  Since 
poles  come  in  conjugate  pairs  for  those  not  on  the  ft  axis  then 
we  can  relate  a-  to  a+  as 


(2.54) 


which  specifies  which  a.  goes  with  which  a+  except  for  the  case 
of  multiple  poles  at  sa  in  which  case  we  make  the  identification 
of  the  a-  set  to  the  a+  set  with  the  value  of  n«  in  each  case 
the  same  so  that  we  require 

V “ na+  (2*55) 


Equations  2.54  and  2.55  define  a one  to  one  correspondence  be- 
tween a-  and  a+  index  sets,  unless  we  have  mode  degeneracy  in 
which  case  we  also  require  a-  and  ot+  correspond  to  conjugate 
modes  with  conjugate  coupling  coefficients  as  well. 

From  the  conjugate  symmetry  requirement  the  natural  modes 
and  coupling  coefficients  can  be  made  to  have  the  same  conjugate 
relations.  Thus  we  set 
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With  these  relations  we  can  now  write  the  surface  current  den- 
sity response  functions  as 
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'Vhere  the  case  where  fia  factors  out  a simple  delay  is  also  in- 
cluded- Note  now  that  each  term  in  the  sum  has  the  conjugate 
symmetry  of  equation  2.4,  is  real  on  the  ft  axis,  and  corresponds 
to  real  valued  time  function.  Similar  properties  apply  to  the 
surface  charge  density  response  functions.  Taking  the  simpler 
form  with  fia  factored  as  in  equations  2.34  with  only  n0  * 1 we 
ban  write 
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The  same  idea  is  readily  applied  to  the  more  complex  forms. 

If  rig  factors  as  in  equation  2.21  then  the  s dependence  of 
each  term  is  a delay  with  a term  of  the  form  (s  - sa)“n  where 
n’  might  be  na  or  other  more  complicated  exponents.  Converting 
this  term  to  the  time  domain  through  the  well  known  Laplace 
transform  pair  gives 
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For  the  extremely  interesting  case  of  n'  *1  this  is 
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This  is  a damped  sinusoidal  waveform  like  that  which  suggested 
looking  at  the  natural  frequencies  in  the  first  place.  The 
damping  constant  is  just  fta  = Re[sa)  (plus  or  minus  as  one 
wishes)  and  the  radian  oscillation  frequency  is  just  wa  « Im[sa] 
where  a world  be  taken  as  a+  or  aQ  for  this  frequency  to  make 
«o  1 0.  Since  fla  < 0 for  all  poles  of  interest  then  each  term 
in  the  response  function  goes  to  zero  for  large  time  as  we  would 
expect.  Writing  out  the  real  and  imaginary  parts  gives 
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If  the  oa  terms  do  not  factor  as  above  then  we  need  their  in- 
verse transforms  to  convolute  with  terms  similar  to  these. 


Consider  then  the  surface  current  density  delta  function 
response  (for  n factored) . From  equation  2.57  we  have 
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For  the  surface  charge  density  for  n„  ■ 1 we  have  the  delta 
function  response  as 
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From  the  way  that  the  terms  are  split  up  we  have  separated  the 
damped  exponential  parts  of  the  response  from  the  damped  oscil- 
latory parts.  Depending  on  the  real  and  imaginary  parts  of  the 
natural  modes,  natural  frequencies,  and  coupling  coefficients 
the  initial  "phase  angle"  of_^the  sinusoidal  oscillation  may  vary 
yith  direction  of  incidence  ei,  polarization  p,  and/or  position 
r*  on  the  object.  If  one  wished  he  might  consider  the  real  and 
imaginary  parts  of  a natural  mode  as  separate  modes  and  consider 
the  coupling  to  each  part  with  different  distributions  perhaps 
over  the  object. 

With  the  incident  waveforms  reintroduced  we  can^also  di- 
rectly write  down  response  waveforms  for  cases  that  n0  factors 
as  before.  Consider  the  case  that  no  waveform  singularities  are 
at  any  of  the  sa.  In  particular  consider  the  step  function  re- 
sponse for  all  na  « 1.  For  the  surface  current  density  the  ob- 
ject and  waveform  parts  from  equations  2.43  and  2.46  are 
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For  the  surface  charge  density  the  object  and  waveform  parts  of 
the  step  response  functions  from  equations  2.50  and  2.53  are 
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In  the  time  domain  only  one  term  is  added  to  the  delta  function 
response,  namely  a step  function.  The  coefficients  of  the 
damped  sinusoidal  (and  simple  exponential)  terms  are  rather 
simply  altered. 

The  response  functions  for  other  kinds  of  incident  wave- 
forms can  be  readily  found  for  the  object  part  by  introducing 
coefficients  fp(sa)  in  the  time  domain  waveforms  in  the  same 
forms  as  introduced  in  the  Laplace  versions  such  as  in  equations 
2.41  and  2.48.  The  waveform  part  can  be  more  difficult  if  only 
( because  of  the  many  possibilities  one  might  choose  for  incident 

waveforms.  Different  incident  waveforms  give  different  types  of 
time  domain  waveforms  when  combined  with  the  terms  arising  from 
the  object  poles.  Note  that  the  response  to  an  incident  wave- 
form cannot  always  be  simply  split  into  object  and  waveform 
parts.  As  a trivial  example  suppose  £p(s)  itself  has  a pole  at 
some  sa.  In  such  a case  the  contribution  from  sa  to  the  re- 
sponse has  a higher  order  pole  than  the  delta  function  response. 
However  the  response  associated  with  sa  is  then  easily  treated 
separately  and  the  same  type  of  time  domain  function  as  in  equa- 
tion 2.59  results. 

In  this  section  we  have  tried  to  give  some  insight  into  the 
power  of  the  singularity  expansion  method  for  representing  solu- 
tions in  both  frequency  and  time  domains,  at  least  for  finite 
sized  objects.  While  many  variations  of  this  problem  have  been 
considered,  complicating  the  notation  somewhat,  the  solution  of 
specific  problems  may  be  expressed  somewhat  simpler  with  the  no- 
tation adapted  to  the  results  of  the  problem  at  hand.  Specific 
convenient  choices  of  the  o index  set  are  needed  for  each  prob- 
lem and  the  possibilities  for  na  can  be  limited.  The  natural 
modes  and/or  coefficients  may  be  expressible  as  purely  real 
quantities  in  some  cases,  even  for  natural  frequencies  off  the  0 
axis.  A basic  question  concerns  the  coupling  coefficients  na 
and  the  other  entire  functions  wp.  These  need  to  be  optimally 


33 


chosen  for  the  different  problems  at  hand  depending  on  early 
times,  late  times,  etc. 

There  are  still  types  of  problems  which  have  singularity 
expansions  with  terms  such  as  branches  which  have  not  been  con- 
sidered here.  Infinite  sized  objects  have  object  responses  with 
such  terms.  However  the  reader  should  have  a general  idea  by 
now  of  what  the  singularity  expansion  method  is  all  about,  tfhat 
has  been  done  for  expanding  in  terms  of  natural  frequencies  for 
object  response  plus  another  term  for  waveform  response  can  be 
carried  over  to  natural  frequencies  and  branches  for  the  object 
response  with  perhaps  increased  complexity  associated  with  the 
branches . 
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III. 


Properties  of  the  Singularity  Expansion  of  the  Object 
Response  for  Finite  Size  Objects  Viewed  from  the  Finite 
Matrix  Formulation  of  Integral  Equations 


There  are  two  somewhat  complementary  ways  to  view  some  of 
the  questions  regarding  the  form  of  the  singularity  expansion 
for  objects  of  finite  size.  One  approach  involves  consideration 
of  the  properties  of  an  integral  equation  formulation  for  the 
continuous  object  geometry.  Various  types  of  integral  equations 
such  as  those  classed  as  electric  field  formulation,  magnetic 
field  formulation,  extended  boundary  condition  formulation,  etc. 
cai  be  investigated  to  obtain  mathematical  theorems  concerning 
the  properties  of  the  singularity  expansion  for  various  kinds  of 
object  classes.  This  might  be  termed  the  continuous  integral 
operator  approach.  Drs.  Marin  and  Latham  are  presently  using 
the  magnetic  field  integral  equation  to  consider  the  question  of 
there  being  only  poles  in  the  finite  s plane  in  the  object  re- 
sponse for  finite  size  perfectly  conducting  object  from  this 
viewpoint.  Clearly  this  kind  of  approach  is  needed  for  consid- 
ering many  such  questions  so  as  to  establish  general  theorems 
applying  for  all  frequencies  and  time  and  stated  exactly  in  an 
analytical  form.  This  could  be  considered  a viewpoint  which  is 
based  on  the  continuous  nature  of  the  object  geometries.  Eigen- 
function expansions  can  also  be  used  for  such  considerations  but 
the  tabulated  cases  of  such  expansions  are  limited.  Eigenfunc- 
tion expansions  are  possible  for  general  kinds  of  object  geom- 
etries but  they  must  be  numerically  calculated. ? However  the 
analytic  properties  of  such  eigenfunction  expansions  can  still 
be  used  to  investigate  the  singularity  expansion.  Viewed  an- 
other way  such  eigenfunction  expansions  are  representations  of 
the  integral  operators  defined  over  the  object  characteristics. 

The  second  and  complementary  approach  might  be  termed  the 
discrete  approach.  This  refers  to  zoning  the  object  into  many 
discrete  zones  and  treating  each  zone  as  a position  with  a par- 
ticular current  density  etc,  associated  with  it.  The  integral 
equation  then  takes  the  form  of  a vector-matrix  equation  which 
has  the  general  form 

(g_  ~ (s) ) • (J(s))  = (I_(s) ) n,m  - 1,  2,  ••*,  N (3.1) 

lu  n 

where^N  is  some  typically  large  integer.  Here  the  vectors  (Jm) 
and  (In)  each  have  N components  and  are  not  the  same  as  th<g 
three  component  space  vectors.  The  index  n refers  to  the  r co- 
ordinates for  the  "incident"  quantities  shown  here  in  general  as 
(in)  which  might  come  from  an  electric  or  magnetic  field  that  is 
incident  and  therefore  specified.  Note  that  in  forming  the  lin- 
ear equations  summarized  as  equation  3.1  each  spatial  zone  can 
have  one,  two,  or  three  nonzero  current  density  components.and 
this  influences  how  the  indices^n  and  m are  set  up.  The  (Jm) 
refers  to  the  unknown  current  (r*  coordinates)  with  one,  two,  or 
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three  components  in  each  zone.  The  matrix  (gn,m)  has  elements 
which  couple  the  n zone-component  with  the  m zone-component;  it 
is  basically  a spatially  discrete  form  of  the  integral  operator. 
While  equation  3.1  is  shown  in  a form  suggesting  that  it  is  in- 
tended for  solving  for  the  current  density  this  is  not  necessar- 
ily the  case;  it  could  be  solving  for  another  quantity  from 
which  the  current  density  would  be  obtained  by  a subsequent  cal- 
culation. For  our  present  discussion,  however,  consider  it  the 
current  density. 

Considering  the  current  density  components  in  each  zone  as 
our  unknowns  is  not  the  only  way  to  obtain  a set  of  linear  equa- 
tions from  a given  integral  equation.  The  current  density  can 
be  expanded  in  a more  general  set  of  expansion  functions  and  a 
matrix-vector  equation  formed  to  obtain  the  coefficients  in  this 
expansion.  Of  course,  only  a finite  set  is  used  to  obtain  a fi- 
nite N x N matrix  and  the  sets  of  functions  involved  should  be 
in  some  sense  complete  in  the  \imit  of  large  N.  This  general 
approach  is  often  termed  the  method  of  moments. 8 Both  the  zon- 
ing approach  and  more  general  function  expansions  are  valuable 
from  the  viewpoint  of  numerical  calculations.  In  terms  of  the 
singularity  expansion  method  there  is  clearly  much  work  to  be 
done  in  refining  the  numerical  techniques  to  find  natural  fre- 
quencies, modes,  etc.  most  efficiently  and  most  accurately. 

Much  that  has  been  done  for  other  numerical  problems  can  likely 
be  applied  here. 

Besides  the  practical  aspect  of  numerical  computations  the 
matrix-vector  formulation  can  be  used  as  a theoretical  technique 
for  establishing  some  of  the  general  characteristics  of  singular- 
ity expansions.  If  in  the  limit  of  large  N the  solution  for  the 
current  density 


(Jm(s))  » (g  I"1  • (I)  (3.2) 

m n,ro  n 

converges  to  the  exact  solution  of  the  integral  equation  then  by 
understanding  general  properties  of  (Jm)  we  can  find  general 
properties  for  the  continuous  case,  i.e.  J(?',  s) . Only  where 
the  matrix  elements  are  not  uniquely  defined  or  the  inverse  ma- 
trix 


F 5 5 'Vn,!'1  <3-3> 

does  not  exist  is  (Jm)  not  uniquely  defined  or  non  existent. 

Let  us  then  look  at  the  singularity  expansion  characteris- 
tics of  the  approximate  numerical  solution  (Jra) . In  the  limit 
of  large  N for  convergent  matrix-vector  formulations  the 


singularity  expansion  of  (Jm>  will  be  that  of  s) . Con- 

sider the  case  that  some  zoning  technique  is  used  to  convert  the 
integral  equation  to  a matrix-vector  equation  by  dividing  the 
body  surface  or  volume  as  appropriate  into  discrete  zones  of  fi- 
nite linear  dimensions.  The  number  of  zones  is  then  less  than 
or  equal  to  N.  Clearly  then  we  are  only  considering  finite 
sized  objects  here  because  an  infinite  surface  area  or  infinite 
volume  cannot  be  divided  up  into  a finite  number  of  areas  and/or 
volumes  with  finite  linear  dimensions  unless  some  very  strange  (,V 

extensions  to  infinity  with  finite  volume  and/or  surface  area  j 

are  included.  Let  u$  consider  only  finite  sized  objects  here  so  ' ■ 

that  all  d^stanc§s  Ir^yi-  r»l  between  z$nes  are  finite.  The  co-  *1 

ordinates  r and  r'  go  over  the  rg>nand  rp^nd  are  centered  on  ^ 

each  zone  in  some  sense.  The  incident  quantities  and  resulting 
current  density  are  evaluated  in  an  appropriate  average  sense  at 
each  zone  center. 

1 

Let  f(r,  s)  be  based  on  an  incident  delta  function  plane 
wave.  It  might  come  from  the  electric  and/or  magnetic  field  or 
some  spatial  projection  of  these  on  the  object.  Since  this  is 
an  entir$  function  of  s (analytic  in  the  entire  finite  s plane 
for  all  r on  the  finite  size^object)  then  the  discrete  formula- 
tion (In)  has  each  element  (in(s))  as  an  entire  function  of  s. 

In  this  formulation  let  the  current  density  be  the  sum  of  all 
terms  representing  charge  motion,  including  oE  (conduction  cur- 
rent density)*  s(e  - e0)fi  (displacement  current  density),  and  ! 

V x [ (y  - v>0)  h]  (magnetization  current  density).  Split  electric 
and  magnetic  fields  each  into  the  sum  of  incident  plus  scattered 
parts.  The  incident  parts  are  given  in  the  problem  definition. 

The  scattered  parts  can  then  be  represented  as  integrals  over 
the  current  density  with  kernels  involving  the  free  space 
Green's  function  which  uses  y « s/c  instead  of  some  propagation 
constant  involving  local  medium  parameters  which  may  vary  with 
position. 

Thus  we  have  a pair  of  volume  integral  equations  (which  may 
reduce  to  surface  integral  equations  say  for  perfectly  conduct- 
ing bodies)  which  equate  scattered  electric  and  magnetic  fields  . . 
(which  could  be  thought  of  as  a six  component  vector  or  even  a j ; ' 

four  tensor  in  relativistic  formulation)  to  integrals  involving  co- 
incident plus  scattered  fields.  Move  the  terms  involving  inci-  ftf'' 
dent  fields  to  one  side  of  the  equations  and  use  these  to  form  w V 
the  N component  vector  (in)  . Note  that  o,  e,  and  y are  assumed  4 
to  be  single  valued  analytic  functions. of  s except  possibly  for 
poles  which  can  be  removed  by  multiplying  through  the  equations 
by  an  appropriate  zero  to  make  each  in(s)  an  entire  function. 

The  scattered  parts  are  used  to  form  (gn,m)  and  (Jm>  which  -an  Is 
appropriate  scattered  field  vector  or  the  current  density  vector 
less  that  part  directly  proportional  to  the  incident  fields. 

Note  that  a,  c,  and  y enter  the  coefficients  on  the  scattered 
side  and  are  analytic  single  valued  here  then  as  well  while  the 
exponential  terms  use  y ■ s/c  which  is  an  entire  function. 
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Through  thin  procedure  we  can  construct  an  integral  equation  re- 
lating current  density  and/or  surface  current*  density  as  appro- 
priate to  the  incident  fields  involving  only  analytic  Bingle 
valued  functions  as  long  as  the  media  are  suitably  well  behaved. 
By  various  manipulations  such  a pair  of  integral  equations  can 
be  converted  into  various  more  desirable  forms.  Then  the  in- 
tegral equations  can  be  converted  into  a single  matrix-vector 
equation  which  can  have  various  forms. 

One  advantage  of  a zoning  formulation  for  these  considera- 
tions is  the  somewhat  physical  picture  one  can  associate  with 
the  discrete  zones.  Infinitesimal  differential  line,  surface, 
and  volume  elements  become  discrete  ones  of  small  size.  As  long 
as  complex  radian  wavelengths  have  magnitude  large  compared  to 
zone  size  then  the  interaction  between  adjacent  zones  and  one 
zone  on  itself  are  described  by  statics.  We  can  start  thinking 
of  the  zoned  object  as  a big  circuit  with  simple  elements  con- 
necting adjacent  zones  but  more  complex  ones  connecting  distant 
zones.  In  each  zone  there  are  a few  equivalent  sources  associ- 
ated with  the  discrete  source  elements  in. 

The  kernels  of  the  various  integral  equations  are  based  on 
the  scalar  Green's  function  (for  free  space)  as 


G(s,  r-r'l) 
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(3.4) 


with 


(3.5) 


This  sca^a£  Green's  |u§ction  is  an  analytic  function  of  s for 
finite  |r-r'  | with  |r-r'|  jrf  0.  This  function  forms  the  basis 
for  the  various  kernels  used  in  the  various  types  of  integral 
equation  formulations  of  general  interaction  or  scattering  prob- 
lems. The  dyadic  Green's  function  can  be  written  in  the  form9#10 
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Note  that  the  dyadic  Green's  function  is  formed  from  the  scalar 
Green's  function  by  multiplication  by  Y“2  (which  con£a^ns  s“2) 
and  spatial  derivatives  as  in  VV.  Thus  for  finite  | r— r • | (but 
not  zero)  the  dyadic  Green's  function  is  also  an  analytic  func- 
tion of  s,  except  at  s * 0.  However  in  going  to  static  limits 
the  s“2  term  combines  with  other  terms  in  s to  avoid  an  object 
pole  at  s = 0 as  discussed  before.  Similarly  terms  like  VG  etc. 
appear  in  the  kernels  and  the  same  analyticity  properties  with 
respect  to  s apply. 

^ Now  in  forming_^the  elements  of  the  matrix  (gn,m)  for  tn  1* 
rm  the  coordinates  rn  and  rm  are  used  in  the  various  Green's 
function  type  terms.  In  such  cases  then  the  matrix  clement  must 
be  an  analytic  function  of  s except  possibly  for  poles  which 
don't  bother  us.  The  spatial  derivatives  can  be  replaced  by  fi- 
nite difference  operators  based  on  the  spacing  between  nearby 
jones^in  the  zoning  system  defined  to  segment  the  body.  For 
rg  * rh  (including  the  case  g * h)  the  problem  is  basically  a 
static  one.  The  scalar  Green's  function  can  be  written  as 

G(s,  |r-r'|)  - — p™  - h + I + °<s3>  <3-8> 

4TT|r-r'  I 

Taking  th§  fir^t  few  terms  as  needed  then  the  matrix  elements 
gn,m  for  r^  * rm  (including  gn,n)  can  be  defined  and  they  too 
are  analytic  functions  of  s.  Thus  the  matrix  (gn,m)  can  be  de- 
fined as  an  analytic  function  of  s except  for  possible  poles  of 
finite  order.  Note  that  finite  zone  size  is  important  here  be- 
cause zone  size  contributes  to  the  matrix  elements  and  vector 
components  and  we  want  to  avoid  infinite  values  for  these. 

One  should  be  careful  of  various  approximate  formulations 
so  that  the  resulting  matrix  elements  are  analytic  in  s without 
branch  cuts.  For  example  it  is  common  in  various  thin  wire  type 
formulations  to  obtain  terms  containing  functions  like  ln(yw) 
where  here  w is  the  wire  radius.  Such  terms  result  from  inte- 
grating the  Green's  function  to  remove  a coordinate  such  as  the 
azimuthal  angle  around  the  wire,  perhaps  approximating  the  re- 
sult and  then  integrating  say  along  the  wire  in  the  process  of 
finding  an  asymptotic  form  for  small  wire  radius.  Such  thin 
wire  formulations  are  important  and  perhaps  even  useful  in  de- 
veloping singularity  expansions  but  for  finite  size  bodies  they 
do  have  this  limitation  which  needs  to  be  recognized.  For  our 
present  purposes  we  do  not  use  such  types  of  formulations. 

The  cofactor  matrix  D for  (gn,m)  is 


D = 


<dn 

n,m 


( (-1) n+mdet  G J 
v n ,m' 


(3.9) 
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where  Gn/Itl  is  the  N - 1 by  N - 1 matrix  formed  by  deleting  the 
nth  row  and  mth  column  from  (gn,m)  * The  determinant  of  (gn,m) 
we  indicate  by 

4 5 aet<Vm>  <3-10) 


Note  that  the  transpose  is 


(3.11) 


The  solution  of  our  matrix  equation  can  then  be  written 

<VS>>  - <fm,n(s))  • (in(s))  (3-12> 

where 


1, . 


<fm,n>  “ r<dm,n> 


d „ - (-l)n4®det  G 
m,n  n,m 


(3.13) 


Now  since  the  gn/m  ate  analytic  functions  of  s except  pos- 
sibly for  poles,  then  the  same  is  true  of  both  the  dro/n  and  A 
since  the  determinant  is  a polynomial  function  of  the  matrix 
elements.  Therefore  since  the  In(s)  are  entire  functions  of  s 
and  the  fm.n(s)  are  analytic  except  for  poles,  then  the  Jm(s) 
must  also  be  analytic  functions  of  s except  for  poles.  Further- 
more any  poles  of  the  gn,m  elements  are  associated  with  poles  in 
the  medium  parameters  or  with  the  Green's  function  with  various 
operators  on  it.  Poles  associated  with  the  Green's  function  are 
at  b ■ 0 and  do  not  concern  us  because  they  give  no  resulting 
poles  with  non  zero  coefficients  based  on  physical  grounds  if 
the  media  are  passive.  Poles  in  the  medium  parameters  o,  e,  y 
can  be  troublesome  so  for  our  present  discussion  we  restrict 
ourselves  to  the  case  of  no  such  medium  poles.  Such  polee  in 
passive  medium  parameters  could,  however,  be  included  in  the  ex- 
pansion technique  if  desired.  With  this  restriction  then  the 
poles  in  (Jm(s))  are  the  zeros  of  A and  the  order  of  each  pole 
is  less  than  or  equal  to  the  order  of  the  corresponding  zero  of 
A. 


The  order  of  the  zero  at  sa  must  be  finite  since  if  A and 
all  its  derivatives  are  zero  at  Sq  then  as  an  analytic  function 
A must  be  identically  zero  implying  no  unique  solution  for  (3m) 
violating  the  uniqueness  theorem  for  the  solution  of  Maxwell's 
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equations.  Furthermore  the  zeros  of  A in  the  finite  s plane 
must  be  isolated,  i.e.  only  a finite  number  are  allowed  in  a fi- 
nite region  of  the  s plane.  This  again  is  a property  of  ana- 
lytic functions.  As  an  example  suppose  that  there  were  a line 
of  zeros  for  A.  Then  by  analytic  continuation  A would  be  zero 
in  the  entire  s plane,  again  not  allowed. 

This  result  is  of  course  of  fundamental  importance  to  the 
singularity  expansion  for  finite  objects.  It  relies  on  restric- 
tions on  the  form  of  the  medium  parameters  and  the  convergence 
of  the  particular  matrix-vector  formulation  of  various  possible 
integral  equations.  The  response  has  only  pole  singularities  of 
finite  order  and  isolated  in  the  finite  s plane.  Since  we  also 
restrict  the  object  to  be  passive  then  these  poles  cannot  be  in 
the  right  half  plane.  Furthermore  as  discussed  previously  (and 
in  appendix  A)  any  poles  on  the  iu>  axis  must  have  zero  coupling 
coefficients  and  so  we  can  state  for  our  purposes  for  passive 
objects  that  all  poles  in  the  finite  s plane  have  negative  real 
part  (fta  < 0) . Stated  another  way  how  car.  the  (Jm)  have  branches 
(removing  multiple  values)  when  the  elements  gn.m  are  chosen  to 
have  no  branches  because  they  are  single  valued? 

Another  way  to  view  this  question  of  only  poles  for  the 
singularities  in  the  finite  s plane  concerns  the  numerical  re- 
sults per  se.  Matrix-vector  formulation  of  various  integral 
equations  has  been  common  for  a long  time  and  has  been  ured  to 
obtain  many  accurate  solutions  for  finite  size  objects.  For 
formulations  which  have  used  analytic  matrix  elements  and  inci- 
dent vector  components  (except  for  poles  in  the  finite  s plane) 
the  numerical  results  obtained  (usually  on  the  iw  axis  and  some- 
times converted  to  time  domain)  must  be  representable  by  a sin- 
gularity expansion  which  has  only  poles  in  the  finite  s plane  by 
our  previous  discussion.  Thus  in  cases  where  accurate  results 
have  been  obtained  such  results  are  accurately  representable  by 
such  a singularity  expansion.  The  accuracy  of  such  a singular- 
ity expansion  for  finite  size  objects  is  t.ius  related  to  the  ac- 
curacy of  the  corresponding  matrix-vector  representation  of  the 
integral  equation  from  which  it  is  derived. 

Knowing  that  the  inverse  matrix  (fwi,n)  ^as  only  poles  in 
the  finite  s plane  this  still  leaves  open  the  question  of  sin- 
gularities at  infinity  associated  with  an  entire  function  such 
as  e~&to.  one  might  argue  that  this  matrix  is  like  an  admit- 
tance matrix  for  a large  circuit  composed  of  inductors,  resis- 
tors, and  capacitors  and  like  such  cases  might  have  only  ratios 
of  polynomials  in  s for  matrix  elements.  However  for  large  N 
the  degree  of  such  polynomials  could  get  arbitrarily  large. 

Thus  for  completeness  at  the  moment  let  us  write 


(f  (s)) 
m,  n 


[E 


-n_ 


(s-sa) 


(fm  n) 

m,  n 


(f  _ (s) ) 
m , n 


(3.14) 
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where  the  (im,n)a  are  constant  matrices  which  are  the  coeffici- 
ents of  each  pole  of  order  na  > 1 and  (fm,n(s))e  is  an  entire 
matrix  function  of  s.  An  interesting  question  for  further  in- ‘ 
vestigation  concern*;  the  properties  of  (fm,n(s))e  for  some  ma- 
trix size  N by  N and  in  the  limit  of  large  N;  perhaps  it  is  zero 
for  many  cases  of  interest. 

Written  in  matrix-vector  form  the  current  density  expansion 
as  in  equation  2.20  for  the  delta  function  response  with  polar- 
ization p is  then 


(Jm(s))  -2X(ei,s)<v, 

P a 


-n„ 
(S-S)  ° 

m a ° 


+ (Wm(eirs))  (3.15) 


At  each  sa  we  can  expand  each  na  in  a power  series  since  each  rja 
is  an  entire  function.  The  successive  terms  in  such  an  expan- 
sion produce  poles  of  order  n<j  - 1,  na  - 2,  etc.  associated  with 
each  na  2.  Note  that  for  eac£  sa  there  may  be  several  n0 
while  each  na  can  have  several  va  (degeneracy)  each  of  which  has 
an  na.  Expanding  the  inverse  matrix  (fm,n(s))  as  in  equation 
3.14,  substituting  this  in  equation  3.12,  and  multiplying  by 
(s-sa)N«  we  then  obtain  for  each  polarization 


E S0 <«!>*»> (vm) 

no  No 

(all  modes) 


a 


(f  n) 
m,n  „ 


• {1nieVBa)\ 

nc’Ma 


(3.15) 


where  Na  is  the  maximum  pole  order  n«  at  the  sa  of  interest. 

Note  that  the  number  of  independent  va  for  this  case  need  not  be 
infinite  because  the  coefficient  matrices  (fm,n)a  come  from  the 
inversion  of^ (gn,m(s) ) near  the  sa  where  A * 0.  Similarly  we 
can  expand  (In)  near  sa  to  obtain  the  coefficients  of  the  next 
higher  order  pole  contribution  (order  Na  - 1)  as 


E "°(Sl'8“,(Va  ♦ E [is  «.«•!->]  «y 


a a 


a a 


(f_,  n)  * (r(evBj)  + (f  n) 

m,n_  n l a m,n_ 

P 

a a 


n »N 
a a 


• aT[(Jn«cl's>>  ] 

P^S=S. 


(3.17) 
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Knowing  the  v0  and  the  choice  for  ?ia(ei,s)  then  the  terms  for 
na  “ Na  can  be  calculated  to  leave  this  last  equation  as  one  for 
the  na  » Na  - 1 terms.  This  can  be  repeated  to  obtain  all  the 
terms  for  the  poles  of  all  orders  at  sa.  If  the  va  are  degener- 
ate for  a given  na  and  are  constructed  to  give  an  orthogonal  set 
for  this  case  in  the  sense 


<v 


o' 


<v 


for  a ^ a1 


(3.18) 


then  the  coupling  coefficients  have  the  relation  for  no  as 

[<V0-  ( Vj'1* Va  • <*■,»>,  * <intVs«>> <3-191 

For  purposes  of  calculation  with  the  inverse  matrix  suppose 
A has  a zero  at  sa  of  ord  ;7a,  the  maximum  pole  order  of 
(fm,n(s))  at  sa.  Say  that  near  sa 


N«  N +1 

A(s)  * <s-sa)  An  + (s-sa)  An  +i  + •• 

a a 


1 -N„  “N  +1 

TOT  * <S"Sa)  6-N  + (s~sa)  fi-N  +1  + 

a a 


(3.20) 


with  the  coefficients  related  as 


-N 


N 


-v1 


v+1 

a 


N_ 


(3.21) 


and  so  on  by  inverting  the  series  for  A(s)  to  one  for  A~l(s) . 
Then  the  (fm,n)ct  can  b®  calculated  from 


(£m,n<»»  n<»» 


(3.22) 
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( 


I*a“na 


(f  „) 
m,n  „ 


n_ 


^ j 5 M If  (d_  _) 

-Na+£  m,n  , 


(d  ) 
ni/  n 


1 [d* 

*157 


<dm  „(8>> 

m,n 


s*s_ 


(3.23) 


and  so  on  by  collecting  terms  in  expanding  the  series  for  A(s) 
around  sa  into  a series  for  1/A(s).  Note  that  there  is  a dif- 
ferent set  of  (djn,n)£  coefficient  matrices  at  each  s0  so  they 
can  also  use  the  general  a subscript  set. 

Instead  of  using  the  inverse  matrix  approach  we  can  set  up 
an  eigenvector  equation,  factoring  out  ths  coupling  coefficient, 
by  first  expanding  (gn,m(s>)  in  a Taylor  series  around  sQ  as 


00 


(9n,m(s)) 

* 2 (s-s  )*  (g_  ) 

* a n ,m  . 

Il'[fpr(9n,m(8,) 

(3.24) 

(g„  m)  “ 

n,m  £ 

f 

Note  again  that  there  is  a separate  set  of  coefficient  matrices 
for  each  sa  so  the  (gn(m)i  can  also  use  the  a set.  Substitute 
this  series  together  with  equation  3.15  into  equation  3.1  to  ob- 
tain 


£*0 


V 


<V»> 


2XCV»>(y,  <B'sa> 
{ a Ot 


-n 


(Wm(ei,s)) 


(3.25) 


Since  the  right  side  has  no  poles  at  sa  then  collect  the  coef- 
ficients of  each  pole  on  the  left  side  and  equate  them  to  zero 
beginning  with  the  highest  order  pole  at  sa  for  na  » Na  as 


(gn  m) 
n,m  , 


(v) 
in  ^ 


(0)  - (g  _ ( s ) ) • (V) 
n,m  a n 


n*N 
a a 


(3.26) 
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where  the  coupling  coefficient  is  assumed  non  zero  and  factored 
out  and  (0)  is  an  N component  vector  with  all  zero  elements. 
There  may  be  several  natural  modes  for  na  * Na  perhaps  associ- 
ated with  symmetry.  If  there  is  mode  degeneracy  for  na  * Na 
then  construct  an  orthogonal  set  which  spans  the  space  of  such 
modes  so  that  the  property  in  equation  3.18  applies.  Going  to 
the  na  * N<*  - 1 term  we  have  (for  Na  >.  2) 


n?  • «V0  + Jr <v»>c 

a a a a a 


<V  - (0) 

a 


At  this  point  let  us  consider 
these  terms. 

Suppose  that  we  restrict 


n a(®l's)  H ca(®l) 
a A T(s)  a A 


(3.27) 

the  form  of  na  to  simplify  some  of 
fja  to  have  the  form 

(3.28) 


where  we  have  factored  out  the  complex  frequency  dependence  as  a 
common  factor  for  all  a.  This  common  factor  T(s)^is  specified 
to  be  an  entire  function  in  the  complex  s plane;  T(s)  may  even 
have  zeros  in  the  finite  s plane  if  we  wish  to  remove  some  poles 
before  making  a singularity  expansion.  The  new  or  modified  set 
of  coupling  coefficients  ca(ei)  are  specified  to  be  constants 
and  are  the  same  type  of  constant  coupling  coefficients  as  used 
in  section  2 as  for  example  in  equation  2.21.  One  choice  cf 
T(s)  is  clearly  esto  as  used  in  section  2. 

Now  multiply  the  matrix-vector  equation  3.1  by  T(s)  tc  ob- 
tain 


(gn#m(s))  • [f(s)  (3ro(s))  j - T(s)  (In(s))  (3^‘29) 

The  right  hand  side  is  still  an  entire  function.  Then  restrict 
for  our  present  analysis  the  coupling  Coefficients  to  have  the 
form  in  equation  3.28.  Solve  for  the  ca  and  divide  the  results 
by  T(s) . 

Equation  3.25  can  now  be  rewritten  as 
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l£«0  * a a ' 


T(s) (In(s)) 


(3.30) 


The  b derivatives  of  the  coupling  coefficients  as  in  equation 
3.27  now  do  not  appear.  Equation  3.26  can  then  be  generalized 
for  na  ■ 1,  2,  •••,  Na  as 


(v  ) 
m _ 


(0) 


n =N  -l 
a a 


(3.31) 


from  which  natural  modes  can  be  calculated  in  the  form  of  orth- 
ogonal natural  modes  constructed  for  each  for  which  there  is 
modal  degeneracy. 

Similarly  using  the  expansion  of  the  inverse  matrix  we  can 

write 


«(S)(5m(s).  - +<fn,n<s>>e)  * 


-n 


a a 


(3.32) 


Around  sa  make  the  Taylor  expansion 


v i 


T(s)(i  (S))  - y)(s-s  J1  (T) 
n a n . 

1*0  * 


(TJ 


1 [d* 


(T(S)  (ln(8)  ) 1 


s«s. 


(3.33) 


where  the  coefficient  vectors  (Tn) £ can  also  use  the  a index 
set.  Then  for  the  pole  at  sa  of  order  M we  have 
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N0-M 


Ecrt(e.)(v  ) - 2 n* 
- a 1 m„  M m*n  „ 


(3.34 


na=M 


n =M+l 
a 


from  which  we  can  obtain  both  the  natural  modes  (constructed  as 
orthogonal  set^.  and  the  coupling  coefficients.  Note  that  the 
dependence  on  e^  is  contained  in  the  vectors  (Tn)t.  If  we  let 
M ^ Na  then  we  can  relate  the  entir§  functions  as  well.  Alter- 
nately having  found  all  the  ca  and  va  these  can  be  subtracted 
from  equation  3.32  to  leave_the  entire  function  T(s) (Wm(ei,s) ) 
expressed  in  terms  of  T(s)(In(s))  together  with  the  inverse  ma- 
trix coefficients. 

With  the  natural  modes  known  equation  3.34  gives  an  expres 
sion  for  the  coupling  coefficients  fgr  a «*  o'  as 

°«’(^  * [rv”\-  • • 


|na»M+i 


(3.35 


The  coefficient  matrices  for  the  inverse  matrix  can  be  related 
to  the  Taylor  expansion  of  (gn,m(s))  through 


(g_  ^.(s))  • (f  , Js»  = (6  J = (f  _ » ( s ))  . (g  , (s» 

m ,m  nrm  nfiu  •'ro  fm 


(3.36 


(6  * !5D(s"8a,A^n  m.)  ! * I m>  + (*„.  Js»  ! 

n'm  Uto  a n'm  l)  (<jTl  m 'mq  m 'n  e) 

where  now  the  coefficient  matrices  (fm’^nJq  aPP*y  for  the  par- 
ticular 8a  of  interest;  the  remaining  terms  at  other  s0  can  be 
included  with  the  remainder  function  for  present  purposes. 
Equating  tv  resulting  singular  terms  to  zero  gives  for  the 

hi*rhfci-t  or  3.  pole 


(g_  ) • (f_,  _)  - (0  ) • (f  _, ) • (g,  ) 

n,m  m /in  „ nfni  nfm  „ m #m 

O N O 

a a 


(3.37 


where  (On#ii  *.s  an  N x N matrix  with  zero  elements.  Prom  this 
result  the  columns  (fixed  m)  o£  (fm'  m)No  must  be  eigenvectors 
of  (gn.m')o  and  the  rows  (fixed  n)  of  (gn/m’)o  must  be  eigenvec 
tors  of  the  transpose  of  (gn/m')o*  Also  the  rows  of  (fn/m')Na 


are  eigenvectors  of  the  transpose  of  (gm',m)o  and  the  columns  of 
<gm',m)o  are  eigenvectors  of  (fn#m')N0*  For  poles  of  order  n' 
in  equation  3.36  we  have  the  general  result 
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) • 
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,n 


) 

n'-Jt 


«L  J - Zj  -.> 
n'm  q-iin'  n'm  n'-i 


) 

1 


(3.38) 


where  1 < n'  <_  Na. 

Consider  the  special  but  practically  important  case  that 
the  pole  at  sa  is  simple  (Na  = 1) . Then  the  natural  modes  are 
the  eigenvectors  from 

(g_  .(si)  • (v  ) = (0)  (3.39) 

n,m  a m 


The  coupling  coefficients  are 
<=.<*!»  " [<Va  ’ * <fm,n 

’ K ’ ‘XJVV.  * 
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(3.40) 
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where 


(f  n)  * lim  ( s— s ) (f  ( s) ) 

m'n  1 8+s  a m'n 
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B+B_ 


8~sa 

ITsT 


(d  (s  )) 
m,n  a 


^7(am,n(8a,, 


(3.41) 


with  the  relation 

(g*  «•(*,>)•  <fm»  J = (°n  J = (fn  m.)  • (gm.  _(s  ))  (3.42) 

rtfin  c.  m*  ,m  ^ n,m  n/in  ^ in  a 

In  another  form  the  natural  modes  and  coupling  coefficients  are 
found  from 


( 
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2c«<ei)<vm)  “ (fm  n)  * (T  ) 

n7»l  ° 1 ma  m'n  1 n 


<fm  n>  • <*„<*«> >1 

m,n  ^ u n a 


which  for  the  case  of  no  modal  degeneracy  reduces  to 


ca(e,MvJ  - <f  ) • [T(SfJ  (i  (srt))l 

ct  jl  m ^ ^ cx  n cx 


(3.43) 


(3.44) 


There  is  an  alternate  way  of  calculating  the  coupling  co- 
efficients. Again  let  Na  “ 1 at  the  sa  pole.  From  equation 
3.42  the  columns  of  (fm,n)i  are  eigenvectors  of  (gn,m(sa))*  As- 
sume no  degeneracy  of  the  natural  modes  at  sa . Then  the  columns 
of  (fm.n)i  are  all  the  same  except  for  a scalar  factor  since 
(vnOa  is  unique  except  for  a scalar  factor.  Call  these  constants 
ym  so  we  can  write 


(f  _)  ■ (y, (v  ) , y-(v  ) , 
m,n  , i m ^ t ro  ' 
1 a a 


“n'V  W > i3'iS) 

a a 


Then  we  have 


(v  ) • (f  _)  » I(w  ) • (v  )](y  ) 

m a m,n  ^ m m nfl 


(3.46) 


defining  a new  N component  vector  so  that  (fm*n)i  can  be  written 
in  dyadic  form  as 


(f  _> 

ro,  n 


(V  ) 

m a n a 


(3.47) 


where  the  outer  vector  product  or  dyadic  product  is  used  here. 
Thus  we  also  have 


<*n> 


r<%>  • <vj(v 

L a a j a 


(3.48) 


Now  from  equation  3.42  the  rows  of  (fm,n)i  are  eigenvectors  of 
the  transpose  of  (gn,m(sQ))  so  that  we  have 
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(g_  _ (s) ) 
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(3.49) 
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Prom  equation  3.40  we  can  then  write  the  coupling  coefficient  as 


ca<«l>  ' <V  • 

a 


(3.50) 


giving  a simpler  form  of  the  result  in  terms  of  the  eigenvector 
of  the  transpose  of  (gn,m(sa)).  Note  that  since  we  have  assumed 
only  one  independent  (vnOct  there  is  only  one  independent  (um)a. 

Now  similar  statements  can  be  made  about  the  eigenvectors 
of  (fm,n)l  from  equation  3.42.  Specifically  the  columns  of 
(9n,m)  must  be  eigenvectors  of  (fm,n)l  and  the  rows  of  (gn,m) 
must  be  eigenvectors  of  the  transpose  of  (fm,n)l»  From  equation 
3.47  the  eigenvectors  of  (fm,n)l  are  orthogonal  to  (Mn)a  an<3  the 
eigenvectors  of  the  transpose  of  (fm,n)l  are  orthogonal  to  (vm)a> 

We  can  calculate  the  transpose  eigenvector  (yn)a  from  equa- 
tion 3.49  except  for  a scale  factor.  To  find  this  scale  factor 
we  can  use  equation  3.46  to  give 


[<H  • ‘ Vj  [<  va  • <vn>j  = <va-<w 


“5-l<V  • <Vn‘»«»  • <S„>  (3-51> 


If  we  have  found  a (Un)a  from  equation  3.49  then 
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Another  powerful  result  comes  from  taking  the  constant  term 
in  equation  3.36  as 

<6n  J - * Cf  . _)  + (g_  m.)  • 

n^iA  n,m  . in  *m « n»in  m n oi 

i i o e 
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(3.53) 


Dot  multiply  on  the  right  by  the  eigenvector  (vjn)0  of  (gm*  m)Q 
to  make  the  second  term  vanish  giving 
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(3.54) 


(v  ' 
n . 


,m^ 


This  defines  a matrix  which  maps  (\>m)a  back  into  itself.  Having 
some  (tin) a solving  equation  3.49  we  also  have  from  dot  multiply- 
ing on  the  left  of  equation  3.53  the  result 


(>0.) 


m' 


(fm,  m) 
m ,m 


(3.55) 


Thus  multiplying  equation  3.44  on  the  left  by  (wA)0  * (9n,m)l  * 
we  have  after  rearranging  terms 

ca(n>  ■ [(llA,0'  * <V”>’a] ‘ ‘V  <intso)  1 1 (3'56) 


so  that  we  can  construct  (yn)a  as 

(“n»0  - [<“A>0  • • ‘V„] 


giving  the  additional  result 


(u  ) * (g  ) • (v  ) **  1 

nJ0  '9n,m,1  m'0 


(3.57) 


(3.58) 


Thus  (yn) a can  be  explicitly  constructed  from  (gn#m)o 
through  its  eigenvector  and  the  eigenvector  of  its  transpose 
plus  another  matrix  (gn,m)i  which  also  comes  directly  from 
(9n,m(s)).  Thus  after  one  finds  so  es  a zero  of  A(s)  then  if 
this  zero  of  A(s)  is  simple  and  the  natural  mode  nondegenerate, 
both  the  natural  mode  and  the  coupling  coefficient  can  be  found 
from  (9n.m(a))  through  the  special  formulas  above  without  having 
to  calculate  the  coefficients  in  the  expansion  of  the  inverse 
matrix.  If  the  zero  of  A(s)  is  not  simple  but  of  second  order, 
third  order,  etc.  then  the  formulas  are  somewhat  more  complex. 
Thus  given  an  integral  equation  for  the  currents  on  the  object 
and  some  matrix-vector  representation  of  this  which  uses  ana- 
lytic elements  with  at  most  poles  and  assuming  the  approximate 
solution  converges  in  the  limit  of  large  N,  then  there  are  var- 
ious techniques  to  calculate  the  natural  frequencies,  natural 
modes,  and  coupling  coefficients.  Furthermore  the  coupling  co- 
efficients can  have  many  forms,  all  of  which  give  the  correct 
contributions  at  the  poles  but  can  have  various  behavior  as 
|s|  -►  «;  this  also  affects  the  form  of  the  entire  function  which 
is  an  additional  term  in  the  solution. 
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Let  us  now  list  some  of  the  alternate  forms  that  the  cou- 
pling coefficients  can  take.  While  all  forms  give  the  correct 
damped  sinusoid  parts  of  the  waveform  at  late  times  there  is 
still  the  question  of  convergence,  particularly  at  early  times. 
We  need  to  know  more  general  results  for  the  behavior  of  the  en- 
tire functions  (or  behavior  at  infinity)  for  various  types  of 
objects  as  they  enter  into  the  singularity  expansion.  Of  course 
convergence  to  the  correct  result  can  be  determined  in  part  by 
calculating  the  entire  functions.  One  can  observe  the  result 
for  the  matrix  inversion  times  the  incident  vector  for  particu- 
lar cases  of  ei  and  p in  frequency  and/or  time  domains  and  com- 
pare to  the  sum  of  the  first  several  poles  to  see  where  there  is 
or  is  not  convergence  and  what  it  takes  to  remedy  the  situation. 

Consider  then  several  types  of  coupling  coefficients. 

Type  Is  Factor  out  the  turn  on  time  of  the  object. 

This  type  of  coupling  coefficient  is  defined  by 


T(s)  = e ° (3.59) 

so  that  the  coupling  coefficients  factor  as  restricted  by  equa- 
tion 3.28  in  the  form 


na(ers)  = e ^(e^  (3.60) 

This  form  is  used  in  many  cases  in  section  2.  In  appendix  B the 
perfectly  conducting  sphere  is  shown  to  have  this  form  of  result 
with  no  entire  function  in  addition  to  the  pole  expansion. 
ThereTs  an  entire  function  e*sto  which  is  common  to  all  cou- 
pling coefficients.  However  whether  we  need  an  entire  function 
as  a separate  term  in  the  sum  is  not  at  all  clear  for  more  gen- 
eral objects.  In  time  domain  as  discussed  in  section  2 the  re- 
sulting pole  terms  go  over  to  damped  sinusoids  (including  cases 
of  no  oscillation)  for  simple  poles  with  powers  of  t appearing 
for  higher  order  poles.  For  t < to  there  are  no  currents  on  the 
objects  which  is  physically  correct.  At  t » to  all  modes  turn 
on  all  over  the  object.  Before  any  fields  can  reach  a particu- 
lar position  on  the  object  (associated  with  propagation  by  the 
shortest  possible  path  which  is  not  always  a straight  line)  the 
modes  with  the  transformed  entire  function  must  all  add  to  sero 
(if  the  sum  converges)  at  that  position,  so  this  is  a test  for 
this  type  of  coupling  coefficient  representation. 


Type  2:  Factor  out  the  time  the  incident  wave  function  first 
turns  on  at  the  observer  position. 

This  type  of  coupling  coefficient  is  defined  by 
Tf*n  = 

(3.61) 

fla<V8>  5 e_St  ca<®l) 


where 


(3.62) 


with  r*  as  the  particular  observer  position  on  the  object  where 
the  current,  charge,  etc.  is  to  be  calculated.  For  this  type  of 
definition  then  Ca(ei)  is  not  the  same  for  the  whole  object  be- 
cause the  definition  is  changed  to  apply  to  each  observer  posi- 
tion separately.  However  the  natural  frequencies  and  natural 
modes  are  still  the  same  for  the  whole  object.  For  t < t'  at 
the  observation  position  the  fields,  etc.  must  be  2ero  and  for 
most  positions  t'  > tQ  for  general  objects.  This  form  of  cou- 
pling coefficient  then  may  have  some  advantages  for  early  time 
calculations.  For  the^. numerical  problem  the  observation  posi- 
tion is  discrete,  say  r^.  Thus  we  would  use 


t' 


c 


(3.63) 


Note  the  tn  is  the  first  time  that  the  vector  component  In  in 
the  incident  vector  (In)  turns  on.  This  ca  is  as  easy  to  calcu- 
late as  the  one  in  equation  3.60  but  it  is  calculated  as  many 
times  as  there  are  observer  positions.  In  fact  for  t > t'  the 
two  forms  give  exactly  the  same  time  domain  waveform  for  simple 
poles.  For  t0  < t < tn  type  1 gives  a non  zero  waveform  for  the 
particular  pole  contribution.  Thus  we  can  calculate  the  type  1 
coefficient  but  just  wait  until  time  tn  to  turn  it  on. 

To  see  this  in  general  form  take  an  arbitrary  time  ta  and 
make  * 


T(s) 


(3.64) 


Then 
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(3.65) 


s°,taV  • 'W> 

a 

which  when  combined  with  the  simple  pole  (s  - sa)“l  gives  a time 
domain  form  as 


L_1[na(Si,s)(s-Sa)  n“]=eS“ta<V  * 


«~<t“t. ) 


u(t“ta) 


) 

! X 


s t 

(sa) ) e ° u(t-ta) 


(3.66) 


Thus  for  a simple  pole  the  chosen  starting  time  only  shifts  the 
turn  on  in  the  unit  step.  The  waveform  it  multiplies  stays  the 
same.  Therefore  for  t > ta  all  simple  pole  waveforms  are  the 
same  independent  of  ta.  For  higher  order  poles  a term 
(t  - ta)na“l  enters  in  but  this  can  be  expanded  to  leave  t001""1 
as  the  leading  power  invariant.  Furthermore  considering  all  the 
pole  contributions  at  sa  for  na  * 2,  3,  the  results  can  be 
manipulated  to  cancel  some  terms  but  we  do  not  go  into  this  here. 

Type  3:  Factor  out  the  time  that  resultant  fields  can  first 
exist  at  the  observer  position. 

One  can  choose  ta  as  in  equation  3.64  such  that^xt  is  the 
first  time  that  any  field  can  reach  the  observer  at  ra.  This 
can  be  longer  than  the  time  ta.  It  can  be  calculated  from  geo- 
metrical optics.  For  example  if  the  gbserver  is  on  the  surface 
of  a perfectly  conducting  object  and  ei  is  such  that  it  must 
point  through  the  object  to  reach  the  observer  then  the  wave 
must  come  around  the  object  and  arrive  at  a time  greater  than  t^. 
As  another  example  a dielectric  object  with  propagation  velocity 
less  than  c may  still  have  the  wave  propagating  through  the  ob- 
jecc  reach  an  observer  on  a "shadowed"  side  first;  this  time  is 
still  greater  than  t/*.  Also  for  example  the  expansion  of  the 
response  of  the  perfectly  conducting  sphere  need  not  be  started 
at  time  t « to  but  t0  can  be  replaced  by  a more  appropriate  ta 
in  all  the  step  functions. 

Clearly  if  one  choos®-*  ta  to  be  greater  than  the  time  a 
signal  first  begins  at  an  observer  then  there  must  be  an  entire 
function  contribution  foi  t < ca  because  the  pole  contributions 
would  all  be  zero  Such  a choice  would  not  seem  very  useful. 
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Type  4:  Expand  the  inverse  matrix  in  poles  but  leave  the  entire 
function  for  the  incident  wave  as  a coefficient. 

This  approach  is  attractive  for  the  case  that  the  expansion 
of  the  inverse  matrix  (fm,n(s))  as  in  equation  3.14  needs  no  ad- 
ditional entire  function  (fm,.n(s))e  in  the  expansion  or  at  least 
that  this  entire  function  is  known  and  preferably  has  a simple 
form.  Define  a time 


(3.67) 


While  the  time  t©  is  the  turn  on  time  this  new  time  ti  is  the 
time  when  the  incident  wave  has  reached  every  position  on  the 
object  neglecting  any  scattered  fields  from  the  object.  As  such 
ti  might  be  called  the  turn  off  time.  Suppose  we  write  the  in- 
cident vector  as 


(I  (e. ,s) ) 
n 1 
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") 


(3.68) 


so  that  the  time  delay  in  each  component  is  explicitly  displayed. 
Consider  a special  excitation  function  consisting  of  the  nth 
component  of  this  vector  being  as  abov§  but  all  other  components 
zero.  Then  we  would  calculate  an  ?ian(el'8?  with  e“s-n  as  a fac- 
tor. Repeat  this  for  all  n and  add  up  the  results  on  the  basis 
that  the  equations  are  linear  and  superposition  can  thus  be  ap- 
plied. We  would  then  calculate  our  coupling  coefficients  for 
the  case  of  simple  poles  as 


/ — st' \ 

na(Srs)  « (iin)  • (I^e^s))  » (pn)  • (^(e^e  nf  (3.69) 


For  higher  order  poles  derivatives  of  the  incident  vector  with 
respect  to  s also  come  in.  Now  in  the  time  domain  we  have  for  a 
simple  pole 

C<*i't)  - (Un>  • (In(«wt))  - (u  ) • (b(e1)6(t-t')>  (3.70) 

ax  n _ nx  n _ nx  n 
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The  pole  gives  a damped  sinusoid  e*ot  for  a simple  pole  and 
a factor  tna~1  for  higher  order  poles.  This  is  convoluted  with 
Ha  where  each  element  of  the  incident  vector  makes  its  contribu- 
tion at  a time  tn>  As  discussed  before  for  t < t#  where  t«  can 
be  arbitrary  the  time  domain  waveform,  at  least  for  sisqple  poles. 


is  the  same  regardless  of  ta.  This  applies  to  our  case  of  con- 
sidering each  component  of  the  incident  vector  as  separately  non 
zero  and  adding  up  the  results  at  the  end.  Thus  for  t < t0  the 
pole  contributions  are  zero.  For  to  < t < ti  the  convolution  as 
above  is  required  which  consists  of  turning  on  the  contribution 
of  each  component  of  the  incident  vector  at  time  tn  and  in  time 
domain  we  have 

L-1[natS1.sHs-sar1]  (bn(S1)u(t-t;).*°<t  tn)) 
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For  ti  < t this  result  goes  to  the  form  as  in  equation  3.66  giv- 
ing the  same  late  time  behavior  as  the  other  types  of  coupling 
coefficient  definition.  This  type  4 coupling  coefficient  gives 
a type  of  early  time  behavior  which  is  different  from  the  three 
previous  types  md  thus  gives  another  form  of  pole  contribution 
from  a convergence  viewpoint. 

We  have  shown  several  possible  ways  to  define  coupling  co- 
efficients. Any  of  these  or  combinations  of  them  can  be  used  to 
obtain  the  best  representation  in  the  sense  of  the  fewest  number 
of  terms  required  in  some  region  of  frequency  or  time  of  inter- 
est. There  are  clearly  many  other  cases  one  might  consider. 

The  ones  discussed  here  are  some  that  rather  directly  follow 
from  physical  considerations  and/or  give  simpler  results  consis- 
tent with  the  requirement  of  having  the  correct  coefficients  at 
the  poles.  Note  that  types  1 through  3 for  the  coupling  coeffi- 
cients can  all  be  calculated  on  a common  basis.  A ca  can  be 
found  for  any  choice  of  T(s)  in  the  form  of  a time  advance.  In 
particular  T(s)  can  be  1.  The  resulting  ca  then  applies  to 
every  position  on  the  object.  It  is  just  the  turn  on  time  in 
th3  unit  step  which  is  shifted  in  types  1 through  3 based  on  a 
physical  time  of  interest  either  for  the  whole  object  or  for  a 
local  observation  position  on  the  object.  Type  4 gives  more 
complicated  waveforms. 

The  considerations  about  the  expansions  in  this  section 
have  been  based  on  incident  delta  function  plane  waves.  They 
can  be  carried  over  to  other  incident  waveforms  directly  by  the 
techniques  outlined  in  section  2.  The  type  of  the  incident  wave 
may  influence  what  one  considers  as  the  best  form  for  the  cou- 
pling coefficients.  This  is  because  the  early  time  convergence 
of  the  expansion  will  be  affected  by  the  high  frequency  content 
of  the  incident  waveform.  As  an  example  suppose  we  have  an  in- 
cident s£ep  function  waveform.  Then  ire  can  invert  a term  of  the 
form  na(ei,s)/(sa(s-sa) ) into  the  time  domain  and  add  a static 
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term  or  we  can  invert  a term  of  the  form  na(ei,s)/[s(s-sa) ) into 
the  time  domain  which  has  a form  like  the  time  integral  of  the 
first  form.  This  should  give  different  convergence  characteris- 
tics for  early  times  and  thus  another  form  to  be  looked  at  in 
combination  with  types  1 through  4 coupling  coefficients. 

The  remaining  entire  function  in  the  general  expansion  is 
clearly  a subject  of  much  interest.  Depending  on  which  form  of 
coupling  coefficient  is  used  this  entire  function  clearly  has  a 
form  which  is  different  for  different  forms  of  coupling  coeffi- 
cients. If  ta  is  chosen  as  in  type  3 but  made  later  than  the 
time  a signal  first  begins  (resultant  field  signal)  at  the  ob- 
server then  such  an  additional  function  must  be  non  zero  to  give 
the  only  possible  fields  at  early  times  before  ta.  If  ta  is 
chosen  less  than  the  first  time  resultant  fields  reach  the  ob- 
server then  the  pole  terms  must  either  sum  to  zero  (if  conver- 
gent) or  have  their  sum  cancelled  by  such  a function.  If,  in 
specific  problems  being  calculated,  the  sum  does  not  go  to  zero 
for  such  early  times  then  this  other  function  must  be  non  zero 
for  such  a case.  An  optimum  choice  of  coupling  coefficients 
might  be  one  in  which  the  remaining  function  is  identically  zero. 
This  would  have  the  beginning  time  for  each  mode  no  later  than 
the  first  time  resultant  fields  reach  the  observer,  perhaps  even 
at  exactly  this  time.  The  "best"  coupling  coefficients  may  have 
more  complicated  forms  for  na  than  those  used  here.  Perhaps  us- 
ing geometrical  diffraction  theory  to  consider  asymptotic  forms 
for  ]s|  -*■  one  can  investigate  the  properties  of  the  remaining 
entire  function  and/or  impose  tighter  restrictions  on  the  "best" 
form  for  the  coupling  coefficients. 
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Let  us  now  consider  some  of  the  implications  of  the  singu- 
larity expansion  method  for  some  general  classes  of  objects. 

Host  of  the  attention  in  this  note  has  been  given  to  objects 
with  finite  linear  dimensions.  This  is  clearly  an  important 
class  of  objects  for  electromagnetic  interaction  and  scattering. 
The  perfectly  conducting  sphere  in  appendix  B is  an  example  of 
this  class  of  objects  and  can  be  used  to  suggest  various  general 
results  for  this  class  of  objects.  In  fact  some  of  the  results 
and  conjectures  discussed  in  other  sections  were  originally  sug- 
gested to  me  from  this  example . The  possible  factorisation  of 
the  coupling  coefficients  na(ei,s)  into  ea(ei)e“sto  for  per- 
fectly conducting  objects  of  finite  linear  dimensions  is  a good 
example  of  such  a conjecture.  The  fact  that  the  perfectly  con- 
ducting sphere  has  cnly  poles  in  the  singularity  expansion  while 
the  perfectly  conducting  infinite  length  circular  cylinder  has 
branch  cuts  in  its  singularity  expansion  (associated  with  the 
cylindrical  Hankel  functions  in  its  eigenfunction  expansion) 
suggested  that  the  expansion  only  in  poles  is  associated  with 
the  finite  dimensions  of  the  object. 

Tnere  are  other  finite  sized  perfectly  conducting  objects 
which  can  be  studied  analytically  for  their  singularity  expan- 
sions. Some  examples  might  be  the  prolate  spheroid,  oblate 
spheroid,  circular  disk,  etc.  Such  examples  have  less  symmetry 
than  the  sphere  and  one  would  then  expect  less  degeneracy  of  the 
natural  modes.  With  the  circular  disk  an  edge  would  be  intro- 
duced and  one  could  see  how  this  affected  the  natural  modes. 

For  finite  sired  perfectly  conducting  objects  with  an  axis 
of  symmetry  (objects  of  rotation)  one  can  base  a cylindrical  (1 1 , 
z)  coordinate  system  on  this  axis  and  decompose  the  natural 
modes  based  on  cos(m$')  and  sin(m$')  for  integer  m while  the 
coupling  coefficients  have  factors  cos  (m<j»i ) and  sin(m<f>i)  based 
on  one  cf  the  angles  of  incidence.  The  integral  equation  over 
the  surface  reduces  to  a one  dimensional  integral  equation  for 
each  m making  the  numerical  solution  simpler  and  the  indexing  of 
the  natural  modes  also  simpler.  Ors.  L.  Marin  and  R.  W„  Latham 
(private  communication)  are  already  putting  together  numerical 
techniques  to  handle  this  case. 

As  suggested  by  Prof.  C.  Taylor  (private  communication)  one 
can  also  look  at  thin  wire  approximations  to  simplify  tho  singu- 
larity expansion  analysis  for  such  cases  and  perhaps  obtain  ap- 
proximate analytic  expressions  for  the  natural  frequencies,  nat- 
ural modoo,  and  coupling  coefficients.  This  would  have  the  ad- 
vantage of  determining  the  approximate  values  of  these  quanti- 
ties and  suggesting  an  appropriate  indexing  system  for  the  fre- 
quencies, mod as , and  coefficients  for  some  rather  complex  object 
shapes,  such  aa  tliin  wire  models  of  aircraft  structures.  Then 


in  more  detailed  calculations  of  "fatter"  structures  one  can  use 
the  thin  wire  results  to  help  locate  all  the  natural  frequencies 
etc.  because  one  expects  the  singularity  expansion  quantities  to 
vary  somewhat  continuously  as  the  object  shape  and  dimensions 
are  changed.  Furthermore  one  could  develop  numerical  techniques 
in  which  the  approximate  thin  wire  results  are  used  as  first 
terms  in  an  iterative  solution  for  the  corresponding  fat  objects. 

While  we  have  been  viewing  the  singularity  expansion  method 
vom  the  viewpoint  of  the  interaction  or  scattering  problem 
there  is  no  reason  to  expect  this  method  to  be  limited  to  such 
problems.  This  method  is  fundamentally  based  on  the  expansion 
of  analytic  functions  of  the  complex  frequency  s in  terms  of 
their  singularities  in  the  complex  s plane.  For  example  antenna 
problems  in  transmission  and  reception  can  be  considered  from 
this  viewpoint.  Prof.  S.  W.  Lee  (private  communication)  has 
looked  at  some  features  of  the  cy ! indr ical  antenna  and  this 
method  appears  to  give  son.a  insight  here  and  can  even  be  used  to 
relate  this  antenna  problem  to  the  interaction  problem  for  a fi- 
nite length  perfectly  conducting  cylinder. 

The  results  for  finite  sized  perfectly  conducting  objects 
can  be  applied  to  other  types  of  objects  as  well.  Consider  an 
aperture  in  a perfectly  conducting  plane  as  shown  in  figure  4.1A. 
By  the  Babinet  principle  this  can  be  related  to  a complementary 
perfectly  conducting  disk. il* 12  Essentially  by  interchanging 
the  roles  of  the  electric  and  magnetic  fields  (rotating  the  po- 
larization) and  including  a plane  wave  term  for  the  reflection 
from  the  infinite  plane  the  solution  for  the  aperture  scattering 
can  be  found  from  that  for  the  disk  scattering.  Thus  one  can 
describe  the  deviation  of  the  currents,  fields,  etc,  from  the 
continuous  plane  case  by  means  of  the  natural  frequencies,  modes, 
and  coupling  coefficients  of  the  complementary  disk.  Then  these 
results  can  be  applied  to  define  natural  frequencies,  modes,  and 
coupling  coefficients;  the  modes  can  be  formulated  for  change  in 
fields  in  the  aperture  and/or  changes  in  surface  current  and 
surface  charge  densities  on  the  remainder  of  the  perfectly  con- 
ducting plane.  Thus  it  is  quite  possible  to  define  the  singular- 
ity expansion  for  the  lack  of  an  object,  i.e.  a hole,  at  least 
in  the  case  of  a perfectly  conducting  plane. 

Similar  conclusions  apply  to  a protrusion  on  a perfectly 
conducting  plane  as  shown  in  figure  4.1B.  This  follows  from 
image  considerations.  With  the  image  of  the  protrusion  included 
then  the  incident  field  can  be  split  into  symmetric  and  antisym- 
metric parts  with  respect  to  the  perfectly  conducting  plane. * 

The  interaction  of  each  part  with  the  equivalent  object  with  a 
symmetry  plane  in  place  of  the  perfectly  conducting  plane  can 
then  be  studied  separately.  However,  due  to  the  reflection  of 
the  incident  wave  at  the  perfectly  conducting  plane  only  an 
antisymmetric  field  distribution  can  exist  and  thus  contribute 
to  the  result.  Thus  the  change  in  the  fields,  currents,  etc. 
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THESE  TWO  PROBLEMS 
ARE  RELATED  THROUGH 
THE  BASINET  PRINCIPLE. 


THE  CHARGE  DUE 
TO  THE  APERTURE 
HAS  ONLY  SYMMETRIC 
NATURAL  MODES. 


EQUIVALENT  DISK 

A.  APERTURE  IN  A PERFECTLY  CONDUCTING  PLANE 


INCIDENT  FIELDS  ARE  ONLY 
ON  THIS  SIDE  OF  THE  PLANE. 


PROTRUSION 


THE  ELECTROMAGNETIC 
DISTRIBUTIONS  HAVE  ONLY 
AN  ANTISYMMETRIC  PART. 

THE  CHANGE  DUE  TO  THE 
PROTRUSION  HAS  ONLY 
ANTISYMMETRIC  NATURAL  MODES. 


\ I 

IMAGE  OF  PROTRUSION 


PERFECTLY  CONDUCTING 
INFINITE  PLANE 


B.  PROTRUSION  ON  A PERFECTLY  CONDUCTING  PLANE 


BOTH  SYMMETRIC 
ANO  ANTISYMMETRIC 
ELECTROMAGNETIC 


OBJECT 


ONLY  ANTISYMMETRIC 

ELECTROMAGNETIC 

DISTRIBUTIONS 


PERFECTLY  CONDUCTING 


TWO  MIRROR  OBJECTS  INTERACTING 
TO  SPLIT  NATURAL  FREQUENCIES. 

C.  TWO  MIRROR  OBJECTS 


IMAGE 




INFINITE  PLANE 


OBJECT  NEAR  PERFECTLY  CONDUCTING 
INFINITE  PLANE  WITH  SHIFTED 
NATURAL  FREQUENCIES. 


FIGURE  4.1  SOME  "OBJECTS"  WITH  A PLANE  OF  SYMMETRY 
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can  be  described  by  the  natural  frequencies  etc.  of  the  protru- 
sion with  its  image  but  only  the  antisymmetric  natural  modes  are 
needed.  This  is  in  contrast  to  the  aperture  problem  in  a per- 
fectly conducting  plane  where  only  the  symmetric  part  of  the  in- 
cident wave  interacts  with  the  aperture  making  the  aperture  nat- 
ural modes  also  symmetric. 

This  reasoning  can  be  carried  a little  further  in  consider- 
ing two  finite  size  objects  which  are  mirror  images  of  each 
other  with  respect  to  a symmetry  plane  as  shown  in  figure  4. 1C. 
Such  a pair  of  objects  can  be  considered  as  one  object  for  an 
expansion  in  natural  -'v  ' -a:. facies,  natural  modes,  and  coupling 
coefficients.  If  the  objects  are  far  apart  then  the  scattering 
from  one  will  not  be  vary  large  at  the  second  when  compared  to 
the  incident  fie 2c . For  large  separation  there  is  not  very  much 
interaction  between  the  two  and  we  can  treat  them  as  two  sepa- 
rate objects  with  the  same  natural  frequencies  and  the  same  nat- 
ural modes  except  that  the  modes  on  the  two  bodies  would  be  mir- 
ror images  of  one  another.  Now  natural  modes  can  be  multiplied 
by  any  scalar  merely  in  changing  their  normalization.  Thus  for 
the  two  objects  considered  as  one  we  can  define  natural  modes  as 
symmetric  and  antisymmetric  modes  by  taking  sums  and  differences 
of  the  mirror  modes.  Now  the  symmetric  and  antisymmetric  parts 
give  an  exact  division  of  the  natural  modes  on  an  object  with  a 
symmetry  plane;  there  is  no  interaction  between  the  two.  Even 
for  the  two  objects  close  together  this  is  the  case.  As  the  ob- 
jects are  brought  together  one  may  typically  expect  a splitting 
of  the  natural  frequencies  in  two,  one  with  the  symmetric  and 
the  other  with  the  antisymmetric  modes.  This  is  analogous  to 
the  energy  level  splitting  in  quantum  mechanics,  say  as  two 
identical  atoms  are  brought  together.  For  a single  finite  size 
object  near  a perfectly  conducting  plane  only  antisymmetric 
modes  contribute  and  so  only  natural  frequencies  associated  with 
the  antisymmetric  modes  are  present.  As  the  object  nears  the 
perfectly  conducting  plane  one  might  typically  expect  a shift  of 
the  natural  frequencies  and  not  a splitting  unless  some  symmetry 
in  the  object  is  destroyed  in  the  process  or  there  is  some  other 
degeneracy  in  the  natural  frequencies. 

These  results  with  perfectly  conducting  infinite  planes 
suggest  yet  further  results  might  foe  obtained  for  various  per- 
turbations on  infinitely  large  perfectly  conducting  objects. 

The  perturbation  is  regarded  as  the  "object"  and  what  is  calcu- 
lated is  the  change  in  the  electromagnetic  quantities  associated 
with  the  introduction  of  this  "object."  If  the  perturbation  is 
of  finite  size  then  we  might  expect  its  singularity  expansion  to 
comprise  of  natural  frequencies,  natural  modes,  and  coupling  co- 
efficients. However  more  development  is  needed  to  understand 
this  thoroughly.  One  might  even  extend  this  perturbation  con- 
cept to  perturbations  on  finite  sized  objects.  Particularly  if 
the  perturbation  is  small  compared  to  the  object  dimensions  then 
one  could  make  a singularity  expansion  for  the  object  and  use 
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the  results  to  define  the  incident  fields  for  the  perturbation 
which  might  in  turn  be  approximately  solved  by  another  singular- 
ity expansion. 

Of  course  there  is  the  question  of  the  singularity  expan- 
sion for  infinite  or  semi  infinite  objects.  We  know  that  the 
perfectly  conducting  circular  cylinder  has  branch  cut  contribu- 
tions in  its  singularity  expansion.  There  are  many  other  shapes 
besides  spheres  and  cylinders  which  can  be  treated  from  the 
viewpoint  of  eigenfunction  expansions. 13  This  can  be  used  to 
help  divide  up  the  terms  in  the  singularity  expansion  by  treat- 
ing each  term  in  the  eigenfunction  expansion  separately  for  its 
singularity  characteristics.  Prof.  Garbacz  has  developed  a 
method  for  calculating  eigenfunctions  associated  with  general 
geometries  for  lossless  objects.?  Perhaps  these  expansions  can 
be  used  by  studying  each  term  to  determine  its  singularity  ex- 
pansion and  thereby  aid  in  developing  or  even  indexing  the  terms 
in  a singularity  expansion  for  such  objects.  In  studying  infi- 
nite objects  such  as  general  cylinders  (say  irregular  but  of 
some  maximum  "radius")  or  semi  infinite  objects  such  as  general 
cones  (say  contained  within  some  maximum  "half  cone  angle”)  one 
may  find  some  general  properties  of  the  singularity  expansion 
associated  with  certain  features  of  the  general  geometries. 

This  in  turn  may  give  some  guidance  on  how  to  approach  the  sin- 
gularity expansions  (and  index  them)  for  specific  cases  of  such 
objects. 

Note  that  some  objects  which  are  finite  in  size  may  have 
properties  of  infinite  bodies  such  as  branch  cuts  in  the  singu- 
larity expansion.  For  example  take  a perfectly  conducting  ob- 
ject of  finite  size  located  between  two  infinite  parallel  per- 
fectly conducting  plates.  This  can  be  replaced  by  an  equivalent 
problem  involving  an  infinite  humber  of  images  extending  infi- 
nitely far  away.  This  is  basically  a segmented  infinite  object. 
For  the  case  of  a thin  wire  of  finite  length  between  parallel, 
plates  some  frequency  domain  results  (on  the  iw  axis  of  the  s 
plane)  exhibit  peculiar  step  and  slope  discontinuities.14  This 
may  be  associated  with  new  terms  such  as  branch  contributions  in 
the  singularity  expansion.  Perhaps  the  case  of  infinitely  re- 
peated objects  (and/or  images)  in  one,  two,  and  three  directions 
can  be  specially  treated  so  as  to  obtain  some  general  results 
for  the  singularity  expansions  for  such  problems.  Since  trans- 
lational symmetry  is  present  for  such  repeated  objects  perhaps 
group  theory  considerations  can  be  applied  to  obtain  general  re- 
sults for  this  type  of  problem. 

Other  important  classes  of  objects  involve  lossy  media  of 
infinite  size,  such  as  half  spaces  with  finite  non  zero  conduc- 
tivity. Objects  of  finite  size  may  be  clcse  enough  to  such 
media  to  affect  their  response  character is vies,  thereby  altering 
their  singularity  expansions.  Furthermore,  infinite  objects 
such  as  wires  may  be  in  proximity  to  semi  infinite  lossy  half 
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spaces;  this  could  introduce  yet  additional  features  in  the  sin- 
gularity expansion. 

Another  whole  class  of  problems  concerns  the  analysis  of 
experimental  data,  say  from  tests  using  EMP  simulators.  Using 
numerical  Laplace  transforms,  numerical  Hilbert  transforms,  etc. 
features  of  the  singularity  expansion  of  the  experimental  data 
can  be  found  within  the  limits  of  accuracy  of  the  experimental 
data  and  the  numerical  techniques.  It  would  seem  that  various 
approaches  to  this  problem  are  possible  depending  on  the  kind  of 
experimental  data  and  type  of  object  being  considered. 

Clearly  there  are  numerous  topics  in  the  singularity  expan- 
sion method  involving  classes  of  objects,  numerical  techniques, 
etc.  which  need  extensive  development.  In  this  note  we  have  for 
the  most  part  considered  finite  size  objects.  Even  for  this  im- 
portant though  limited  class  of  objects  much  needs  to  be  done. 
For  example  the  natural  modes  come  from  the  coefficients  of  the 
poles  in  the  expansion.  The  natural  modes  may  be  orthogonal 
over  the  volume  or  surface  of  the  object.  Such  is  the  case  for 
the  sphere  but  what  about  in  general?  Perhaps  the  topology  of 
the  object  can  be  used  to  help  index  the  singularity  expansion 
quantities  and  can  be  used  to  identify  whether  or  not  certain 
kinds  of  terms  are  present.  The  object  symmetries  have  much  in- 
fluence on  the  degeneracy  of  natural  frequencies  and  modes. 

Group  theory  should  then  oe  useful  in  understanding  the  degener- 
acy and  splitting  up  the  resulting  modes  as  well  as  indexing  the 
natural  frequencies,  natural  modes,  and  coupling  coefficients. 

A finite  size  object  need  not  be  perfectly  conducting. 
Suppose  that  it  is  composed  of  linear  passive  media.  Prof. 

C.  T.  Tai  has  suggested  (private  communication)  that  the  types 
of  theoretical  considerations  applied  in  circuit  theory  can  be 
applied  to  the  singvilarity  expansion  to  obtain  new  general  re- 
sults. Considerations  like  conservation  of  energy  are  important 
here  in  constraining  the  allowable  forms  of  the  solutions.  Per- 
haps some  properties  of  the  pole  expansion  can  be  deduced  such 
as  the  permissable  order  of  the  poles.  The  perfectly  conducting 
sphere  has  only  simple  poles.  It  seems  safe  to  conjecture  that 
this  is  true  of  all  finite  size  perfectly  conducting  objects;  as 
of  yet  I have  not  found  any  case  to  the  contrary . The  case  of 
lossy  objects  may  admit  more  general  pole  types.  The  properties 
of  the  coupling  coefficients  also  need  investigation  for  these 
more  general  finite  size  objects.  Whether  or  not  the  linear 
passive  media  are  also  reciprocal  should  also  have  important  im- 
pact on  the  properties  of  the  singularity  expansion. 

Numerical  calculations  for  objects  such  as  finite  length 
perfectly  conducting  cylinders  can  be  used  to  test  various  nu- 
merical techniques  for  calculating  the  singularity  expansion 
quantities  and  suggest  improvement  on  them.  Dr.  F.  Tesche  is 
already  doing  such  calculations  showing  some  cases  of  very  rapid 


63 


convergence  for  the  step  function  response.  Furthermore  such 
calculations  can  test  conjectures  of  a general  nature  and  if 
calculations  are  performed  for  various  different  objects  in  the 
class  of  interest  and  the  conjecture  proves  correct  in  all  cases 
then  one  has  high  confidence  in  the  general  validity  of  the  con- 
jecture. Consideration  of  various  example  objects  has  been  very 
useful  to  me  in  suggesting  new  results  and  testing  old  conjec- 
tures. It  would  seem  that  one  important  driving  force  in  the 
future  development  of  this  technique  will  be  the  calculation  of 
the  object  response  of  various  important  objects  of  practical 
interest. 

An  interesting  and  important  question  for  finite  objects 
concerns  the  uniqueness  of  the  form  of  the  singularity  expansion. 
Clearly  the  natural  frequencies  and  modes  for  finite  Js0|  are 
well  defined , but  th$re  are  rjme  possible  alternatives  in  defin- 
ing na(ei,s)  and  ca(ei),  the  forms  of  the  coupling  coefficients. 
Specifically  fja  can  be  an  entire  function  of  s (no  singularities 
in  the  finite  s plane) . This  gives  some  flexibility  in  choosing 
the  form  of  the  coupling  coefficients.  Of  course  the  choice  is 
not  completely  arbitrary.  The  final  resulting  current  density , 
charge  density,  etc.  are  unique  quantities  and  all  exact  repre- 
sentations of  them  must  amount  to  the  same  thing.  The  individ- 
ual terms  in  a series  expansion  can  be  altered  as  long  as  the 
sum  remains  the  same.  This  then  raises  the  question  of  what  is 
the  "best-  form JLn  which  to  express  the  coupling  coefficients. 

A form  like  cu(ei)e~8to  clearly  has  much  to  recommend  it  for  its 
simplicity.  However  there  may  be  other  terms  needed  if  such  a 
form  is  used.  This  problem  is  associated  with  the  time  during 
which  the  incident  delta  function  wave  is  sweeping  over  the 
body.  Since  the  object  response  must  be  zero  before  a field  can 
reach  any  particular  point  on  the  object  (with  this  time  calcul- 
able from  geometrical  diffraction  theory  considerations)  and 
since  to  is  the  time  the  first  point  on  the  object  is  excited, 
then  all  the  terms  in  the  expansion  must  sum  to  zero  (if  the  sum 
converges)  for  times  between  to  and  the  time  that  an  excitation 
can  reach  a point  of  interest  on  the  object.  In  order  to  best 
define  the  expansion  for  early  times  so  as  to  obtain  the  most 
rapid  convergence  then  some  other  definition  of  the  coupling  co- 
efficients may  be  appropriate.  The  possible  alternative  forms 
of  the  coupling  coefficients  seems  to  me  to  be  an  issue  of 
fundamental  importance  in  the  whole  theory  and  practical  utility 
of  the  singularity  expansion  for  finite  size  objects,  especially 
for  early  times.  For  infinite  bodies  the  terms  analogous  to 
these  coupling  coefficients  may  also  have  similar  questions  as- 
sociated with  them.  Much  research  is  needed  then  on  coupling 
coefficient  representation  both  in  terms  of  general  considera- 
tions and  specific  examples. 

In  past  years  there  has  been  some  consideration  of  the  nat- 
ural frequencies  of  some  simple  objects  and  to  some  extent  the 
natural  modes  have  also  been  investigated  for  such  objects.  Let 


us  mention  some  examples.  Thin  wire  natural  frequencies  have 
received  some  attention. 15  Natural  frequencies  have  been  dis- 
cussed in  the  context  of  antenna  resonances. 15  Natural  frequen- 
cies and  modes  of  a sphere  have  been  discussed; 17  this  forms  a 
starting  point  for  our  discussion  of  the  singularity  expansion 
of  the  perfectly  conducting  sphere  in  appendix  B.  Prolate 
spheroidal  geometry  has  alsc  been  considered  for  natural  fre- 
quencies and  modes. I8  Since  the  perfectly  conducting  sphere  has 
shown  so  many  interesting  results  it  would  seem  a good  idea  to 
look  at  the  prolate  and  oblate  spheroids  to  see  to  what  extent 
the  general  form  of  the  results  carries  over  to  these  geometries. 
For  example  spheroidal  geometries  can  allow  one  to  look  at  some 
analytic  results  for  forms  of  the  coupling  coefficients.  These 
and  other  investigations!  even  though  limited  in  sane  respects, 
at  least  solve  some  portions  of  terms  in  the  singularity  expan- 
sion for  some  objects.  As  such  they  can  shed  some  light  on  some 
details  of  the  singularity  expansion  for  such  objects  and  give  a 
start  for  obtaining  the  full  singularity  expansions.  They  also 
give  some  guidance  about  what  problems  can  be  profitably  consid- 
ered for, examples  to  develop  general  results  for  the  solution 
representations  in  the  singularity  expansion  method. 


V.  Summary 


This  note  is  intended  to  introduce  a new  way  of  looking  at 
many  kinds  of  EMP  interaction  problems,  although  it  has  bearing 
on  scattering  problems  as  well.  From  an  EMP  time  domain  inter- 
action viewpoint  this  approach  has  the  potential  for  directly 
calculating  t a amplitudes,  frequencies,  damping  constants,  and 
phunea  of  the  dampW  sinusoidal  oscillations  that  are  commonly 
seen  as  major  portions  of  interaction  waveforms  on  systems  under 
test.  The  idea  is  to  then  construct  a large  portion  or  even  all 
of  such  waveforms  as  a sum  of  such  damped  sinusoids. 

The  general  technique  can  aptly  be  called  the  singularity 
expansion  method  because  it  is  based  on  representing  the  func- 
tions of  the  complex  frequency  s in  terms  of  their  singularities 
in  the  complex  s plane.  In  the  time  domain  the  individual  terms 
are  the  inverse  Laplace  transforms  of  the  singularity  terms. 

While  for  general  objects  we  can  expect  branch  cut  contributions 
the  results  for  finite  size  objects  using  well  behaved  media  in- 
clude only  poles  for  the  singularities  in  the  finite  s plane. 

This  simplifies  the  form  of  the  terms  considerably  and  allows 
one  to  factor  the  terms  into  natural  frequencies,  natural  modes, 
and  coupling  coefficients.  The  natural  frequencies  and  modes 
are  independent  of  the  incident  wave  parameters  while  the  inci- 
dent wave  parameters  enter  into  the  coupling  coefficients  for 
the  delta  function  response.  The  incident  wave  can  also  have 
singularities  in  the  finite  s plane  but  these  can  be  separated 
out  so  that  the  response  can  be  generally  written  as  the  sum  of 
an  object  part  and  a waveform  part. 

There  are  various  matrix  techniques  for  solving  integral 
equations  numerically.  In  this  note  we  have  considered  these 
from  a general  viewpoint,  not  specifying  which  integral  equation 
is  being  approximated.  This  shows  some  general  ways  to  calcu- 
late natural  frequencies,  natural  modes,  and  coupling  coeffici- 
ents. The  actual  numerical  procedures  that  cne  could  use  are 
numerous  and  need  to  be  considered  for  various  problems  to  de- 
termine the  most  efficient  and  accurate  techniques. 

There  are  various  theoretical  problems  associated  with  the 
convergence  of  the  matrix  representations  which  need  to  be  con- 
sidered for  the  integral  equations  for  finite  objects.  Prefer- 
ably continuous  operators  over  the  body  geometry  can  be  developed 
to  analytically  represent  the  terms  in  the  singularity  expansion. 
The  question  of  the  singularities  at  infinity  or  additional  en- 
tire functions  for  finite  objects  needs  treatment.,  Of  course 
the  completeness  of  the  singularity  expansion  with  some  allow- 
able chosen  form  for  the  coupling  coefficients  can  be  readily 
checked  for  any  given  boundary  value  problem  by  comparison  to 
the  solution  by  standard  numerical  techniques.  This  method  can 
be  used  to  determine  better  forms  of  coupling  coefficients. 
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The^e  has  been  some  work  done  in  the  paBt  on  the  natural 
frequencies  of  objects  and  less  work  done  on  natural  modes. 

This  work  can  serve  as  useful  starting  points  for  the  singular- 
ity expansions  of  some  classical  geometries  which  can  be  used 
for  test  problems.  There  are  the  essentially  new  questions  of 
the  pole  order  and  coupling  coefficients  for  finite  size  bodies. 
This  note  has  included  a common  object  of  past  investigations., 
the  perfectly  conducting  sphere.  There  are  analytic  forms  for- 
coupling  coefficients  and  there  are  only  simple  poles  in  the  fi- 
nite s plane  with  no  additional  entire  function  required  for  the 
delta  function  response.  Investigation  of  other  common  objects 
considered  previously  should  also  give  some  valuable  insight 
into  appropriate  forms  for  coupling  coefficients  and  questions 
such  as  pole  order. 

I hope  that  this  note  has  given  the  reader  some  insight  into 
what  the  singularity  expansion  method  is  all  about,  particularly 
with  regard  to  finite  size  objects.  It  appears  to  be  quite 
powerful  for  some  kinds  of  EMP  interaction  problems.  Several 
investigators  are.  already  performing  sor  2 studies  of  both  gen- 
eral problems  and  specific  examples  using  this  type  of  expansion. 
I would  hope  then  that  the  near  future  will  see  some  significant 
additions  both  to  the  theory  of  the  method  and  problems  solved 
using  the  method. 


"Off  with  her  head  I"  the  Queen  shouted  at  the  top  of  her 
voice.  Nobody  moved. 

"Who  cares  for  you?"  said  Alice  (she  had  grown  to  her  full 
size  by  this  time).  "You're  nothing  but  a pack  of  cards!" 

At  this  the  whole  pack  rose  up  into  the  air,  and  came  fly- 
ing down  upon  her ; ... 

"Wake  up,  Alice  dear!"  said  her  sister.  "Why,  what  a long 
sleep  you've  had!" 

"Oh,  I've  had  such  a curious  dream!"  said  Alice,  and  she 
told  her  sister,  as  well  as  she  could  remember  them,  all  these 
strange  Adventures  of  hers  that  you  have  just  been  reading 
about;  . .». 


(Lewis  Carroll,  Alice  in  Wonderland) 


Appendix  A:  The  Natural  Frequency  s » 0 

As  mentioned  in  section  2 it  is  possible  for  finite  sized  " 
bodies  to  have  static  current  and/or  charge  distributions.  How- 
ever these  do  not  couple  to  the  incident  wave  and  thus  do  not 
start  or  stop  at  a time  such  as  to.  These  kind  of  static  solu- 
tions can  be  added  as  separate  terms  to  any  response  of  the  ob- 
ject to  the  incident  wave.  This  is  not  the  same  term  as  the 
static  response  of  the  Object  as  in  equations  2.46  and  2.52. 

The  natural  modes  for  s = 0 are  similar  to  others  on  the  iu 
axis  in  that  they  correspond  to  no  power  going  to  infinite 
radius.  For  static  current  and  charge  distributions  the  fields 
at  large  r for  di.pole  and  higher  terms  decay  like  r”3  (no  radi- 
ated power) ; the  r~l  electric  monopole  term  is  a radial  electric 
field  (no  radiated  power) . For  other  poles  on  the  iw  axis  not 
at  s » 0 any  radiation  field  like  r~l  at  infinity  would  consti- 
tute radiated  power  and  thus  damping  the  mode  (i.e.  making 
ft  < 0) . Thus  such  modes  when  expanded  over  some  sphere  contain- 
ing the  object  in  terms  of  divergenceless  spherical  wave  func- 
tions must  give  no  terms  behaving  like  r~l  at  infinity,  whereas 
all  such  fur.  tions  do  for  s j*  0 on  the  iw  axi-  Thus  the  fields 
for  such  modes  are  contained  in  some  volume  of  finite  dimensions 
and  might  logically  be  called  the  cavity  modes. 

What  are  some  oZ  the  characteristics  of  these  static  natu- 
ral modes?  Since  we  are  dealing  with  the  case  of  s * 0 the  wave 
equation  reduces  to  the  Laplace  equation  and  the  electric  and 
magnetic  fields  are  decoupled.  Thus  we  first  distinguish  be- 
tween electrostatic  natural  modes  and  magnetostatic  natural 
modes.  As  shown  in  figure  ALA  there  are  various  types  of  ex- 
amples as  in  1 a perfectly  conducting  object  with  net  charge  Qj.? 
thic  gives  rise  to  an  electrostatic  natural  mode  surface  charge 
density  Pej.  A3  in  2 an  insulating  dielectric  can  have  a net 
charge  Q2  giving  rise  to  the  same  kind  of  field  at  large  r as  in 
1;  the  charge  distribution  throughout  the  body  is  not  constrained 
by  the  surface  shape  as  in  1 and  static  electric  fields  in  the 
body  are  possible  so  that  the  distribution  is  arbitrary  to  some 
extent.  As  in  3 conductors  (perfect  or  imperfect)  can  be  com- 
bined with  insulating  dielectrics  to  allow  not  only  a net  charge 
but  allow  a volume  charge  distribution  in  some  parts  but  not  in 
others.  Of  course  all  the  above  mentioned  cases  can  be  combined 
together,  say  as  laid  out  in  figure  A1A  as  multiple  objects  in  a 
volume  of  finite  dimensions  so  that  the  whole  ensemble  has  a net 
charge  Qi  + Q2  + • • • and  an  associated  volume  and  surface  charge 
distribution.  It  is  not  necessary . for  there  to  be  a net  charge, 
or  even  for  there  to  be  any  fields  for  large  r.  The  charged  ob- 
jects could  be  contained  in  a closed  conducting  shell  and  the 
net  charge  of  all  (including  the  shell)  made  zero;  the  natural 
electrostatic  mode  would  still  have  a non  trivial  charge  distri- 
bution. 
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In  figure  A1B  we  illustrate  some  magnetostatic  natural 
modes.  Excluding  magnetic  charge  from  Maxwell's  equations  wa 
have  no  magnetic  monopole  term  to  give  an  r~l  magnetic  fiold  at 
large  r.  However,  we  can  do  next  best  and  get  a magnetic  dipole 
term  from  a closed  perfectly  conducting  loop  as  in  1.  If  the 
loop  is  perfectly  conducting  then  the  total  magnetic  flux 
through  the  loop  cannot  change  because  this  would  imply  an  elec- 
tric field  tangential  to  the  perfect  conductors  which  is  impos- 
sible by  hypothesis.  A perfectly  conducting  object  as  in  2 can 
have  a surface  current  density  3s2  when  immersed  in  a magnetic 
field  such  as  from  the  loop  1.  A permeable  object  as  in  3 can 
have  a magnetization  (with  an  equivalent  volume  current  density 
distribution)  induced  by  a magnetic  field  from  the  loop  1.  One 
might  also  consider  permanent  magnets  but  since  we  wish  to  remain 
with  linear  Maxwell's  equations  we  may  wish  to  exclude  such 
things. 

Practically  speaking  we  do  not  normally  deal  with  perfect 
conductors  but  they  still  make  a useful  idealization  for  many 
problems.  Thus  the  magnetostatic  natural  modes  still  are  useful 
concepts.  In  a practical  case  these  would  not  be  exactly  at 
s = 0 but  have  fi  slightly  negative.  Of  course  for  the  case  of 
superconductors  such  magnetostatic  modes  do  exist  and  have  been 
observed;  they  are  even  quantized.  This  leads  to  another  phe- 
nomenon in  which  the  magnetostatic  fields  are  excluded  from 
superconductors  except-  for  a thin  surface  layer.  This  is  analo- 
gous to  the  case  of  magnetostatic  fields  excluded  from  perfectly 
conducting  objects  as  discussed  above.  However,  there  may  be 
cases  of  highly  conducting  objects  for  which  it  is  useful  to 
think  of  a magnetostatic  mode  with  magnetic  field  penetrating 
what  is  thought  of  as  a perfect  conductor  for  purposes  at  hand. 

Having  considered  the  electrostatic  and  magnetostatic  natu- 
ral modes  there  is  no  reason  why  one  cannot  combine  them  and 
ha^e  both  associated  with  some  object  or  collection  of  objects 
contained  in  some  volume  of  finite  dimensions.  Note  that  the 
static  natural  modes  can  easily  be  degenerate.  Por  example  one 
can  change  total  charge  (Qi,  Q2»  etc.)  on  each  of  several  dis- 
crete conductors  as  well  as  change  the  currents  (Ii,  l2r  etc.) 
circulating  around  perfectly  conducting  loops  (or  equivalently 
separate  holes  through  perfectly  conducting  objects) . All  of 
these  apply  at  the  natural  frequency  u * 0.  Stated  briefly  the 
natural  mode3  for  s «*  0 are  any  electrostatic  and/or  magneto- 
static modes  which  are  not  associated  with  any  incident  field. 


Appendix  B . Example  to  Illustrate  the  Singularity  Expansion: 
TnePerl'ectly  Conduct JLng~5phere  * c 

At  an  aid  to  understanding  the  form  of  the  actual  singular- 
ity expansions  we  consider  an  exapple  chosen  so  that  the  various 
terms  in  the  expansion  may  be  more  readily  expressed  in  terms  of 
common  functions.  For  this  purpose  we  choose  the  perfectly  con- 
ducting sphere.  This  example  shews  all  simple  poles  in  the  ex- 
pansion and  poles  on  both  the  negative  ft  axis  and  in  conjugate 
pairs  with  ft  < 0;  the  internal  cavity  inodes  and  an  electrostatic 
mode  with  poles  on  the  iu>  axis  which  have  zero  coupling  coeffi- 
cients are  not  included.  The  natural  modes  are  degenerate  and 
we  choose  those  for  the  surface  current  and  charge  densities  to 
correspond  to  an  appropriate  set  of  spherical  harmonics.  The 
coupling  coefficients  have  a time  advance  which  factors  out 
leaving  dependence  only  on  direction  of  incidence  and  polarisa- 
tion. Again  this  example  is  directed  toward  explicitly  exhibit- 
ing the  form  of  the  singularity  expansion  and  the  kinds  of  gen- 
eral results  that  ensue.  For  numerical  purposes  the  singularity 
expansion  may  not  be  the  mose  useful  in  the  case  of  the  sphere. 
However  it  can  be  used  to  more  readily  suggest  general  results 
which  also  apply  to  more  complex  shapes.  Other  interesting  re- 
sults can  be  found  from  considering  other  spherical  problems  in- 
volving surface  resistance  and  finite  volume  conductivity,  per- 
mittivity, and  permeability  but  wu  do  not  go  into  these  varia- 
tions in  this  note. 

Consider  then  the  problem  of  a plane  wave  incident  on  a 
perfectly  conducting  sphere  as  illustrated  in  figure  Bl.  We 
gave^a  spherical  (r,  0,  $)  coordinate  system  and  unit  vectors 
er,  ee,  which  can  also  be  listed  with  a prime  to  indicate  the 
object  coordinates.  Let  the  sphere  have  radius  a and  let  the 
incident  plane  wave  be  described  as  in  section  2.  The  unit  vec- 
tors for  the  incident  plane  wave  are  illustrated  in  figure  Bl. 

As  in  an  earlier  note*  the  unit  vectors  for  the  plane  wave  are 
expanded  as 


ej  » sin  ( 6^  cos  (♦j)^  + sinfe^  sint^e^  + cos(61)ez 


e2  ■ -cos  (8^)  cos  (^J^  - cos  (0^  sin(^1)ey  + sin(ex)ez 
e3  « •in(41)ex  - cosC^Jey^ 


where  9i  is  the  angle  of  ex  with  respect  to  the  reference  z axis 
(arbitrarily  chosen)  and  4i  is  the  orientation  of  the  projection 
of  ex  on  the  x,  y plane  with  respect  to  the  x axis  (arbitrarily 
chosen)  v The  second  unit  vectoj  B2  is  chosen  in  a plane  parai-  ~ 
lei  to  ex  and  the  s axis  while  eg^is  then  parallel  to  the  x,  y *- 

plane . v J 
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FIGURE  81. 
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The  cartesian  (x,  y,  z)  and  spherical  (r,  0,  0)  coordinate 
systems  ave  related  as 

x = r sir. (6)  cos($) 

y « r sin(0)  sin(4>)  (B2) 

7,  » r cos(0) 

and  similarly  for  primed  coordinates.  Expanding  the  plane  wave 
unit  vectors  in  spherical  coordinates  gives 

e^  ® [cos (0^) cos (0)  + sin(0^) sin(0)coa($-$^) ler 
+ I-cos (9^)  sin(0)  + sinfO^ces^cos^-^)  )e0 

- sin  (0^) sin 

«2  * [sin (0 cos (0)  - cos^)  sin (0) cos ]er 

- [sin(01)  sin(6)  + cos  (0^  cos  (0)  cos  ($-4^)  ]e0  ,B3} 

+ co&(01)  sinW-^e^ 

e3  ■ -sin(O)  sin(^-^1)er 
-cos(0) sin e0 
-cosU-^Je^ 

where  we  can  expand  some  of  the  terms  in  the  forms 
cos  (4-^)  ■ cos (♦j)  cos  (♦)  + sin($3)  sin(t) 

(B4) 

sin($-$^)  - co8(41)sin(4)  - sin($3)cos(t) 

Having  the  direction  of  incidence  and  two  polarizations  ex- 
pressed in  spherical  coordinates  we  can  go  on  to  express  the 
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response  of  the  perfectly  conducting  sphere  to  the  two  delta 
function  plane  waves  3p  as  in  equations  2.17.  This  can  be  fol- 
lowed by  finding  the  delta  function  responses  for  the  surface 
current  density  and  surface  charge  density. 

For  the  incident  delta  function  plane  wave  we  first  need 
spherical  harmonics  and  vector  wave  functions  in  which  to  ex- 
press the  expansion  in  spherical  coordinates.  In  spherical  co-‘ 
ordinates  we  have  the  common  differential  operators  as 

VF  = ®r  If  F + «e  r F+%  r "sTnfST  I?  F 

!r(r2pr)  +F“sIFW  |g-(8in(e)p0) 

* «r[r“Tin ffiT  !s(8in(0) F^)  - F^i^m 
+ e6  [r  sin(9)  Of  Fr  “ r Ir  *rFV] 

+ %[f  lr(rFe)  “?  f©  Fr] 

wherj  ? is  a general  vector  and  F a general  scalar.  Other  oper- 
ator such  as  Laplacian  (V2j?  and  V2F)  can  be  constructed  using 
the  three  in  equations  B5.  These  operators  are  suggestive  of 
ones  that  could  be  defined  to  operate  with  respect  to  the  sur- 
face coordinates  0,  $ on  a unit  sphere.  Of  course  these  operate 
on  scalar  and  vector  quantities  which  are  functions  of  these 
surface  coordinates  (or  considered  as  only  functions  of  such  co- 
ordinates) . Using  a subscript  s to  denote  these  operators  we 
have 

'sF-efl^F  + e*  sin  (S)  Tf  F 


+ r "sfnTtn  I*  P* 


■I—  F 

*0j 


(B5) 


V*  - slnTST  l?(8in(6)F0>  + iTnTST  I?  F* 

V*  -=  8r[ilW  llf<8in<0)  V ’iT5W  h Fe] 

(B6) 

i a „ 

+ ®6  sin(6)  7F  Fr 


" % Iff  Fr 

Note  that  in  the  process  of  removing  all  r factors  and  der iva- 
tives  with  respect  to  r the  Va  operator  leaves  the  units  un- 
changed instead  of  multiplying  the  units  of  F or  F by  meter”  *■. 

Let  us  now  consider  the  spherical  harmonics.  The  scalar 
spherical  harmonics  can  be  written  as 

V e(8,»)  - m-  0.  1,2,  •••,  n (B7) 

n,m,_ 


The  subscript  § meaning  even  or  odd,  indicating  that  one  is  to 
be  chosen  corresponding  to  whether  cos(m$)  or  sin(m^)  respec- 
tively is  intended.  The  Legendre  functions  are  given  their 
standard  definition  such  that  for-l<£<l  the  P"?  (£)  have  the 

definition*^ 


„ , m/2  .m 

p”(C)  s (-i)m(i-e2)  * — p (o 

n d?m  n 


P.K)  S Pn(t)  ' ^r(C2-l)n 

n n -»n_,  j».n 


(B8) 


2 nl  d5 


For  convenience  a subscript  o (for  symmetry)  can  be  used  to  in- 
dicate e or  o or  as  an  index  for  sums  over  both.  Using  hue 
Kronecker  delta  notation  defined  by 


for  oL  - o 2 
for  / o2 


(B9) 


we  can  write  the  spherical  harmonics  as 


Yn,n.o<*'»  » P™(cos(0) ) l«e  acos(ra$)  + 0sin(m*)l  (BIO) 

where  we  use  letters  e,  o for  the  arguments  with  obvious  meaning. 
If  desired  + and  - or  +1  and  -1  could  be  used  to  denote  even  and 
odd  respectively  being  some  kind  of  parity  value. 
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Having  considered  the  scalar  spherical  harmonics  we  now 
need  the  vector  spherical  harmonics , three  kinds  of  them.  In  a 
manner  similar  to  previous  usage?®  we  define  three  types  of  vec- 
tor spherical  harmonics.  This  definition  differs  slightly  from 
one  of  our  previous  notes4  but  our  present  definitions  seem  more 
natural.  The  first  kind  have  only  an  r component  and  are  simply 
defined  as 


n 9 hi  , < 


(Bll) 


( 


The  second  kind  has  only  6 and  $ components  and  is  defined  as 


o 


(6,$) 


V ¥ 
s n,m. 


o 


<8,D 


3 v 

W Y 


n,mro 


(e'*>  + *<fc  slrTTST  5* 


n,m#o 


(6,D 


(B12) 


which  can  be  written  out  as 


$ G(S,<t>)  *eQ 

n,m,0 


dP*™(cos  (6) ) 
n 

M 


(cos  (ra$) 
(sin  (m$) 


pJJ(cosO)) 
sin (6) 


n(n-m+l)  Cl(cos(9)l  (n+1)  (n+m)  Pn-l(cos(e)  * 1(  cos(m^)) 
- . 2n+l  sin  (6)  2 n+1  sin (6)  j j sin(m<J>)l 


+ e.m 
9 


F®(cos(0)) 
sin (0) 


(B13) 


The  third  kind  also  has  only  6 and  $ components  and  is  defined 
as 


L m n(9,D  = V 
n,m,o  s 


(eY  m 
7 n,m,o 


sin  (01 


d 

IT# 


¥ _ ja,$)  - e. 

n,m#0  6 


a 

T5- 


n,m,a 


(e,*>  (Bi4) 


which  can  be  written  out  as 


i 
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<e.*> 


P“(COS(6)) 


n,m,( 


j-sin (»♦) J + dPn(cos(e) ) <co>(w^) 
jcos(m$)  j <16  jsin(m^) 


t m Pn(cos(9))  |-8in(M») 
e0m  sTnX§l  jcos(m^) 


f n(n-m+l)  Pn+l(cos(0))  <n+l)  (n+m)  Pn-l(cos(0) } ] {cos(m*)l 

* l ilnTO + 1 "2n+I »T5n^T J ( s in  (m* ) I 


(B15) 


Some  useful  relations  hold  among  the  three  Kinds  of  vector 
spherical  harmonics  as 


<L  m • e_.  x t rt(e,<|>) 

n,m,o  r n,m,o 


- ”®r  x Sn,m,o(e'*,) 

3n,B,o<6'*>  * VSr  ■ 


(B16) 


e * IV  x ? (e.A) ] 

r 1 s n,mra' 


S « ff<0'*>  " v8  x «<0^) 

n,m.a  s 

These  also  have  various  relations  to  the  scalar  spherical  har- 
monics as 


m O(0'*>  * «rYn  m „<0^> 

n,m,o  r n,nt#o 


” «r  « IV,  x I V„,»,o 


-«r  * Vn.m,0«8<*> 


(B17) 
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Note  that  the  three  types  of  vector  spherical  harmonics  are  mu- 
tually orthogonal  at  each  point  on  a unit  sphere  for  the  same 
set  of  indices.  They  also  are  mutually  orthogonal  in  an  inte- 
gral sense  on  the  unit  sphere  for  any  combination  of  index  sets 
as 

£ £\.m.olt'*'  • * 0 <Bi8> 

£ £\^^  • Sn.,m.,0.(e.*)sin(8)dM8  - 0 

For  the  same  kinds  of  vector  spherical  harmonics  we  have  orth- 
ogonality relationships  on  the  unit  sphere  as 

£ • **•,*•, O' 


[1+  [6o  - 6 IS  J JJL.  4aiSl|  fi  fi  , 6 . 

e,o  oto  ofm  2n+l  (n-m) 1 n,n  m,m'  a, o' 


L L tn.m.o10-*'  * 3n.(n,,j0.(8,8)sin(e)a*d8 


(B19) 
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“ ” + 'se,o-so,o T5TT  TS^ff  {nrn's«,»'4o,o- 

As  a next  st^p-we  need  vector  wave  functions  for  spherical 
coordinates.  Opepart  of  these  functions  comes  from  spherical 
Bessel  functions  which  are  functions  of  Yr  “ ikr  * sr/c.  It  is 
from  these^tnat  the  radial  and  complex  frequency  dependences  are 
formed.  , 'These  are  commonly  expressed  as 
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jn(kr)  , ynl*r> 

hll>  (Kr)  - 3n(kr>  + iynIltrt 

*V  14 


h*25 (Kr)  » “ iyft(kr) 

n 

where  the  jn  ajre  used  for  cas®®d°5  a) 8 ^^used*  for  Outgoing  B 
^.r^uS/i^Kinion  condition  .t  infinity.  For  C * 0 

•.•a  Kbvb19 


in(5>  ” Tin+UVf,l+0<E*>1 

h(J)<5)  - »«>«*>» 

n 

where  the  double  factorial  is  defined  by 

<m(m-2)-**<4)  (2)  for  ra  even 
mU  = * • ' <3)  <1)  for  m odd 


(B21) 


(B22) 


ui  » on  * c—i)  u * i 

«*■  for  present f0n°* 
tions  with  argument  C * H ana  we  ocu. 


i^Uo  = in«c> 


(B23) 


£i2>  5 V° 

The  first  kind  is  used  to  expend  the  IT 

in(Kr) ) and  the  second  to  expand  the  h«(l)  can  be  formed 

as  a linear  combination  of  the  fn  *°*0  the  complex  con- 

these  functions  to  be  real  ^iii  apply  end  *ero»  of  the 

jugate  relationship  of  equation  2^4v  11PP  y thsn 

functions  will  have  complex  conjugate  symmetry.  u. 

constrain  t ♦ 0 the  asymptotic  for©® 
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iiessMa mw 


< 

! 


i 


V‘>  ■ •(T^TnT1140"2'1 

/• 

kn(C>  - c"n_1(2n-l)  I!  Il+OU*)] 
from  which  we  can  make  the  identification 
in«*  * in3nCO  - injn(“U) 

kn(c)  * i’""2hn2) “ i"n"2h^2>  (-ic) 
- -i“nh^2> (C) 


(B24) 


(B25) 


For  these  functions  we  need  a Wronskian  relation  with  re- 
spect to  the  argument  t as 


W{i  U).kw  (c)}  i i (C)k'(O  - i‘(C)k  (C)  - -c“2  (B26) 

n n n n n n 

where  the  prime  with  the  Bessel  function  indicates  differentia- 
tion with  respect  to  the  argument.  Another  related  expression 
is 


in(c)UknU)1'  ‘ k„<«)  lCin(C)  1 ' - -c”1  (B27) 

which  in  Wronskian  form  is 

wUin<0,Ckn(C>)  - -1  (B28) 

These  are  useful  for  simplifying  coefficients  in  the  field, 
current,  and  oharge  expansions. 

These  spherical  Bessel  functions  can  be  written  as  combina- 
tions of  polynomials  and  exponentials  giving  them  a simpler  form 
for  their  exact  representation.  This  is  important  for  finding 
the  poles  for  our  singularity  expansion.  From  a standard  refer- 
ence^ the  spherical  Hankel  functions  can  be  written  for  n « 0, 
1,  2,  •••  as 


(B29) 


n 


•,»>«>  ■ f» W £ f^±4>fT, -i«rB 

B-0 


hu,(E)  . 

$■0 


From  these  we  can  construct  the  other  spherical  Bessel  functions 
as 


Ve>  “ l[k»1,(£)  + hn2l(£)] 

V£>  * -tK1,(£>  ‘ hn2><£>3 

With  C “ i£  we  then  can  write  our  k functions  as 


k (c)  - i“n“V2)  (£) 

n n 


(n+g)  1 
6 1 (n-V) ! 


(2C) 


-e 


(B30) 


(B31) 


with  the  resulting  simplification  of  the  expression  this  last 
form  for  the  spherical  Bessel  functions  for  outgoing  waves  has 
some  useful  advantages.  Next  consider  the  spherical  Bessel 
functions  which  are  analytic  at  C * 0;  these  can  now  be  written 
as 


inu)  - i"jn<c> 


<*>  * •>*«  «>] 


V'  (n+6)  1 , „ ,-B  . , ,.R+1  e~£  (n+B) I ,..,-6 

tza i B-i(n-srr(  2t)  * ( 11  tv  A 67ri“B)T('c) 

p*0  p— 0 


cosh(U^2{(-3)B  + (-i)n"1]_^±|^T(2c)"0'"1 
+ sinh(C)^2t(-l)a  + (-i)n+i]]rin±|M_.{2C)-P”1 


/v . 


(B32) 


\ 
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The  finite  sums  in  equations  B31  and  B32  can  be  expressed  as 
ratios  of  polynomials  in  r For  the  kn  functions  the  numerator 
,*nd  denominator  both  hav<_-  zeros  (in  conjugate  complex  pairs  ex-' 
cept  on  the  Rett)  axis).  For  the  in  functions  we  only  have 
zeros  in  the  complex  t plane;  the  in  are  then  entire  functions. 
All  the  inverse  powers  of  c cancel  when  the  exponentials  are  ex- 
panded as  power  series  in  t. 

Having  the  Bessel  functions  in  the  forms  we  desire  we  now 
consider  the  spherical  vector  wave  functions.  These  are  closely 
related  to  the  spherical  harmonics.  As  a building  block  we  have 
the  spherical  scalar  wave  functions  as 


~U> 


(yr) 


-(f) 

"R,m,0 


(yr 


(B33) 


which  car.  be  written  out  as 


<yr) 


f<£>  (Yr)P™(e,<|>) 
n n 


(cos  (md> ) ) 
| sin(m<j>)  j 


m = 0,  1 , 2,  • * • 


n (B34) 


where  i «=  1,  2 refer  to  in  and  kn  respectively.  Coefficients 
times  this  when  summed  over  all  possible  indices  satisfy  the 
scalar  wave  equation  which  fcr  each  function  we  can  write  in  op- 
erator form  as 

[V2  - y2)3^  (yr)  = 0 (B35) 

II  p III  f \J 


From  the  solution  of  the  scalar  wave  equation  one  con- 
structs as  usual  the  solutions  of  the  vector  wave  equation,  and 
these  are  of  three  kinds.  The  first,  kind  have  zero  curl  but 
non  zero  divergence  and  are  defined  by 


is'1’  jlT?) 

n,m,  j 


= \ vs _(yr) 
y n,m,a 


= [sr  nwr 


f(£i'(yijp  (0,<J>)  + •— — H-l  5 
n w ' ::,m,ov  yr  wn,m,c' 


f (yr) 


(B36) 


v/hich  can  be  written  out  in  components  as 
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n,m. 


dl)  oh 

"<”'o 


«r  (^>  «fr— <»»tOM<») 

yr  <T5  |sin(m+) 


(B37) 


n«”'o 


(Tr)  P£(eo,|9|)  „|-.in(«4)l 
yr  SlnTH  )cos(m$)  | 


The  prime  is  used  to  indicate  derivatives  of  the  functions  with 
respect  to  the  argument  (yr  here) . The  second  kind  have  zero 
divergence  but  non  zero  curl  and  are  defined  by 


a(yr)  = V x 
n,m,o 


: > K 


-r  x VE*1*  (yr) 
n,m,cr 


« -e  x v 5(i)  (yr) 
r sn,m,o  7 


'*Sr  x 


fn  <Yr)l4r  x 3n,m,o(e'*>> 


The  components  are 


M<*>  e (yr)  - 0 

n»M»0 


*e  e 

n,m,0 


(l)  p^tcosje))  hsinM, 


- f^'  (yr) 


cos  (m$) 


(B38) 
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(B39) 


M<£)  <yr) 

♦n,m4 


U)  i.-'dPn(cos(e))  (cos (m$) 


(Yr) 


| sin (m$) 


The  third  kind  also  have  non  zero  curl  but  zero  divergence  and 
are  defined  by 


Kl)m  o(Y?)  = b V x 0(y?) 

n,m,o  y n,mf  a 


7 ° - [v  * N!i.o(6'»)]] 


= -f  v2e<*>  (0,<i>)  + ^(r • V) (e,$) 

Y n,m,a  Y n#iu,o 


+ | (0,<J>) 

Y n,m,o 


■erYrHn;^a(e^)  + 


3 n.U) 


(e,*)  +1  V=i'l  „(0,d») 


TTyrT  v"n,m,o'°'?;  * y v=n,m,0 


-yrf^£)  (yr)  + yrf^£)"(yr)  + 2f^£)  * (yr)J  Pn#m,a  Or*) 


f^£) (yr)l' 


[f„ (yr)  f**' (yr)) 

+ |Yrl  -yr-  J +2  -Tr  “} 


ni)(yr)h  (e 

yr  (®n,m,0 


(^r)  * [yrf^  (yr)]’ 

= n(n+l)-iL_—  ?nfm,a(e,9)  ^•L-^ L 5n,m,a<0'<»  <B40) 


where  we  have  used  the  differential  equation  for  the  spherical 
Bessel  functions  and  some  vector  calculus  identities.lt  The 
components  are 


N<£)  (yr) 
n,m,Q 


f ( ^ tyr) 

n(n+l)  n Pm(oos(e)) |c?9  *♦> 

yr  n |sin(m<j>) 


“e*’  e(yh  . 

n'm'o 


Yrfn;<Yr>,  <(c08(e))(cos(m^) 


sin  (m$) 


(B411) 
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rr 


»***  e(Y?) 

n>m-o 


[yrf^1  <T»>]'  pS(cos<6>)l-«ln<Bt)j 
yr  slnTST  (cos(m$)  j 


Note  that  all  three  kinds  of  vector  wave  functions  satisfy 
the  vector  wave  equation  in  Laplacian  form  which  we  can  summer- 
ize  as 


(v2  - *2lpj  * A 


(B42) 


However  from  the  operator  identity 


Vxyx  = w • - v4 


and  noting  from  their  definitions  that 


(B43) 


(B44) 


VV  • L » i v t V2H]  = yVS  * y2L 


we  can  write  a curl  curl  wave  equation  for  only  the  second  and 
third  kinds  of  vector  wave  functions  as 


[VxVx  + y ] = 0 


(B45) 


The  three  kinds  of  vector  wave  functions  have  some  interrelations 


«,<*?>  ■ -Y?  x I*1*  ( yr) 

nfm,o 


R<‘>  (Yr)  - -i  V X S<*>  (y {) 

n / m / o Y n * in*  o 


„(T?)  . i V x fi<^a<y?> 


n,m,0 


(B46) 
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C.o'-'  --*«[*“  C,o‘*] 


Since  the  vector  spherical  wave  functions  have  their  6,  <>  de- 
pendence expressed  in  terms  of  the  vector  sphericaL  harmonics 
then  these  vector  wave  functions  have  certain  orthogonality 
properties  on  the  unit  sphere  based  on  those  for  the  spherical 
harmonics  (equations  B18  and  B19) . However  since  two  spterical 
vector  harmonics  are  used  for  the  I and  N functions  the  orth- 
ogonality relations  for  the  vector  wave  functions  on  the  unit 
sphere  are  not  as  simple  and  convenient. 


Returning  to  our  dyadir  nlane  wave  from  equation  2..  13  in 
Laplace  form  for  propagation  in  the  direction  Si  as 


X x -ve1*r 
II  = Ie 


(B47) 


where  I is  the  unit  dyadic  which  can  be  expressed  many  ways  such 
as 


t / 1 0 °\ 

I «=  (6h  . ) = ( 0 1 0 ) 

bl'b2  \0  0 1/ 


= exex  + eyey  + 6zez 


~ elel  + e2e2  + e3e3 


erer  + e0e6  + e<pe<p 


(B48) 


This  dyadic  plane  wave  is  expanded  in  our  spherical  vector  wave 
functions  as^O 


oo  n 


!1=EZ  £ (2n+l)  0(9^r Oj) (Yr) 

n=0  m=0  c=e , o 

+rerir[sn,m,o<ei'»i>i!i”i,<,'1fJ>  -5n,m,o,ei'*l>S^m,o<^>]|  <B49> 

where  0i  and  $1  are  angles  giving  the  direction  of  ei  as  used 
previously.  Note  that  for  n = 0 the  summation  is  not  extended 
over  the  M and  ft  functions  which  are  identically  zero. 
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This  dyadic  plane  wave  expansion  is  related  to  integral 
representations  for  the  spherical  vector  wave  functions  for 

i - 1 as20 


\ n+1  rir  r2v  -Ye,*r. 

— JaLe  *»,..,c.(9l'9l,8ln(9l,d9ld9l 


'o  ■'o 


fit  f2v  Ye,*r 
hi  I 8 *n,m,a<9l>»l>8in‘9l>d*ld9l 


'o  •'o 


/iv  / ..n  /.it  /*2n  -Ye,  *r 

C-o^^jT  /D  8 1 


-it  r2n  -Ye,  *r 

/ X 6 5n,m,o(W8in(9l>ddldel 


/*2tt  Ye,*r 

-hf  J 8 m,=  ;91'9l»sin<9l)d*ld9l 

•'o  •'o 


'O  O 


-^ff  "Tei'r3n.n,o«9l'9l>8in<9l>d9ld9l 


-Ye,8^ 


'o  "o 


.n  f2 n ye-  *r 


i r » Ye , r 

= ^X  X 6 ^Vn.a<V*l>8i"(9l>d*ld9I 

Thus  the  spherical  vector  wave  functions  for  Jt  «=  1 can  be  con- 
sidered as  weighted  integrals  over  plane  waves  travelling  all 
possible  directions  of  propagation  ei. 

Having  an  expansion  for  the  dyadic  plane  wave  we  can  find 
the  two  delta  function  plane  waves  for  p «=  2,  3 as 


u = e 8 t. 

P p 1 


(B51) 


-►  -»• 


For  these  we  need  ep  expressed  in  terms  of  eri*  ©01/  «$i  unit 
vectors  which  we  can  see  by  referring  to  figure  B1  are 


-&■  *> 
el  = er. 


e2  “ "e0. 


(B52) 
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Ben  considering  the  teres  in  the  expansion  coe«i=ients  we  have 


c-  * P * ° 

n,m,n 


V8  e(ei^l)B"^I 
n#m»0 


Y e(^ 
n |i#»  _ 


dP^coste^)  (cos(m^)| 
r$!>  * £L^  ”|sin(»*l)j 


v*  e(W 

n,m,Q 


-si57®ir4llt0,<(ei'+1> 

pm(cos(61))  i-sinOn+j)) 
!“  nsin(e1i“"m(co8(m*1)  j 


(B53) 


(0,  -° 


e3-?  ev°l' vl' 
n,m,c 


V3  ete 
n,m,0 


. a \ -JL-  Y 

l'*l)  * nAe0  1 1 

P^tcoste^)  1-sin  (m^)} 
* " -^sinCSp  m(cos(m^1)  j 


e3  * R e 

n»n»»o 


and  for  completeness 


e • f 

e). 


n,m. 


e(0l'*l)5=\ 


dP^tcos^))  | coe  (m<>1)  | 
- -TJ  {sin (»<!>!>! 

Icon  (m^)* 
sin(m$1)| 


(B54) 
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^•5  e(ei'*l'*° 

n,m. 


Vs  Vi"" 

n'm'o 


For  p ■ 2,  3 the  delta  function  plane  waves  (transformed)  can  be 
written  as 


4 4 


s*- v L ♦■t.-.-e..**] 


» n 

21  _ 
n=l  m=Q  o=e,o 

-ye.T  00  n 


(B55) 


»3’S3°  1 -£  £ E [bAt»,oRn*i,0(T*>  -aA,m,oSn^,o(lr?>] 

n®l  m-0  o=e,o 


where 


,,n+l  (211*1)  (n-m)i  „ P^tcoste^l-sintm^) 
i- n (n+1)  7n+mT7  sxnT^T  jcos(m$^) 


n,ra. 


o,m 


(B56) 


r2-a  i (2n+l)  (n-m)  i dP^cos^jcos (m^)] 

o,m*  nTn^TT  (n+m) ! tfl5^  |sin(mt^)j 


n,m, 


The  prime  is  used  with  these  coefficients  to  differentiate  these 
from  the  a0.  Note  that  we  have 


[■-  “Ye1*r-|  ^ -Yej-r 

[e2e  J * e3e 


“Yel*r I * 

le3e  J “ "e; 


-ye^r 


(B57) 


which  is  associated  with  the  curl  relations  between  the  M and  S 
functions.  Furthermore  any  divergenceless  electric  field  expan- 
sion (E)  can  be  converted  to  a magnetic  field  expansion  (H)  by 
dividing  by^the  ^ave  impedance  ? of  the  medium  and  changing 
ig  -N  and  ^ -*■  +M.  To  go  from  H to  E multiply  by  Z and  change 
M -*•  +N  and  N -M. 
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V 


Now  define  two  sets  of  coef f icients  for  q ■ 1,  2 as 


l,n,m,0,p 


2,n,m,o,p 


e'  _ _ for  p - 2 
n,m,o 

b ' for  p - 3 

n,m,o 


b ' for  p » 2 

n,m,  o 

-a'  for  p = 3 

n,m,0  r 


(B58) 


Then  for  p * 2,  3 we  can  write  our  unit  incident  plane  wa”e  as 

“p  " V Yei  C ^ £ 2 K,n,m,o,p?1n!i,,c(',?) 

n=  1 m=0  o~e,o 

+ A2,n,m,0,pNn,m,o 

Our  incident  plane  wave  electric  field  is  written  as 
iinc(r,s)  = EQtf2(s)u2  + f3(s)u3I 


(B59) 


(B60) 


In  the  presence  of  our  perfectly  conducting  sphere  of  radius  a 
we  have  a scattered  electric  field  as 


Kc{*'s) 


Eo[i2(s)5iSC>  + l: 


(B61) 


where  for  p = 2,3  the  scattered  electric  field  response  func- 
tions are 


oo  n 


3‘sc)  = E £ E k(sk  „ A2l  „<t?>  „ J ‘2i  „(yJ)1 

p l i|n^ni|0 n ^ rti  * o 2 ^ o ^ p n ^ in p o J 

n*!  m»0  o*e, o 


(B62) 


Constraining  the  tangential  electric  field  to  be  zero  on  r * a 
requires  for  the  tangential  components 
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* 


1 


i . fxi  n , „ a*1*  (T«er)  +a<«>  „ a<2>  (t.SjI  -s 

r ^ x ^ n $ si  f (#  ip  n § in  i o r x § ix  § hi  |0|p  n i iAf  o x j 


s « fA, a*1’  (yae  i +*<«>  „ a'2!  a(T«ar)1  -« 

r l 2 1 n |Hi| cf i p n,ro,o  r 2|]i|iii|0 iP  n^m^Q  r j 


(B63) 


( 


This  gives  equations  for  the  coefficients  as 


- (sc/  m A » 

l,n,m,o,p  kn(Ya)  l,n,m,o,p 


„c)  _ . 

2,n,m,o,p  [yakn ( ya) J r 2,n,m,o,p 


(B64) 


The  surface  current  and  charge  densities  (equations  2.18, 
2.19,  2.25,  and  2.26)  are  written  as 

E r ~(J  ) ~(J  ) 1 

3g(r\s)  *2^[*2<s)$2  8 (r * ,s)  + f3(s)$3  8 (r',s)J 


(B65) 


+ r~  ~(PS>  + - ~(PS)  + 1 

ps(r',s)  *e0E0[f2(s)U2  8 (r',s)  +f3(s)U3  & (r*,s)J 


r.(ps> 


To  find  the  surface  current  and  charge  densities  we  need  to 
evaluate  just  outside  r * a the  expressions 


^8(r' ,s) 


P8(r\s) 


e^.  x H(r • ,s) 


= eo®r  * S C? * , s) 


(B66) 


The  surface  current  density  response  functions  are  then 

«?•..)  -Si « E t E 

n»l  m=0  o*e,o 

"AlSn'in  cj  o^r2i  0<Yae')  + A,  rt<Yae*) 

X § xi  f m i o i p n / in  i o xr  2 ^ n i in  i c i p n p iti  p 0 xr 


+ *iSV*  a oBn2i  „(va?i)l 

2 1 n 1 m p 0 $ p n p m § 0 xr  J 
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y y j.A 

r &£'oo£'J  Y 


tyai. (ya) 1 ' 


Y« 


in(Ya)  (Yakn(ya)]' 
k(ya)  ya 


-IS  (0  • A * ) 

Jwn,ra,o  * 7 ' 


r Iyain(ya)]’  n± 

+ A2#n#m,a,pl1n(Ya)  ‘ [yak^ya)  ] T kn(Ya)J  ^n,m,or (e ' »*' * 


(B67) 


From  the  Wronskian  relations  for  the  spherical  Bessel  functions 
this  reduces  to 


(3  ) 


- ' ■£  £ oSoh-— <e,-v,7^t 

‘ A2,n,m,o,p^n,m#o *0  * ya[yakn(ya) f ’ | 

The  surface  charge  density  response  functions  are  then 


(ya) 


0*Ps)  (r*  ,s)  - £ / 2 2 |Ai  n m o o®r  * K^m  o(Ya«i) 

p ™“"  ( if  o r n,m,a  r 

n-1  m*0  o=e ,o 

x*  (sc)  ii  . m ( 2 ) (vaeM 

• 1, n,m,o,per  Mn,m,d'Y  er' 


od  n 


+ A2,n,m,o,p?r,Ri“,o<lri,SP 


+ A*sc)  e'  • $(2)  (yae* ) 

2,n,m,d,p  r n,m,d  1 r 


/ 


? ? ? A2,n,m,d,p[n^n+^“ 

n*l  m=0  o=e ,o 


i-(Ya) 


ya 


lyai  (ya)l'  k lyan 

'Ty5irTyiTTr  n(n+1)_ilYir“Jyn,m#a(e,'<*>,)  <B6» 


k (ya) 


(B68) 
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Note  thrt  only  the  q * 2 terms  contribute  to  the  surface  charge 
density.  Using  a Wronskian  relation  we  have 


.<PJ  - 
Up  * <r’,s) 


n**l  ra-0 


Y.  -n(n+l)A2 
o-e,o 


i n #m , o f p 


y (e  * .<!>•) * — = 

n'm'a  ('ra)S[Yakn(Ya)]* 


(B70) 


Note  that  since  n * 1,  2,  2,  •••  there  is  no  pole  at  y*  * 0 in 
either  the  surface  current  density  or  surface  charge  density  re- 
sponse functions. 

Now  that  we  have  explicit  representations  of  the  response 
functions  for  surface  current  and  charge  densities  in  terms  of 
known  functions  we  can  identify  various  terms  with  the  terms  in 
the  singularity  expansion.  Let  us  start  with  the  natural  fre- 
quencies. These  are  the  zeros  defined  by 

’ [sa  t*nK  f)]’  - 0 <B71> 

There  are  then  two  classes  of  natural  frequencies  which  can  be 
labelled  by  q « 1,  2 depending  on  which  of  equations  B71  they 
satisfy.  Clearly  n is  another  index  and  for  each  n there  are 
some  number  of  natural  frequencies  which  we  index  by  n'.  Thus 
the  index  set  a as  applied  to  labelling  the  natural  frequencies 
can  be  written  as  q,  n,  n'  and  we  have 


1 


»n,n' 


0 , 


s = s , 

O q,n,n' 


,n,n’ 


S2 


»n,n' 


0 

(B72) 


Since  from  equation  B31  we  have 


kn(0 


(n+p)  1 
B!  (n-Bil 


(20 


-e 


(B73) 


Then  let  us  write  the  spherical  Bessel  function  terms  in  equa- 
tions B68  and  B65  as 
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where 


t - -2. 
o c 


(B75) 


Since  the  B functions  are  all  ratios  of  polynomials  in  ya  then 
we  can  make  a pole  expansion  of  them.  Note  that  tQ  is  just  the 
turn  on  time  when  the  incident  wave  first  touches  the  sphere  so 
that  the  coupling  coefficients  are  factored  as  in  equation  2.21; 
the  perfectly  conducting  sphere  is  then  an  example  of  this  fac- 
toring. 

The  ratios  of  polynomials  are  written  as 


b2,„«>  = 


B3.n<c» 


. I,  V tn-fB)  1 .-el 

Ch  i 


.(  . V'  (n+BH  ,,.,-B  „„V'„  (n+6) ! ,,.,-B-ll  1 

I Bl  n-6  l<2cl  2t2UBBnri-W<2° 

v B«=0  B=1  ’ 

(BIS) 

B2.n«> 1 

c CC2,n^ 


In  terms  of  the  spherical  Bessel  functions  these  rational  func- 
tions (i.e.  polynomial  ratios)  can  be  written  as 
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Let  us  make  a pole  expansion  of  these  rational  functions  where 
Bl,n  has  the  poles  for  q * 1 and  B2  n an<3  B3tn  both  have  poles 
for  q * 2.  Note  that  Bi,n  has  the  rorm  divided  by  a poly- 
nomial of  degree  n.  B2#n  has  the  form  divided  by  a polynom- 
ial of  degree  n + 1.  B3>n  has  the  form  ^n-1  divided  by  a poly- 
nomial of  degree  n + 1.  Since  n _>  1 then  Bi  n and  B2,n  are  zero 
for  C -*■  0 while  B3  n is  a constant  (for  n = 1)  or  zero  (for 
n 2)  for  C -*■  0.  'Thus  there  are  no  poles  at  s - 0 (consistent 
with  physical  requirements)  . For  s •*  “ all  three  rational  func- 
tions go  to  zero;  thus  there  are  no  poles  at  s * » and  no  con- 
stant terms  in  the  expansions. 

Using  the  rational  functions  C 1 n»  C2,n»  and  C3#n  we  can 
then  write  our  pole  expansions  around  the  sa  = Sq/n/n«  simple 
zeros  of  the  C functions.  Since  the  number  of  zeros  of  a poly- 
nomial is  equal  to  the  degree  of  the  polynomial  then  for  q ■ 1 
we  have  n values  for  sa  and  for  q = 2 we  have  n + 1 values  of  sa. 
Define 


X (b)  = largest  integer  < b (B78) 


Th<»n  have  a rancie  for  our  index  n'  as 


(B77) 
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B2,n(c)  = cTct^TcTP* 
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3'n  r[Ckn(C)]* 


Let  us  make  a pole  expansion  of  these  rational  functions  where 
Bl,n  has  the  poles  for  q = 1 and  B2  n and  B3,n  both  have  poles 
for  q * 2.  Note  that  Bi#n  has  the  form  Cn”l  divided  by  a poly- 
nomial of  degree  n.  B2,n  has  the  form  Cn  divided  by  a polynom- 
ial of  degree  n + 1.  B3#n  has  the  form  Cn“l  divided  by  a poly- 
nomial of  degree  n + 1.  ' Since  n > 1 then  Bi  n and  B2,n  a*®  zero 
for  c -*■  0 while  B3  n is  a constant  (for  n = 1)  or  zero  (for 
n > 2)  for  z,  -*■  0.  ’Thus  there  are  no  poles  at  s = 0 (consistent 
with  physical  requirements) . For  s -*■  “»  all  three  rational  func- 
tions go  to  zero;  thus  there  are  no  poles  at  s » » and  no  con- 
stant terms  in  the  expansions. 

Using  the  rational  functions  Cl  n,  C2.n#  and  C3,n  we  can 
then  write  our  pole  expansions  around  the  s0  » Sqfn,n»  simple 
zeros  of  the  C functions.  Since  the  number  of  zeros  of  a poly- 
nomial is  equal  to  the  degree  of  the  polynomial  then  for  q *=  1 
we  have  n values  for  sa  and  for  q * 2 we  have  n + 1 values  of  sa. 
Define 


X (b)  = largest  integer  <_  b 


(B78) 


Then  we  have  a range  for  our  index  n’  as 

for  <1  rJL 

-X (j)  < n'  < X with  n'  ? 0 if  n is  even 

(B79) 

for  q « 2 

-X(2y4  £ *»'  < with  n’  ft  0 if  n is  odd 

For  the  sa  we  then  automatically  have  the  convenient  relation 
possible  between  n*  arid  -n ' indices  as 


« 
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sq,n,-n'  " Sq,n,n* 


Im[s..  n]  = 0 for  n odd 
x / n / u 


Im[s0  nj  » 0 for  n even 

Z fU  ,\J 


It  is  known*5*  that  the  zeros  of  kn  lie  approximately  on  an  arc 
in  the  left  half  of  the  s plane  joining  sa/c  = -in  to  sa/c  = +in 
and  passing  through  sa/c  = -.66n.  The  zeros  of  tCkn(C)]'  behave 
similarly.  Then  a convenient  way  to  identify  the  sa  with  spe- 
cific n'  is  to  start  with  the  most  negative  value  of  n'  from 
equations  B79  and  assign  it  to  the  sa  with  the  most  negative 
Im[sa]  and  progressively  work  up  to  the  most  positive  Im[sa] . 


B 


Thus  our  pole  expansions  may  be  written  as 
»©  B 
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(B81) 


»(=£) 


B3,n(cJ  ^n+1  s_s2  n n1 

n-.-A(2±i)  2'n’n 


2tn'n‘  with  n’  f 0 if  n is  c*!c 


where 


D1 ,n,n’  ~ afe  C,n,l(i:)] 


-1 


1 »* 


l,n,-n' 


C=s 


l,n,n*c 


“ 1 

“i.n.n'-sLaTS,^5’] 


2 ,nr-n' 


C-* 


2,n,n'c 


96 


it  jMMl itM  Jll  ,1^1. 


(B82) 


( 


= etc 
, n , n ' ~ a [c 


-,-1 

r ^ ,-l 

[dC  Cn,3(?)] 

a”f[c  3c  SvJ 

3 « 

5“s2,n,n’c  f *' 

2 , n , n ' c 


-°3,n,-n' “Sj  n n.a  D2,n,n- 

To  see  that  the  poles  must  be  simple  poles  note  that  the  zeros 
of  the  C functions  are  all  simple  zeros  because  they  are  the 
zeros  of  kn(C)  and  tCkn(C)J'.  Observe  the  differential  equation 
for  the  spherical  Bessel  functions  is1^ 

C2f*A)"<C>  + 2cf nA)'u)  ‘ [C2+n(n+l)]f^)  (?)  = 0 (B83) 


Suppose  C has  a zero  at  ?0  ^ 0.  Then  since  fn  is  analytic  at 
this  zero  we  can  write  a convergent  power  series  expansion  in  a 
neighborhood  cf  C = Co*  If  the  zero  is  higher  than  first  order, 
say  (C  - Co)^*  then  both  fn  and  fn  are  zero  at  Co»  but  this 
forces  fn  to  also  be  zero  at  Co  so  the  zero  had  to  be  at  least 
(C  “ Co) ^ as  the  leading  term  in  the  power  series.  Then  divide 
through  by  C - Co»  but  fn  and  fn  are  still  zero  making  fn  still 
zero.  This  process  continues  to  make  all  terms  in  the  power 
( series  zero  and  the  function  then  identically  zero.  Thus  the 

zeros  are  all  simple  for  C ^ 0.  Similarly  Ckn{C)  satisfies  the 
Riccati-Bessel  equationl9 


-y5 [cf^A)  U)l"  - Cf(A)U)  =0  (B84) 

C+n(n+l)  L n J “ 


Differentiating  gives  a differential  equation  for  Ickn(c)l'  as 

.2 


C +n(n+l) 


KM ««] 


, n (n+1) 

+ ~J 


IC  +n(n+l) ] 


(B85) 


I (OA" 


Clearing  the  denominators  this  equation  has  the  same  form  as 
equation  B83  and  so  for  C ¥■  0 all  zeros  must  be  simple.  There- 
fore the  perfectly  conducting  sphere  has  only  simple  poles  in 
it.o  surface  current  density  and  surface  charge  density  response 
functions. 


To  see  some  of  the  numbers  we  can  write  out  the  first  few 
terms,  say  for  n ■ 1,  from 


( 
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C-j  , - -2c  - 1 


From  these  results  we  can  construct  a table  allowing  q « 3 for 
the  surface  charge  density  as 


*q,l,n'c 


“ D i • 

c q,l»n 


(like  q=2) 


4 + i £ 
2 1 2 

_I  + i H 

2 1 2 

(like  q=2) 


7 + 1 T 

-i  ^ 

-i  3 


Table  Bl.  Pole  expansion  terms  for  n «*  1 

If  one  wishes  these  natural  frequencies  and  D coefficients  can 
be  generated  to  obtain  any  number  of  terms  in  the  expansions. 

For  n = 1,  2,  3 the  zeros  of  the  Cg,n  polynomials  can  be 
found  from  formulas  for  the  zeros  of  up  to  quartic  polynomials. 19 
There  is  disagreement  in  a few  cases  of  the  natural  frequencies 
with  the  numbers  in  Stratton. However,  the  present  results 
appear  to  be  more  accurate  and  are  confirmed  by  Dr.  Marin  (pri- 
vate communication) . The  zeros  have  been  substituted  in  the 
polynomials  to  check  that  in  fact  the  results  are  closest  to  the 
true  zeros  to  the  number  of  places  listed.  Figure  B2  shows  the 
positions  in  the  complex  sa/c  plane  (normalized  s plane) . The 
division  of  the  natural  frequencies  into  q = 1 and  q **  2 varie- 
ties has  a physical  basis  in  that  only  the  q = 2 poles  contrib- 
ute to  the  surface  charge  density.  Another  way  to  view  this  is 
from  a property  of  the  vector  spherical  harmonics  as 


7S  x ■ ve  x IvsYn>m<0ce,4)]  = J ' 

Vs  • * 7.  • <’s  x tSrYn,m,o(e'*>>>  ’ 0 
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FIGURE  82.  NATURAL  FREQUENCIES  OF  THE  PERFECTLY 
CONDUCTING  SPHERE  FOR  USE  WITH  EXTERIOR 
INCIDENT  WAVE 
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Natural  frequencies  for  n = 1,  2,  3 


from  which  we  can  divide  the  surface  current  density  into  solen- 
oidal  terms  (precisely  q = 1)  and  irrotational  terms  (precisely 
q * 2) . (See  ref.  6 for  more  elaboration  of  this  point  in  the 
general  case.)  However,  the  n,  n*  division  may  not  give  the 
best  indexing.  Referring  to  figure  B2  there  are  various  possi- 
ble paths  through  the  complex  s plane  which  one  might  trace  to 
connect  poles  with  the  same  q index.  For  a given  n (and  q) 
vhere  are  many  modal  distributions  generated  by  varying  m and  a, 
all  applying  to  the  entire  set  of  natural  frequencies  generated 
by  varying  n'.  This  is  a very  degenerate  situation  in  both  nat- 
ural frequencies  and  modes.  Perhaps  more  insight  into  the  divi- 
sion of  the  indices  for  the  natural  frequencies  and  modes  can  be 
gained  from  a group  theory  investigation  of  the  symmetry  proper- 
ties. Symmetry  planes  and  axes  can  be  used  to  divide  up  natural 
modes,  and  thereby  natural  frequencies  as  well.  A diagram  as  in 
figure  B2  is  useful  in  that  it  can  suggest  ways  of  grouping  nat- 
ural frequencies,  even  for  objects  more  complex  than  a perfectly 
conducting  sphere.  Note  that  the  pattern  of  the  natural  fre- 
quencies tends  to  fill  up  the  left  half  of  the  s plane.  This 
two  dimensional  pole  distribution  may  be  associated  with  the 
distributed  nature  of  the  body;  we  are  dealing  with  surface  cur- 
rent and  charge  densities.  For  cases  that  the  currents  are 
idealized  as  on  one  dimensional  paths  then  the  pole  distribution 
should  be  much  less  dense  and  localized  to  "discrete  paths”  in 
the  complex  s plane;  in  any  event  there  would  be  one  less  index 
needed  and  not  say  n and  m both  for  the  modes  or  perhaps  not  n 
and  n'  both  for  the  frequencies. 

Next  we  consider  the  natural  modes  of  the  sphere.  This  is 
the  part  of  the  singularity  expansion  where  the  object  coordi- 
nates are  expressed.  For  the  surface  current  density  the  natu- 
ral modes  are  readily  identified  in  equation  B68  as 
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(B88) 


Furthermore  the  surface  charge  density  natural  modes  are  readily 
identified  in  equation  B69  as 

<PJ  - *Pj 

v s (r')  = v (r') 

a ' ' q , n , m , a ' 


for  q = 1 


Yn,m,0<6''*'>  ,or  <!  ' 2 


(389) 


Thus  for  the  surface  charge  density  we  can  drop  the  q index  in 
the  summation  understanding  that  only  q = 2 is  used  fcr  the  nat- 
ural frequencies  and  modes.  Back  in  equation  2.27  we  observed  a 
relation  between  the  natural  modes  for  the  surface  current  den- 
sity and  surface  charge  density  as 
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(B90) 


where  aa  is  an  arbitrary  constant  depending  on  how  one  has  de- 
fined the  natural  modes  since  the  modes  can  be  multiplied  by  any 
non  zero  complex  constant.  For  the  sphere  problem  we  can  write 
this  in  terms  of  the  divergence  on  the  unit  sphere  as 
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Now  we  have 
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so  that  only  q = 2 natural  inodes  have  non  zero  surface  charge 
density.  From  the  fact  that  the  E functions  defined  in  equation 
B34  satisfy  the  scalar  wave  equation  as  in  equation  B35,  then 
noting  that  from  the  separation  equation  for  the  radial  functions 
in  spherical  coordinates  as 

7 art1 2 If  £n*’  <*r>)  - [v2  + <*r>  - 0 (B93) 

equation  B35  can  then  be  written  as 
0 ■ 


[72-Y2)^M(Yr)Yn<n(0(6,*) 


"'"r-  + f7'(Yr)V-Yn,m,a(e'*' 
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from  which  we  find 


V'  • 3 (0* ,4')  = V,2y  (61 ,*') 

s n,m,a  ' s n,m,ov 


- -n  (n+1)  Y „(9',<fr*) 
n,m,a 


(B95) 


Thus  we  have 


a„  = a^  = n (n+1) a 
<x  n 


(B9f) 


where  only  q = 2 is  relevant. 

Having  the  modes  we  now  only  need  the  coupling  coefficients, 
these  are  everything  that  remains  in  equations  B68  and  B69. 

Since  the  poles  are  all  simple  the  basic  equations  we  are  match- 
ing come  from  equations  2.31  giving  the  plane  wave  response 
functions  as 
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The  coupling  coefficients  are  then  for  q = 1,  2 

C | „ _ / ^-i)  (~*1)  A y\  m rt  Y\  r»  * 
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for  q = 1,  p = 2 
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for  q = 2,  p = 2 
for  q = 2,  p **  3 


(B98) 


The  D coefficients  are  evaluated  from  equations  B82  and  a few 
are  listed  in  table  Bl.  The  and  coefficients  are  found 
explicitly  in  equations  B56.  Note  that  the  surface  charge  den- 
sity expansion  U3es  only  q - 2 for  the  ea  and  the  D<*.  While  we 
can  calculate  D3>n  n»  as  well  it  is  simply  D2,»,n'  c/(s2,n,n,a) ; 
using  the  results  for  aa  with  this  and  the  n(n  + 1)  coefficient 
in  equation  B70  for  the  surface  charge  density  one  can  see  that 
the  same  answer  for  the  surface  charge  density  expansion  results, 
thereby  giving  a check. 

Equations  B97  then  explicitly  give  the  singularity  expan- 
sion for  the  surface  current  and  charge  densities.  The  natural 
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frequencies  are  the  zeros  of  the  Ca  in  equation  B76;  the  natural 
inodes  are  ir  equations  B88  and  B89;  the  aa  are  given  by  equation 
B96;  the  turn-on  time  t0  is  given  in  equation  B75;  the  coupling 
coefficients  are  given  by  equation  B98  together  with  equations 
B82  and  B56. 


Now  that  we  have  the  singularity  expansions  for  the  delta 
function  response  functions  we  can  consider  arbitrary  waveforms 
by  taking  their  Laplace  transforms  and  splitting  the  response 
into  a part  associated  with  the  singularities  of  the  perfectly 
conducting  sphere  and  a part  associated  with  the  waveform  singu- 
larities as  developed  in  section  II.  For  convenience  let  the 
incident  wave  be  a step  function.  Then  from  equations  2.43  the 
surface  current  density  expansion  is 
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and  from  equations  2.50  the  surface  charge  density  step  response 
is  written  as 
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The  static  surface  current  density  response  is 
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where  for  m ® 0,  1 we  have 


4VSl,«,o  IorP'2 
1’1’B’*,P  (tV31,b>,o  for  p * 3 


TlSl'  ®p>  •3l,m,o<ei’*l> 


(B102) 


from  which  the  dyadic  surface  current  density  static  response 
function  can  be  written  as 
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Ijfith  the  direction  of  the  static  magnetic  field  (which  is  ei  x 
ep  for  our  plane  wave  problem)  dot  multiplied  on  the  right  we 
obtain  the  static  surface  current  density  response.  The  static 
surface  charge  density  is 
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where  for  m « 0,  1 we  have 
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from  which  the  vector  surface  charge  density  response  function 
can  be  written  as 
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Dot  multiplying  this  by  ep  (the  direction  of  the  static  electric 
field  in  our  plane  wave  problem)  gives  the  static  charge  density 
response.  Thus  for  the  perfectly  conducting  sphere  the  plane 
wave  delta  function  response  can  be  converted  to  the  plane  wave 
step  function  response  of  both  surface  current  and  charge  densi- 
ties by  multiplying  each  term  by  l/sqfn,n*  and  adding  a static 
term  (consisting  of  a few  simple  known  functions)  with  a unit 
step  function  turning  on  at  time  to. 

This  basically  completes  the  singularity  expansion  of  the 
response  of  a perfectly  conducting  sphere  to  an  incident  plane 
wave  for  simple  waveforms.  As  long  as  the  incident  waveform  can 
be  expressed  only  in  terms  of  poles  the  delta  function  response 
as  in  equations  B97  can  be  combined  with  the  waveform  poles  and 
through  partial  fraction  expansion  as  discussed  in  section  2 
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the  response  can  be  split  into  a waveform  part  and  an  object 
part.  Also  as  discussed  in  section  II  the  time  domain  response 
is  easily  and  directly  obtained  from  the  pole  expansion  since 
the  frequency  dependence  of  the  coupling  coefficients  factors 
out  as  a common  delay  term  e“sto.  For  the  case  of  the  step 
function  response  the  Laplace  form  as  in  equations  B99  and  B100 
can  be  immediately  converted  to  time  domain  using  equations  2.64 
and  2.65  respectively  since  the  poles  are  all  simple  poles. 

Note  also  that  the  natural  mode  functions  are  real  and  the  a0 
are  also  real;  both  of  these  can  then  be  moved  out  (together 
with  c)  from  the  Re  function  in  equations  2.64  and  2.65  leaving 
only  the  natural  frequencies,  coupling  coefficients  and  oscilla- 
tory exponentials  for  n'  > 0 as  the  only  complex  terms  inside 
the  Re  function.  Of  course  the  ca  coupling  coefficients  can  be 
written  in  the  form  an,m,oDct  or  bn,m,aDa  and  only  the  Da  coef- 
ficients are  complex. 

As  an  example  of  the  time  domain  response  consider  just  a 
few  terms  in  the  step  function  response,  say  for  p = 2,  n = 1, 

$1  = 0,  and  6i  = tr/2  so  that  we  have  a vertically  polarized  wave 
propagating  parallel  tc  the  x axis.  Then  we  have 
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where  the  last  terms  indicate  the  remainder  terius.  Writing 
these  out  we  have 
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These  can  be  summarized  as 
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For  reference  we  have 
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For  comparison  to  these  results  for  the  step  response  of  a 
sphere  one  can  consider  the  numerical  results  for  the  step  re- 
sponse graphed  in  another  note. 21  Consider  the  case  in  that 
note  that  the  perfectly  conducting  plane  is  infinitely  far  away 
from  the  perfectly  conducting  sphere.  Note  that  the  basic  ring- 
ing period  agrees  closely  with  (in  units  of  ct/a)  and  that 

in  one  period  the  amplitude  of  the  ringing  decays  by  approxi- 
mately and  that  even  the  coefficients  of  the  ringing 

terms  in  equations  B109  and  B110  give  about  the  correct  ampli- 
tudes for  the  oscillations.  In  the  referenced  note  only  the 
total  current  crossing  the  equator  0'  = ir/2  is  considered,  and 
so  the  comparison  has  to  neglect  the  first  term  in  equation  B1Q9 
which  gives  no  contribution  in  this  case.  Note  that  for  times 
-a/c  < t < a/c  we  do  not  expect  the  first  terms  only  to  accu- 
rately describe  the  surface  current  and  charge  densities  because 
the  terms  for  n > 1 have  not  had  a chance  to  decay  to  zero  ampli- 
tude. Even  so,  the  first  few  terms  give  a simple  description 
with  some  features  of  the  surface  current  and  charge  densities 
even  at  such  early  times.  As  time  goes  on  the  first  few  terms 
asymptotically  give  the  exact  results. 
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